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 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,848,051 T→G F257C (TTT→TGT)  emrD → multidrug efflux MFS transporter EmrD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,848,0510TG96.8% 577.9 / 3.6 158F257C (TTT→TGT) emrDmultidrug efflux MFS transporter EmrD
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (102/51);  total (104/54)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 158 total aligned reads displayed.

GCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCT  >  NZ_CP009273/3847967‑3848137
                                                                                    |                                                                                      
gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                                                   >  1:1670594/1‑90 (MQ=255)
        gTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATgg                                                                           >  2:1026202/1‑90 (MQ=255)
               tGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCg                                                                    >  1:741409/1‑90 (MQ=255)
                  gcgcGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGAc                                                                 >  1:1545404/1‑90 (MQ=255)
                    gcgGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGt                                                               >  1:2453824/1‑90 (MQ=255)
                      ggTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTcc                                                             <  1:1129586/90‑1 (MQ=255)
                            aGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATaa                                                       >  2:263377/1‑90 (MQ=255)
                               gCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACg                                                    <  1:697621/90‑1 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACg                                                                >  2:2191433/1‑76 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTcc                                                             >  2:1586482/1‑79 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTcc                                                             <  2:238955/79‑1 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTcc                                                             >  1:238955/1‑79 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTcc                                                             <  1:1586482/79‑1 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCTAATAAACGCt                                                  >  1:2944068/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAAc                                                     <  2:735469/87‑1 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:859126/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:368734/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:3022850/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:525934/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:2671810/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:2471527/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:605305/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:61275/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1027220/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:2586479/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  <  2:2453824/90‑1 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1350311/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1437179/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1492353/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1651300/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1837201/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1339917/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1432963/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1437001/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:789871/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1603667/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1650477/1‑90 (MQ=255)
                                 tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1789339/1‑90 (MQ=255)
                                   aGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtc                                                <  1:2735843/90‑1 (MQ=255)
                                     cAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                              <  1:2117871/90‑1 (MQ=255)
                                       gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCa                                            >  2:905137/1‑90 (MQ=255)
                                       gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCa                                            >  2:1106084/1‑90 (MQ=255)
                                       gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCa                                            >  1:1706834/1‑90 (MQ=255)
                                         aTGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACg                                          >  2:19976/1‑90 (MQ=255)
                                         aTGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACg                                          >  2:16744/1‑90 (MQ=255)
                                            aCGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACg                                                    <  2:2049345/77‑1 (MQ=255)
                                            aCGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTa                                       >  2:2325934/1‑90 (MQ=255)
                                            aCGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTa                                       >  1:1201459/1‑90 (MQ=255)
                                            aCGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTa                                       >  2:1301416/1‑90 (MQ=255)
                                              ggTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATa                                                        <  1:1260240/71‑1 (MQ=255)
                                              ggTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATa                                                        >  2:1260240/1‑71 (MQ=255)
                                              ggTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAt                                     >  1:76683/1‑90 (MQ=255)
                                              ggTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAt                                     >  1:1961431/1‑90 (MQ=255)
                                              ggTCAGTATTTTGGTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAt                                     >  1:1411433/1‑90 (MQ=255)
                                                tCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgt                                   >  2:657782/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGTAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:28441/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTcc                                                             >  2:1319730/1‑63 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCaa                                                          <  1:929912/66‑1 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCaa                                                          >  2:929912/1‑66 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATaa                                                       <  1:639738/69‑1 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATaa                                                       >  2:639738/1‑69 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGc                                                   <  1:2154859/73‑1 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                              <  1:1936017/78‑1 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:1771506/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:238945/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:905466/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:1472467/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:2330838/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:208190/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:2509085/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:471531/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:781139/1‑90 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCG‑TGTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTcc                                             <  1:1038090/78‑1 (MQ=255)
                                                 cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCG‑TGTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTcc                                             >  2:1038090/1‑78 (MQ=255)
                                                    tATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCa                               >  1:287593/1‑90 (MQ=255)
                                                       tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  1:471531/90‑1 (MQ=255)
                                                       tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  1:885640/90‑1 (MQ=255)
                                                       tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  2:2330838/90‑1 (MQ=255)
                                                       tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  2:1911495/90‑1 (MQ=255)
                                                       tttGTTTATTCTGCAGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  1:1790841/90‑1 (MQ=255)
                                                        ttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCa                               >  2:1386987/1‑86 (MQ=255)
                                                        ttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCa                               <  1:1386987/86‑1 (MQ=255)
                                                        ttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCg                           >  2:24165/1‑90 (MQ=255)
                                                        ttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCg                           >  1:2297433/1‑90 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGt                                                               >  2:2658674/1‑49 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCAc                                           <  2:1561447/69‑1 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  1:308675/84‑1 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            >  2:308675/1‑84 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                      >  1:1225235/1‑90 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                      >  2:2886009/1‑90 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                      >  2:2949967/1‑90 (MQ=255)
                                                             tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                      >  2:1925048/1‑90 (MQ=255)
                                                               ttCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATct                    >  1:53604/1‑90 (MQ=255)
                                                                   gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCa                               >  2:1093183/1‑75 (MQ=255)
                                                                   gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCa                               <  1:1093183/75‑1 (MQ=255)
                                                                   gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATctgctg                <  1:2925734/90‑1 (MQ=255)
                                                                    ccGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                              >  1:700080/1‑59 (MQ=255)
                                                                    ccGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCtgctgc               >  2:1794582/1‑90 (MQ=255)
                                                                            ggCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTgctggc       <  1:2586479/90‑1 (MQ=255)
                                                                                 cGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTggctggct  <  1:1492353/90‑1 (MQ=255)
                                                                                    |                                                                                      
GCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCT  >  NZ_CP009273/3847967‑3848137

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGAT  >  NZ_CP009273/3847957‑3848150
                                                                                              |                                                                                                   
TTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTG                                                                                                >  SRR3722030.1697724/1‑100 (MQ=60)
 TTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGG                                                                                               <  SRR3722030.108062/100‑1 (MQ=60)
 TTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGG                                                                                               <  SRR3722030.3045604/100‑1 (MQ=60)
 TTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGG                                                                                               <  SRR3722030.714674/100‑1 (MQ=60)
   GAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCG                                                                                             <  SRR3722030.1049938/100‑1 (MQ=60)
               tacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCG                                                                                 <  SRR3722030.3010873/54‑1 (MQ=60)
               GGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCG                                                                                 >  SRR3722030.752833/1‑100 (MQ=60)
                  gagtataagagacaGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGAC                                                                              <  SRR3722030.782808/86‑1 (MQ=60)
                  GTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGAC                                                                              >  SRR3722030.1570400/1‑100 (MQ=60)
                    GCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGT                                                                            >  SRR3722030.2498640/1‑100 (MQ=60)
                            GCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAA                                                                    <  SRR3722030.1458682/100‑1 (MQ=60)
                             acaGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAA                                                                   <  SRR3722030.2230814/97‑1 (MQ=60)
                              gtcggcagcgtcagatgtgtataagagacagGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAAC                                                                  <  SRR3722030.2707644/69‑1 (MQ=60)
                                atgtgtataagagacagGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGC                                                                <  SRR3722030.1339937/83‑1 (MQ=60)
                                GGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGC                                                                <  SRR3722030.1612193/100‑1 (MQ=60)
                                GGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGC                                                                <  SRR3722030.1146434/100‑1 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCTAATAAACGCT                                                               >  SRR3722030.2998357/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1360604/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.2588293/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.2516684/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.2721060/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1042714/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.3034133/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.3078517/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.426670/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.614688/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1819360/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.62224/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.640135/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1677273/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1629666/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.746814/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.872311/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.662259/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1455674/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                               >  SRR3722030.1459787/1‑100 (MQ=60)
                                 GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCC                                                               >  SRR3722030.242395/1‑100 (MQ=60)
                                   tgtataagagacagGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTC                                                             <  SRR3722030.944158/86‑1 (MQ=60)
                                     aagagacagGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTC                                                           <  SRR3722030.1279034/91‑1 (MQ=60)
                                      cgtcagatgtgtataagagacagGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCC                                                          <  SRR3722030.1107987/77‑1 (MQ=60)
                                      ataagagacagGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCC                                                          <  SRR3722030.649667/89‑1 (MQ=60)
                                       GGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCA                                                         >  SRR3722030.1734673/1‑100 (MQ=60)
                                         GCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACG                                                       <  SRR3722030.708413/100‑1 (MQ=60)
                                          gagacagGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGT                                                      <  SRR3722030.2193460/93‑1 (MQ=60)
                                            GAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGctg                                                    >  SRR3722030.2085663/1‑97 (MQ=60)
                                            GAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTA                                                    >  SRR3722030.1219229/1‑100 (MQ=60)
                                            GAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTA                                                    >  SRR3722030.2189156/1‑100 (MQ=60)
                                             AGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAA                                                   <  SRR3722030.2786372/100‑1 (MQ=60)
                                              GCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAT                                                  >  SRR3722030.1995711/1‑100 (MQ=60)
                                              GCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAT                                                  >  SRR3722030.77865/1‑100 (MQ=60)
                                              GCAGTATGACGGTCAGTATTTTGGTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAT                                                  >  SRR3722030.1433706/1‑100 (MQ=60)
                                               agGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATG                                                 <  SRR3722030.1969715/98‑1 (MQ=60)
                                               agGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCG‑TGTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGT                                                <  SRR3722030.1053733/98‑1 (MQ=60)
                                               CAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATG                                                 <  SRR3722030.1640936/100‑1 (MQ=60)
                                               CAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATG                                                 <  SRR3722030.2155658/100‑1 (MQ=60)
                                                 GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.2995130/1‑100 (MQ=60)
                                                 GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.950619/1‑100 (MQ=60)
                                                 GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.242385/1‑100 (MQ=60)
                                                 GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.2384664/1‑100 (MQ=60)
                                                 GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.2373161/1‑100 (MQ=60)
                                                 GTATAACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                               >  SRR3722030.211157/1‑100 (MQ=60)
                                                    TGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCA                                            >  SRR3722030.291756/1‑100 (MQ=60)
                                                        GGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCctgtctcttatacacatc                                        >  SRR3722030.878545/1‑82 (MQ=60)
                                                        GGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCG                                        >  SRR3722030.2339061/1‑100 (MQ=60)
                                                          TCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGT                                      <  SRR3722030.1186663/100‑1 (MQ=60)
                                                             GTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGctgtctcttataca                                   >  SRR3722030.1422435/1‑86 (MQ=60)
                                                             GTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCtgtctcttat                                   >  SRR3722030.1586706/1‑90 (MQ=60)
                                                             GTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAT                                   >  SRR3722030.1243306/1‑100 (MQ=60)
                                                              gacagTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATC                                  <  SRR3722030.1109514/95‑1 (MQ=60)
                                                              TATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATC                                  <  SRR3722030.1408732/100‑1 (MQ=60)
                                                               ATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCT                                 >  SRR3722030.54426/1‑100 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.2846940/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.313170/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.1850177/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.478703/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.585549/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.1267343/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                               <  SRR3722030.899219/100‑1 (MQ=60)
                                                                 TTTGTTTATTCTGCAGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCGGC                               <  SRR3722030.1820905/100‑1 (MQ=60)
                                                                   TGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTACCTGCTG                             >  SRR3722030.2489492/1‑100 (MQ=60)
                                                                    GTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGctgtctcttatacacatctga                            >  SRR3722030.710904/1‑79 (MQ=60)
                                                                     TTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCC                           <  SRR3722030.17020/100‑1 (MQ=60)
                                                                             GCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCT                   <  SRR3722030.1766624/100‑1 (MQ=60)
                                                                             GCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCT                   <  SRR3722030.2979669/100‑1 (MQ=60)
                                                                                  TTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACctgtctcttatacacatctgacgctgccgacgacatggccag              >  SRR3722030.2022442/1‑58 (MQ=60)
                                                                                  TTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTT              >  SRR3722030.1082274/1‑100 (MQ=60)
                                                                                      GGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTG          <  SRR3722030.2633940/100‑1 (MQ=60)
                                                                                         AGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATG       <  SRR3722030.151991/100‑1 (MQ=60)
                                                                                           CGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTG     <  SRR3722030.1516264/100‑1 (MQ=60)
                                                                                            GTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggtggtcgcc    >  SRR3722030.311296/1‑38 (MQ=60)
                                                                                            GTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGG    <  SRR3722030.2127579/100‑1 (MQ=60)
                                                                                            GTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGG    >  SRR3722030.2718340/1‑100 (MQ=60)
                                                                                              GTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGAT  <  SRR3722030.1371257/100‑1 (MQ=60)
                                                                                              GTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGAT  >  SRR3722030.2594676/1‑100 (MQ=60)
                                                                                              |                                                                                                   
TTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGAT  >  NZ_CP009273/3847957‑3848150

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