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 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 661,892 A→C G594G (GGT→GGG mrdA ← peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273661,8920AC98.2% 182.0 / ‑3.8 56G594G (GGT→GGGmrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (32/23);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

TTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTC  >  NZ_CP009273/661804‑661978
                                                                                        |                                                                                      
tttCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGcaccac                                                                                       <  1:267994/90‑1 (MQ=255)
   ccGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                    <  2:356943/90‑1 (MQ=255)
   ccGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                    <  2:384354/90‑1 (MQ=255)
                  tgttgtTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATgg                                                                     >  2:1386487/1‑90 (MQ=255)
                   gttgttATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGc                                                                    >  2:412543/1‑90 (MQ=255)
                      gttATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGAc                                                                 >  1:214957/1‑90 (MQ=255)
                        tATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAg                                                               >  1:1191691/1‑90 (MQ=255)
                          tCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGcc                                                             <  1:182989/90‑1 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCg                                                                   <  2:1471451/83‑1 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCg                                                                   >  1:1471451/1‑83 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGcc                                                             >  2:1323090/1‑89 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGcc                                                             <  1:1323090/89‑1 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCa                                                            >  1:1537791/1‑90 (MQ=255)
                           cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCa                                                            >  2:1225236/1‑90 (MQ=255)
                                aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  1:1554770/1‑90 (MQ=255)
                                aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  2:253919/1‑90 (MQ=255)
                                  cATATTGTGGTCGAGAATCTGGCGCACCAGTGTCCCAACCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  1:993591/90‑1 (MQ=255)
                                  cATAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  1:921984/90‑1 (MQ=255)
                                  cATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  2:656765/90‑1 (MQ=255)
                                  cATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  2:1406628/90‑1 (MQ=255)
                                    tAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  1:887133/90‑1 (MQ=255)
                                    tAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  2:688772/90‑1 (MQ=255)
                                    tAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  2:1108810/90‑1 (MQ=255)
                                    tAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   >  2:1389007/1‑90 (MQ=255)
                                      aTGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                    >  2:1227619/1‑55 (MQ=255)
                                      aTGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                    <  1:1227619/55‑1 (MQ=255)
                                       tgtgGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAttgt                                                <  1:458079/90‑1 (MQ=255)
                                           gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  2:323892/1‑90 (MQ=255)
                                           gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  1:796173/1‑90 (MQ=255)
                                           gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  1:276483/1‑90 (MQ=255)
                                           gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  1:1376879/1‑90 (MQ=255)
                                               aGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGgc                                        >  1:1486172/1‑90 (MQ=255)
                                                ggATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGgcg                                       >  2:1235937/1‑90 (MQ=255)
                                                    cTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCaaa                                   >  1:521571/1‑90 (MQ=255)
                                                      ggCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAgg                                 <  2:1537791/90‑1 (MQ=255)
                                                         gcATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCgg                              >  1:1393471/1‑90 (MQ=255)
                                                          cATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGt                             >  1:657374/1‑90 (MQ=255)
                                                              aGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                         <  1:1356237/90‑1 (MQ=255)
                                                              aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                         <  2:611802/90‑1 (MQ=255)
                                                              aGTATACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                         <  1:785191/90‑1 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTtgtg                 >  1:892324/1‑90 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTtgtg                 >  2:212594/1‑90 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTtgtg                 >  1:289539/1‑90 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAAGCGGTCATCAGTTtgtg                 >  1:508945/1‑90 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGAAAAGGCg                               <  1:1146850/76‑1 (MQ=255)
                                                                      aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGAAAAGGCg                               >  2:1146850/1‑76 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTg                                                        <  1:761197/48‑1 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTg                                                        >  2:761197/1‑48 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatca                        <  2:304697/80‑1 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatca                        >  1:304697/1‑80 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              <  2:229354/90‑1 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              >  1:229354/1‑90 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              >  1:1009944/1‑90 (MQ=255)
                                                                         cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTc              >  1:562945/1‑90 (MQ=255)
                                                                         cgccgGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              >  2:1333683/1‑90 (MQ=255)
                                                                           ccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCAc            >  2:501596/1‑90 (MQ=255)
                                                                              gACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTa         <  2:657374/90‑1 (MQ=255)
                                                                                     cccaccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTc  <  1:1524693/89‑1 (MQ=255)
                                                                                     acccccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTc  <  1:428476/90‑1 (MQ=255)
                                                                                        |                                                                                      
TTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTC  >  NZ_CP009273/661804‑661978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ 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
TTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  >  NZ_CP009273/661804‑661988
                                                                                        |                                                                                                
TTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAG                                                                                       <  SRR3722078.271266/100‑1 (MQ=60)
   agcgtcagatgtgtataagagacagACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                                    <  SRR3722078.1246959/75‑1 (MQ=60)
            GATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGAC                                                                           >  SRR3722078.217558/1‑100 (MQ=60)
              TCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAG                                                                         >  SRR3722078.1210337/1‑100 (MQ=60)
                 GTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCA                                                                      >  SRR3722078.1495458/1‑100 (MQ=60)
                 GTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCA                                                                      >  SRR3722078.1562856/1‑100 (MQ=60)
                      GTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGC                                                                 >  SRR3722078.1580084/1‑100 (MQ=60)
                          TCCCCCAGAATAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAT                                                             <  SRR3722078.1344241/100‑1 (MQ=60)
                          TCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAT                                                             <  SRR3722078.185161/100‑1 (MQ=60)
                               cgtcggcagcgtcagatgtgtataagagacagGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTG                                                        <  SRR3722078.770979/68‑1 (MQ=60)
                                 GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                      >  SRR3722078.806642/1‑100 (MQ=60)
                                 GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                      >  SRR3722078.279852/1‑100 (MQ=60)
                                 GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                      >  SRR3722078.1399082/1‑100 (MQ=60)
                                  CATATTGTGGTCGAGAATCTGGCGCACCAGTGTCCCAACCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                     <  SRR3722078.1008040/100‑1 (MQ=60)
                                  CATAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCACCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                     <  SRR3722078.934598/100‑1 (MQ=60)
                                    TAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGG                                                   <  SRR3722078.899052/100‑1 (MQ=60)
                                     AATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGC                                                  >  SRR3722078.1510405/1‑100 (MQ=60)
                                       TGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGC                                                <  SRR3722078.463550/100‑1 (MQ=60)
                                          GGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAA                                             >  SRR3722078.527737/1‑100 (MQ=60)
                                               AGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGG                                        >  SRR3722078.1415990/1‑100 (MQ=60)
                                                GGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGT                                       >  SRR3722078.664904/1‑100 (MQ=60)
                                                        ccagTCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGAAAAGGCGGTCATCAGTT                               <  SRR3722078.1164507/96‑1 (MQ=60)
                                                            TCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTG                           >  SRR3722078.904319/1‑100 (MQ=60)
                                                            TCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTG                           >  SRR3722078.293069/1‑100 (MQ=60)
                                                            TCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAAGCGGTCATCAGTTTGTG                           >  SRR3722078.514970/1‑100 (MQ=60)
                                                              AGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                         <  SRR3722078.1378011/100‑1 (MQ=60)
                                                              AGTATACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                         <  SRR3722078.795425/100‑1 (MQ=60)
                                                               GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                        >  SRR3722078.232134/1‑100 (MQ=60)
                                                               GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                        >  SRR3722078.308412/1‑100 (MQ=60)
                                                               GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                        >  SRR3722078.1024790/1‑100 (MQ=60)
                                                               GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTC                        >  SRR3722078.569552/1‑100 (MQ=60)
                                                                               ACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCA        >  SRR3722078.1369668/1‑100 (MQ=60)
                                                                                     CCCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  <  SRR3722078.1549542/100‑1 (MQ=60)
                                                                                     ACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  <  SRR3722078.433640/100‑1 (MQ=60)
                                                                                        |                                                                                                
TTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  >  NZ_CP009273/661804‑661988

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