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 1,812,807 (C)7→8 intergenic (‑50/+55) celF ← / ← chbR 6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,812,8001.C94.5% 193.0 / 4.5 55intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (28/24);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATT  >  NZ_CP009273/1812713‑1812884
                                                                                        |                                                                                    
gctgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcc                                                                                    <  1:390917/90‑1 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:2922934/1‑90 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:2923564/1‑90 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:487187/1‑90 (MQ=255)
    ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccc                                                                                >  2:1010417/1‑90 (MQ=255)
    ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccc                                                                                >  1:161619/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  2:2160544/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  2:1071581/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  1:1177082/1‑90 (MQ=255)
      cccGCCACCAATAGTGACGATTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:2456164/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:924116/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  2:339157/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:1494380/90‑2 (MQ=255)
                 tAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaata                                                                    <  1:1430596/90‑1 (MQ=255)
                    tGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATc                                                                 >  2:1358185/1‑90 (MQ=255)
                     gacgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                >  1:2478365/1‑90 (MQ=255)
                      acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                               <  1:1358185/90‑1 (MQ=255)
                       cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAg                                                              >  2:2667346/1‑90 (MQ=255)
                       cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAg                                                              >  2:2693513/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2950454/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2510078/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2780914/1‑90 (MQ=255)
                         acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                            <  1:100854/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  1:1010417/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  2:2182773/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  1:498124/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  2:2255958/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:2922934/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:1071581/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  2:887396/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  2:1352458/90‑1 (MQ=255)
                                  ttCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                   <  2:513433/90‑1 (MQ=255)
                                  ttCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                   <  2:324875/90‑1 (MQ=255)
                                            tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                         <  1:2160544/90‑1 (MQ=255)
                                                ttCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                     <  2:1133889/90‑1 (MQ=255)
                                                 tctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                    <  1:1103060/90‑1 (MQ=255)
                                                  ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  2:1529093/90‑1 (MQ=255)
                                                  ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  1:1578484/90‑1 (MQ=255)
                                                    cccTTCAGTTCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     >  2:1472424/1‑70 (MQ=255)
                                                    cccTTCAGTTCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:1472424/70‑1 (MQ=255)
                                                      cTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                               >  1:103686/1‑90 (MQ=255)
                                                       ttCAGTACAGAATACTGATATCTGGCATATCTCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAAc                            >  2:506462/1‑90 (MQ=255)
                                                               aGAATACTGATATCTGGCATATCTCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                    <  1:506462/90‑1 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  2:94533/1‑90 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  1:625974/1‑90 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     <  1:487187/90‑1 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  1:1010828/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                      <  2:1697860/52‑1 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                      >  1:1697860/1‑52 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                >  1:1048844/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                >  1:2527543/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa               >  1:1363958/1‑90 (MQ=255)
                                                                       gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa              >  2:705738/1‑90 (MQ=255)
                                                                       gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa              >  1:2374362/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                 <  2:945454/52‑1 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                 >  1:945454/1‑52 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:2428644/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:2898748/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:168942/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  2:971676/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  2:1316956/1‑90 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   >  2:1227639/1‑64 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  1:1227639/64‑1 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc          >  1:1104816/1‑89 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc          <  2:1104816/89‑1 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa         >  1:253439/1‑90 (MQ=255)
                                                                                   atCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGCCAGATATTTGAATTGTCAGGGATAACGACTTACTGATTAAATTCAGTCAAtt  >  1:1736406/1‑90 (MQ=255)
                                                                                   atCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAtt  >  1:1323175/1‑90 (MQ=255)
                                                                                        |                                                                                    
GCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATT  >  NZ_CP009273/1812713‑1812884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
GTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  NZ_CP009273/1812708‑1812903
                                                                                                    |                                                                                                
GTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCC                                                                                                  >  SRR3722030.1194526/1‑100 (MQ=60)
     GCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGAC                                                                                              <  SRR3722030.396768/100‑1 (MQ=60)
           CCCGCCACCAATAGTGACGATTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                        <  SRR3722030.2501045/100‑1 (MQ=60)
           CCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                        <  SRR3722030.1518332/100‑1 (MQ=60)
           CCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                        <  SRR3722030.938267/100‑1 (MQ=60)
                CACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCC                                                                                   >  SRR3722030.2523655/1‑100 (MQ=60)
                      TAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAC                                                                             <  SRR3722030.1453262/100‑1 (MQ=60)
                           ACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAC                                                                        <  SRR3722030.1379327/100‑1 (MQ=60)
                              ACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGA                                                                     <  SRR3722030.102343/100‑1 (MQ=60)
                                TTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATA                                                                   <  SRR3722030.1025698/100‑1 (MQ=60)
                                TTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATA                                                                   <  SRR3722030.505751/100‑1 (MQ=60)
                                     taagagacagCATAATTTCTCCCTTCAGTTCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGA                                                              <  SRR3722030.1495916/90‑1 (MQ=60)
                                     ATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGA                                                              <  SRR3722030.1087632/100‑1 (MQ=60)
                                     ATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGA                                                              <  SRR3722030.2976816/100‑1 (MQ=60)
                                                 TAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTAT                                                  <  SRR3722030.2199267/100‑1 (MQ=60)
                                                 TAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTAT                                                  >  SRR3722030.105205/1‑100 (MQ=60)
                                                      TCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCtgtctcttatacacatctgacgctgccgacgacatggcc                                             >  SRR3722030.1057861/1‑61 (MQ=60)
                                                      TCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGA                                             <  SRR3722030.1119525/100‑1 (MQ=60)
                                                       tgtgtataagagacagATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGAC                                            <  SRR3722030.1245730/84‑1 (MQ=60)
                                                       CTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGAC                                            <  SRR3722030.1604058/100‑1 (MQ=60)
                                                           CTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAC                                        >  SRR3722030.635682/1‑100 (MQ=60)
                                                           CTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAC                                        >  SRR3722030.1026114/1‑100 (MQ=60)
                                                                GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGctgtctcttatacacatctgacgctgcc                                   >  SRR3722030.1725467/1‑72 (MQ=60)
                                                                GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                   >  SRR3722030.1064611/1‑100 (MQ=60)
                                                                GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                   >  SRR3722030.2573810/1‑100 (MQ=60)
                                                                GTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTA                                  >  SRR3722030.1385239/1‑100 (MQ=60)
                                                                  ACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAA                                 >  SRR3722030.2417537/1‑100 (MQ=60)
                                                                    AGAATACTGATATCTGGCATATCTCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTATATTC                             <  SRR3722030.514215/100‑1 (MQ=60)
                                                                     GAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGctgtctcttatacacatctgacgctgcc                              >  SRR3722030.959896/1‑72 (MQ=60)
                                                                     GAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATT                              >  SRR3722030.2472929/1‑100 (MQ=60)
                                                                     GAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATT                              >  SRR3722030.2952208/1‑100 (MQ=60)
                                                                     GAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATT                              <  SRR3722030.494642/100‑1 (MQ=60)
                                                                     GAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATT                              >  SRR3722030.171373/1‑100 (MQ=60)
                                                                       ATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA                            >  SRR3722030.257095/1‑100 (MQ=60)
                                                                       ATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA                            >  SRR3722030.1121306/1‑100 (MQ=60)
                                                                              TATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGCCAGATATTTGAATTGTCAGGGATAACGACTTACTGATTAAATTCAGTCAATT                     >  SRR3722030.1765054/1‑100 (MQ=60)
                                                                              TATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATT                     >  SRR3722030.1343483/1‑100 (MQ=60)
                                                                                       TATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGAT            >  SRR3722030.2225447/1‑100 (MQ=60)
                                                                                                 CCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  <  SRR3722030.1600815/100‑1 (MQ=60)
                                                                                                 CCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  SRR3722030.1079837/1‑100 (MQ=60)
                                                                                                    |                                                                                                
GTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  NZ_CP009273/1812708‑1812903

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: