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 1,877,816 T→C pseudogene (408/1482 nt) yeaV → BCCT family transporter YeaV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,877,8160TC96.8% 192.0 / ‑1.1 63pseudogene (408/1482 nt)yeaVBCCT family transporter YeaV
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (39/22);  total (40/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 = 7.27e-01

GGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTAC  >  NZ_CP009273/1877729‑1877902
                                                                                       |                                                                                      
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  1:1007596/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:223839/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:1434717/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:422686/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:1312493/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:1014166/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  1:421824/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  1:1211225/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                      >  2:817971/1‑90 (MQ=255)
 gCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGc                                                                                     <  2:614999/90‑1 (MQ=255)
      aTTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCa                                                                                <  1:1199026/90‑1 (MQ=255)
       ttAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCaa                                                                               <  2:771095/90‑1 (MQ=255)
             aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCa                                                                         <  1:1491014/90‑1 (MQ=255)
             aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCa                                                                         <  1:920854/90‑1 (MQ=255)
              tGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCAt                                                                        >  1:1363870/1‑90 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                       >  2:1273959/1‑90 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                       >  2:1250490/1‑90 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                       >  1:1085891/1‑90 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                       >  1:560248/1‑90 (MQ=255)
                gCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATgg                                                                      >  2:395467/1‑90 (MQ=255)
                          ttCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTgg                                                            >  2:1245994/1‑90 (MQ=255)
                          ttCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTgg                                                            >  1:396229/1‑90 (MQ=255)
                          ttCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTgg                                                            >  2:160919/1‑90 (MQ=255)
                           tCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGt                                                           <  2:1428191/90‑1 (MQ=255)
                               gtgCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAt                                                       <  1:422686/90‑1 (MQ=255)
                               gtgCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAt                                                       >  2:870612/1‑90 (MQ=255)
                               gtgCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAt                                                       <  1:1434717/90‑1 (MQ=255)
                                  cGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTg                                                    >  1:130704/1‑90 (MQ=255)
                                         aCAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTcc                                             <  2:421824/90‑1 (MQ=255)
                                         aCAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTcc                                             <  1:1527777/90‑1 (MQ=255)
                                            aaGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCTTGa                                          >  2:1413294/1‑90 (MQ=255)
                                              ggTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATc                                        >  1:1083070/1‑90 (MQ=255)
                                                 cTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGcc                                     >  2:219396/1‑90 (MQ=255)
                                                    aGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                                  <  2:714168/90‑1 (MQ=255)
                                                    aGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                                  <  2:1007596/90‑1 (MQ=255)
                                                    aGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                                  <  1:1250490/90‑1 (MQ=255)
                                                        tttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                              >  2:365295/1‑90 (MQ=255)
                                                        tttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                              >  2:1152671/1‑90 (MQ=255)
                                                        tttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                              >  2:463645/1‑90 (MQ=255)
                                                        tttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                              >  2:608717/1‑90 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTg                                           >  1:745108/1‑76 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTg                                           <  2:745108/76‑1 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCgg                             >  2:9206/1‑90 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCgg                             >  1:1166972/1‑90 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCgg                             >  2:82966/1‑90 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                                  <  2:1400755/79‑1 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                                  >  1:1400755/1‑79 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACt                       <  1:1035620/90‑1 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACt                       <  1:219396/90‑1 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACt                       <  1:1312493/90‑1 (MQ=255)
                                                                attattGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTg                      >  1:1131476/1‑90 (MQ=255)
                                                                attattGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAt                                                       >  1:972827/1‑57 (MQ=255)
                                                                attattGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAt                                                       <  2:972827/57‑1 (MQ=255)
                                                                attattGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTg                      <  1:1088161/90‑1 (MQ=255)
                                                                  tattGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGAc                    >  1:612532/1‑90 (MQ=255)
                                                                        tgctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTc              >  1:547562/1‑90 (MQ=255)
                                                                        tgctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTc              <  1:1245994/90‑1 (MQ=255)
                                                                           tATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCt           >  2:552142/1‑90 (MQ=255)
                                                                                   gCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAttt                                                     <  2:1518591/40‑1 (MQ=255)
                                                                                   gCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAttt                                                     >  1:1518591/1‑40 (MQ=39)
                                                                                   gCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTa   >  2:1525970/1‑90 (MQ=255)
                                                                                    cGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGTTGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTAc  <  2:1360314/90‑1 (MQ=255)
                                                                                    cGTCCGCCCGCCAGGCCCCTGGGGAAAACTGGTCGATTTGATGTTTCCGGTCGCCCCCGTCGGCGCCCTGACCCTTTCCCTTgttttttc  >  1:1410042/1‑88 (MQ=255)
                                                                                       |                                                                                      
GGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTAC  >  NZ_CP009273/1877729‑1877902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
CTTATACGCTGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTACCG  >  NZ_CP009273/1877719‑1877904
                                                                                                 |                                                                                        
CTTATACGCTGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCG                                                                                        >  SRR3722078.1022392/1‑100 (MQ=60)
CTTATACGCTGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCG                                                                                        >  SRR3722078.426911/1‑100 (MQ=60)
CTTATACGCTGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCG                                                                                        >  SRR3722078.1230266/1‑100 (MQ=60)
              TCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCAT                                                                          >  SRR3722078.1385805/1‑100 (MQ=60)
               CATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATG                                                                         >  SRR3722078.1102343/1‑100 (MQ=60)
               CATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATG                                                                         >  SRR3722078.566838/1‑100 (MQ=60)
                ATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGG                                                                        <  SRR3722078.1217841/100‑1 (MQ=60)
                       ATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAAC                                                                 <  SRR3722078.1515322/100‑1 (MQ=60)
                       ATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAAC                                                                 <  SRR3722078.933437/100‑1 (MQ=60)
                          GCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGG                                                              >  SRR3722078.401034/1‑100 (MQ=60)
                                  CTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTG                                                      >  SRR3722078.132213/1‑100 (MQ=60)
                                         GTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCC                                               <  SRR3722078.427784/100‑1 (MQ=60)
                                         GTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCC                                               <  SRR3722078.1458062/100‑1 (MQ=60)
                                              GAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATC                                          >  SRR3722078.1099458/1‑100 (MQ=60)
                                                   ACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCAC                                     <  SRR3722078.1552668/100‑1 (MQ=60)
                                                         GTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGctgt                               >  SRR3722078.754556/1‑96 (MQ=60)
                                                         GTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGG                               >  SRR3722078.1185105/1‑100 (MQ=60)
                                                              AGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCAC                          <  SRR3722078.1270289/100‑1 (MQ=60)
                                                               GCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACC                         >  SRR3722078.1423443/1‑100 (MQ=60)
                                                                CCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTG                        >  SRR3722078.1148800/1‑100 (MQ=60)
                                                                CCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCctgtctcttatacacatctgacg                        >  SRR3722078.986761/1‑77 (MQ=60)
                                                                  TTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGAC                      >  SRR3722078.619601/1‑100 (MQ=60)
                                                                   gatgtgtataagagacagTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACC                     <  SRR3722078.915315/82‑1 (MQ=60)
                                                                        GCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTC                >  SRR3722078.554027/1‑100 (MQ=60)
                                                                         CATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCC               <  SRR3722078.1051030/100‑1 (MQ=60)
                                                                         CATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCC               <  SRR3722078.222048/100‑1 (MQ=60)
                                                                         CATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCC               <  SRR3722078.1333455/100‑1 (MQ=60)
                                                                          ATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCC              <  SRR3722078.1104666/100‑1 (MQ=60)
                                                                                  TGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTT      <  SRR3722078.1265697/100‑1 (MQ=60)
                                                                                   GCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGctgtctcttatacacatctgacgctgccgacgatggcggt     >  SRR3722078.1543357/1‑60 (MQ=60)
                                                                                    CTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTAC    >  SRR3722078.478222/1‑100 (MQ=60)
                                                                                    CTATTAACGGCGTCCGCCCGCCAGGCCCCTGGGGAAAACTGGTCGATTTGATGTTTCCGGTCGCCCCCGTCGGCGCCCTGACCCTTTCCCTTGTTttttc    >  SRR3722078.1432921/1‑95 (MQ=34)
                                                                                      ATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTACCG  >  SRR3722078.1350476/1‑100 (MQ=60)
                                                                                                 |                                                                                        
CTTATACGCTGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTACCG  >  NZ_CP009273/1877719‑1877904

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