Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,664,138 A→C 53.1% E130D (GAA→GAC)  hcaF → 3‑phenylpropionate dioxygenase, small (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,664,1380AC53.1% 104.5 / 207.7 162E130D (GAA→GAC) hcaF3‑phenylpropionate dioxygenase, small (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (34/42);  new base C (48/38);  total (82/80)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 165 total aligned reads displayed.

GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC  >  CP009273/2663992‑2664283
                                                                                                                                                  |                                                                                                                                                 
gAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGaaa                                                                                                                                                  <  1:2876566/148‑1 (MQ=255)
 aaaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGacac                                                                                                                                                 <  1:1993404/148‑1 (MQ=255)
 aaaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAAc                                                                                                                                                 <  2:2203629/148‑1 (MQ=255)
  aaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACa                                                                                                                                                >  2:3217254/1‑148 (MQ=255)
   aCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAt                                                                                                                                               >  1:2929910/1‑148 (MQ=255)
     gggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTc                                                                                                                                             <  2:2581265/148‑1 (MQ=255)
        cATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTAt                                                                                                                                          <  1:791586/148‑1 (MQ=255)
           ggCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                       <  2:266773/148‑1 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAACAGAGCGCGATGACACATTCTCTGTTGGa                                                                                                                                    >  2:1518136/1‑148 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGa                                                                                                                                    <  1:1411570/148‑1 (MQ=255)
                gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAc                                                                                                                                  >  1:2106095/1‑148 (MQ=255)
                 ggCAGAAGAGCCGCCGTCACGCACTCGTCACTTAATCAGCAACTGCCAGGTAAGCGAAACCGACATCCCAAATGTATTTGCCGTGAGGGTGAATTATCTACTTTATCGGGCACAAAAAGAGCGCGATGAAACATTATATGTTGGAAcg                                                                                                                                 >  2:2940852/1‑148 (MQ=255)
                    agaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGt                                                                                                                              <  2:1583100/148‑1 (MQ=255)
                    agaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGt                                                                                                                              <  2:592771/148‑1 (MQ=255)
                     gaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGtt                                                                                                                             >  1:3025796/1‑148 (MQ=255)
                     gaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTCTGTTGGAACGCGtt                                                                                                                             >  2:655483/1‑148 (MQ=255)
                      aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                                            >  1:2218116/1‑148 (MQ=255)
                      aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                                            <  2:705540/148‑1 (MQ=255)
                            cgccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACaa                                                                                                                      <  2:1519340/148‑1 (MQ=255)
                              ccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAg                                                                                                                    <  2:1232522/148‑1 (MQ=255)
                               cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGt                                                                                                                   <  2:1263427/148‑1 (MQ=255)
                               cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGt                                                                                                                   <  2:2912104/148‑1 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTGAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                >  1:2453622/1‑148 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                >  1:3118492/1‑148 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                >  1:2354562/1‑148 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                >  1:3118415/1‑148 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                <  2:769558/148‑1 (MQ=255)
                                    cGCACCCGTCACTTCATCAGCAACTGCCAGAGAACCGAAACCGACATCCCAACCGTATGTGCTGGGCGGGTAAATTATCTGCTTTATCGGGCACAACCAGAGCGCGATGACACAGTCTATGTTGGAACGCGTTTCGACAAAGTTCGcc                                                                                                              <  2:2136455/148‑1 (MQ=255)
                                    cGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGcc                                                                                                              >  1:1963986/1‑148 (MQ=255)
                                       aCCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTc                                                                                                           <  1:1036108/148‑1 (MQ=255)
                                        cccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                                          <  2:2729395/148‑1 (MQ=255)
                                          cGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTgg                                                                                                        <  2:2574335/148‑1 (MQ=255)
                                             cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAg                                                                                                     <  2:3170431/148‑1 (MQ=255)
                                                ttAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATg                                                                                                  <  1:2671001/148‑1 (MQ=255)
                                                 tAATCAGCAACTGCCAGATAAGCGACACCGACAGCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                                                 <  2:913665/148‑1 (MQ=255)
                                                 tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                                                 >  2:2632290/1‑148 (MQ=255)
                                                 tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                                                 <  2:920414/148‑1 (MQ=255)
                                                     cAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAc                                                                                             <  1:759884/148‑1 (MQ=255)
                                                     cAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAc                                                                                             >  1:1426265/1‑148 (MQ=255)
                                                        cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                                          >  1:2557316/1‑148 (MQ=255)
                                                        cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                                          >  2:2212224/1‑148 (MQ=255)
                                                               cAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGc                                                                                   <  1:1578315/148‑1 (MQ=255)
                                                               cAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGc                                                                                   <  2:1013795/148‑1 (MQ=255)
                                                                aGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCt                                                                                  <  2:3145407/148‑1 (MQ=255)
                                                                    aaGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGaa                                                                              >  2:714804/1‑148 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACg                                                                            >  2:1763109/1‑148 (MQ=255)
                                                                       cGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                           >  1:508641/1‑148 (MQ=255)
                                                                       cGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                           >  1:1032717/1‑148 (MQ=255)
                                                                         aaaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                         >  2:7949/1‑148 (MQ=255)
                                                                          aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGat                                                                        <  2:1217305/148‑1 (MQ=255)
                                                                          aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGat                                                                        >  1:2859627/1‑148 (MQ=255)
                                                                           aCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGata                                                                       >  1:2569425/1‑148 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGt                                                                   <  1:3255960/148‑1 (MQ=255)
                                                                                cATCCCACCCGTATTTGCTGTGCGGGTAACTTATCTGCGTTATCGGGCACCAAAACAGCGCGATGACTCATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                  <  2:15746/148‑1 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACCAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                  <  2:1692858/148‑1 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                  <  2:1425470/148‑1 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                  <  2:2254950/148‑1 (MQ=255)
                                                                                 aTCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                 <  2:60403/148‑1 (MQ=255)
                                                                                    ccAAATGTATTTGCCGTGAGGGTGAATTATCTACTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTgg                                                              <  2:2064619/148‑1 (MQ=255)
                                                                                      aaaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAt                                                            <  1:2471700/148‑1 (MQ=255)
                                                                                       aaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                                           >  2:3313123/1‑148 (MQ=255)
                                                                                            aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCg                                                      >  1:2633487/1‑148 (MQ=255)
                                                                                              ttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGt                                                    >  2:349207/1‑148 (MQ=255)
                                                                                                   gtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCa                                               >  2:483754/1‑148 (MQ=255)
                                                                                                     gCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACt                                             <  1:2826265/148‑1 (MQ=255)
                                                                                                        ggTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTcc                                          >  2:1501417/1‑148 (MQ=255)
                                                                                                         gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTccc                                         >  2:3197290/1‑148 (MQ=255)
                                                                                                         gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTccc                                         <  2:77770/148‑1 (MQ=255)
                                                                                                             aTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                     <  2:1426265/148‑1 (MQ=255)
                                                                                                              ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                    >  2:847387/1‑148 (MQ=255)
                                                                                                              ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                    >  1:2513438/1‑148 (MQ=255)
                                                                                                              ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                    <  1:472567/148‑1 (MQ=255)
                                                                                                              ttATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                    <  2:935825/148‑1 (MQ=255)
                                                                                                                 tCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTg                                 >  2:1795380/1‑148 (MQ=255)
                                                                                                                  cTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGa                                <  2:1390491/148‑1 (MQ=255)
                                                                                                                  cTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGa                                >  1:876628/1‑148 (MQ=255)
                                                                                                                     cTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                             >  1:2774294/1‑148 (MQ=255)
                                                                                                                     cTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                             >  1:1702673/1‑148 (MQ=255)
                                                                                                                         aTCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt                         >  1:1856474/1‑148 (MQ=255)
                                                                                                                         aTCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt                         >  1:1046828/1‑148 (MQ=255)
                                                                                                                          tCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                        >  1:2166897/1‑148 (MQ=255)
                                                                                                                          tCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                        <  1:390041/148‑1 (MQ=255)
                                                                                                                            gggCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGtt                      >  2:371922/1‑148 (MQ=255)
                                                                                                                             ggCACAAAAAGAGCGCGATGAAACATTCTATGTTGTAACGCGTTTCGAAAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATAACTTCCCAGAACCTGAGTGTAATGTTc                     >  1:587317/1‑148 (MQ=255)
                                                                                                                               cacaAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtg                   <  2:2154707/148‑1 (MQ=255)
                                                                                                                               cacaAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtg                   <  2:1101551/148‑1 (MQ=255)
                                                                                                                                 caAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgat                 >  1:1207953/1‑148 (MQ=255)
                                                                                                                                   aaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga               <  2:163441/148‑1 (MQ=255)
                                                                                                                                     aaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAt             <  2:2106095/148‑1 (MQ=255)
                                                                                                                                        agCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa          >  1:874012/1‑148 (MQ=255)
                                                                                                                                        agCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa          >  1:416435/1‑148 (MQ=255)
                                                                                                                                        agCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa          >  2:44268/1‑148 (MQ=255)
                                                                                                                                         gcgcGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGaa         <  2:2180513/148‑1 (MQ=255)
                                                                                                                                          cgcgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAt        >  2:250317/1‑148 (MQ=255)
                                                                                                                                          cgcgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTAATGTTATGATGAATCGAAt        >  1:1721216/1‑148 (MQ=255)
                                                                                                                                            cgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAttt      <  2:460743/148‑1 (MQ=255)
                                                                                                                                            cgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAttt      >  2:1086394/1‑148 (MQ=255)
                                                                                                                                             gATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCAGATGAATCGAATTTa     >  1:1091818/1‑148 (MQ=255)
                                                                                                                                              aTGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCGTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAt    <  1:655483/148‑1 (MQ=255)
                                                                                                                                                gAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc  >  1:259175/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC  >  CP009273/2663992‑2664283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTG  >  CP009273/2663992‑2664286
                                                                                                                                                  |                                                                                                                                                    
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACA                                                                                                                                                   <  SRR19973070.2895198/150‑1 (MQ=60)
 AAACGGTCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACACGTCAATTAATCAGAAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGAGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACAT                                                                                                                                                  <  SRR19973070.46270/150‑1 (MQ=60)
 AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACAT                                                                                                                                                  <  SRR19973070.2006903/150‑1 (MQ=60)
 AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAT                                                                                                                                                  >  SRR19973070.2948815/1‑150 (MQ=60)
 AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAT                                                                                                                                                  >  SRR19973070.2950054/1‑150 (MQ=60)
        CATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGT                                                                                                                                           <  SRR19973070.797068/150‑1 (MQ=60)
              CTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAC                                                                                                                                     <  SRR19973070.1421205/150‑1 (MQ=60)
              CTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAC                                                                                                                                     >  SRR19973070.1500873/1‑150 (MQ=60)
              CTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAC                                                                                                                                     >  SRR19973070.2120222/1‑150 (MQ=60)
              CTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAC                                                                                                                                     <  SRR19973070.2442984/150‑1 (MQ=60)
                   CAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTT                                                                                                                                <  SRR19973070.306262/150‑1 (MQ=60)
                   CAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTT                                                                                                                                >  SRR19973070.3045070/1‑150 (MQ=60)
                    AGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTT                                                                                                                               >  SRR19973070.2232688/1‑150 (MQ=60)
                      AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                             <  SRR19973070.1763816/150‑1 (MQ=60)
                      AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCG                                                                                                                             <  SRR19973070.2492525/150‑1 (MQ=60)
                                GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTGAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                   >  SRR19973070.2469667/1‑150 (MQ=60)
                                GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                   >  SRR19973070.2370158/1‑150 (MQ=60)
                                GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                   >  SRR19973070.3138614/1‑150 (MQ=60)
                                GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                   >  SRR19973070.2370044/1‑150 (MQ=60)
                                GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                   >  SRR19973070.3138537/1‑150 (MQ=60)
                                  CACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCC                                                                                                                 >  SRR19973070.1977239/1‑150 (MQ=60)
                                       ACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTG                                                                                                            >  SRR19973070.1391069/1‑150 (MQ=60)
                                       ACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTG                                                                                                            <  SRR19973070.1044122/150‑1 (MQ=60)
                                       ACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTG                                                                                                            <  SRR19973070.1043429/150‑1 (MQ=60)
                                                TTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGAC                                                                                                   <  SRR19973070.2688351/150‑1 (MQ=60)
                                                 TAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACA                                                                                                  <  SRR19973070.2160186/150‑1 (MQ=60)
                                                   ATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAC                                                                                                >  SRR19973070.1436089/1‑150 (MQ=60)
                                                     CAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTG                                                                                              <  SRR19973070.765226/150‑1 (MQ=60)
                                                      AGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGG                                                                                             >  SRR19973070.2574024/1‑150 (MQ=60)
                                                               CAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTG                                                                                    <  SRR19973070.1589384/150‑1 (MQ=60)
                                                                AGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGG                                                                                   >  SRR19973070.2895085/1‑150 (MQ=60)
                                                                AGAAAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTAATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGATCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGG                                                                                   >  SRR19973070.381558/1‑150 (MQ=60)
                                                                     AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                              >  SRR19973070.1040022/1‑150 (MQ=60)
                                                                     AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                              >  SRR19973070.512152/1‑150 (MQ=60)
                                                                       CGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGA                                                                            <  SRR19973070.868249/150‑1 (MQ=60)
                                                                        GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTTCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                           >  SRR19973070.1928570/1‑150 (MQ=60)
                                                                        GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                           >  SRR19973070.2878159/1‑150 (MQ=60)
                                                                         AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                          >  SRR19973070.2586215/1‑150 (MQ=60)
                                                                         AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                          >  SRR19973070.1673584/1‑150 (MQ=60)
                                                                               ACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCC                                                                    <  SRR19973070.3276932/150‑1 (MQ=60)
                                                                                     CAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATC                                                              <  SRR19973070.3261550/150‑1 (MQ=60)
                                                                                      AAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCA                                                             <  SRR19973070.2487808/150‑1 (MQ=60)
                                                                                          GTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCG                                                         >  SRR19973070.2650658/1‑150 (MQ=60)
                                                                                            ATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGT                                                       <  SRR19973070.67664/150‑1 (MQ=60)
                                                                                                     GCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTC                                              <  SRR19973070.2844570/150‑1 (MQ=60)
                                                                                                     GCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTC                                              <  SRR19973070.2844555/150‑1 (MQ=60)
                                                                                                        GGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCA                                           >  SRR19973070.2534316/1‑150 (MQ=60)
                                                                                                            AATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAC                                       >  SRR19973070.2529812/1‑150 (MQ=60)
                                                                                                              TTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCT                                     <  SRR19973070.475910/150‑1 (MQ=60)
                                                                                                                ATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGA                                   >  SRR19973070.882791/1‑150 (MQ=60)
                                                                                                                   TGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTG                                >  SRR19973070.2792086/1‑150 (MQ=60)
                                                                                                                   TGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTG                                >  SRR19973070.1714412/1‑150 (MQ=60)
                                                                                                                       TTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACT                            >  SRR19973070.1869188/1‑150 (MQ=60)
                                                                                                                       TTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACT                            >  SRR19973070.1054203/1‑150 (MQ=60)
                                                                                                                        TATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTG                           >  SRR19973070.2181275/1‑150 (MQ=60)
                                                                                                                          TCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGGTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTT                         <  SRR19973070.392899/150‑1 (MQ=60)
                                                                                                                           CGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGTAACGCGTTTCGAAAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATAACTTCCCAGAACCTGAGTGTAATGTTC                        >  SRR19973070.591372/1‑150 (MQ=60)
                                                                                                                               CACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGAT                    >  SRR19973070.1216443/1‑150 (MQ=60)
                                                                                                                                ACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATG                   <  SRR19973070.2778066/150‑1 (MQ=60)
                                                                                                                                  AAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                 >  SRR19973070.480606/1‑150 (MQ=60)
                                                                                                                                   AAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAT                >  SRR19973070.266966/1‑150 (MQ=60)
                                                                                                                                   AAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATTAAT                <  SRR19973070.1213347/150‑1 (MQ=60)
                                                                                                                                      AGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGA             >  SRR19973070.880150/1‑150 (MQ=60)
                                                                                                                                      AGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGA             >  SRR19973070.419478/1‑150 (MQ=60)
                                                                                                                                        AGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAT           >  SRR19973070.2117113/1‑150 (MQ=60)
                                                                                                                                        AGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTAATGTTATGATGAATCGAAT           >  SRR19973070.1733130/1‑150 (MQ=60)
                                                                                                                                          CGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTT         <  SRR19973070.921868/150‑1 (MQ=60)
                                                                                                                                          CGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTT         <  SRR19973070.599174/150‑1 (MQ=60)
                                                                                                                                           GCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCAGATGAATCGAATTTA        >  SRR19973070.1099359/1‑150 (MQ=60)
                                                                                                                                              ATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCTTCTGGAAGATGACAACTGGCGCTTGCTGGAAAGGGATATCGTCATGGATAAATCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGACTCGAATTTATGC     >  SRR19973070.604466/1‑150 (MQ=60)
                                                                                                                                              ATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC     >  SRR19973070.260913/1‑150 (MQ=60)
                                                                                                                                              ATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC     >  SRR19973070.260667/1‑150 (MQ=60)
                                                                                                                                              ATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCGTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC     <  SRR19973070.660158/150‑1 (MQ=60)
                                                                                                                                                GACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT   >  SRR19973070.2552698/1‑150 (MQ=60)
                                                                                                                                                 AAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTG  >  SRR19973070.2346909/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                    
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTG  >  CP009273/2663992‑2664286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑