Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,756,415 Δ57,799 bp [intA]emrB 53 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2756415 2814213 57799 81 [0] [1] 82 [intA]–emrB [intA],yfjH,alpA,yfjI,yfjJ,abpB,abpA,yfjM,rnlA,rnlB,yfjP,yfjQ,yfjR,ypjK,yfjS,yfjT,yfjU,yfjW,yfjX,yfjY,ypjJ,yfjZ,ypjF,ypjA,ypjB,ypjC,ileY,glaH,lhgD,gabD,gabT,gabP,glaR,kbp,yqaE,ygaV,ygaP,stpA,alaE,ygaC,ygaM,nrdH,nrdI,nrdE,nrdF,proV,proW,proX,ygaZ,ygaH,mprA,emrA,emrB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 27564140 (0.000)81 (1.460) 77/488 0.0 100% coding (256/1242 nt) intA integrase IntA
?NC_000913 2814214 = 0 (0.000)intergenic (+60/+4) emrB/luxS multidrug efflux pump membrane subunit EmrB/S‑ribosylhomocysteine lyase

CCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACATTGATCCTCAGGAACATCAGAAAGAACAGGTGAGGAATTCTCAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTA  >  NC_000913/2756172‑2756662
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
ccAAGCCGTTAACTGATACGGCAATCCAAGCCGCCAAACCTCAAGATGCCGATTACCAGCTTTATGACGGTGAAGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGAAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGCACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgact                                                                                                                                                                                                                                                    >  2:289605‑M1/1‑243 (MQ=255)
 cAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATGGGCGAAAGACActgactaga                                                                                                                                                                                                                                                 <  2:365257‑M1/251‑10 (MQ=255)
  aaGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgact                                                                                                                                                                                                                                                    >  1:93027‑M1/1‑241 (MQ=255)
         tAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAAGCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcag                                                                                                                                                                                                                                          <  1:197137‑M1/250‑17 (MQ=255)
           aCTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTACGAGCCGAATCTAAAGTTTTATTGGCGAACGACActgactagctgtgcagttc                                                                                                                                                                                                                                       >  2:416859‑M1/1‑232 (MQ=255)
              gATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctg                                                                                                                                                                                                                                    <  1:373821‑M1/251‑23 (MQ=255)
                     aaaTCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActtactagatgtgcagttccttcaacttt                                                                                                                                                                                                                             >  2:102593‑M1/1‑222 (MQ=255)
                         cAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGGGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAGTTCCGTTACTATCGGCCTTTGACCAAGCAGCAAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTGTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctct                                                                                                                                                                                                                          <  1:225765‑M1/250‑33 (MQ=255)
                           aaGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttc                                                                                                                                                                                                                       >  2:339261‑M1/1‑216 (MQ=255)
                            aGCCGCCAAACCTACAGACGCCGATTACCAGCTTTATGACGGTGACCGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATTCACGTAAACTCAGAGCCGGATCTAAAGTTTTATTGGCTAAAGACActgactagatgtgcagttcctgcaacttctctttcg                                                                                                                                                                                                                      >  2:109139‑M1/1‑215 (MQ=255)
                                  cAAACCTAAAGCTGCCGATTACCAGCTTTATGACGGTGACGGGCTTAATCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCCTTGACCCAGCAGCGAACCCAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAAATCCGAGCCGACTCTAAAGTTTTATTGGCGAAAGCCActgactagatgtgccgtttctgcaacttctctttcggccgtg                                                                                                                                                                                                                >  2:78677‑M1/1‑209 (MQ=255)
                                  cAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtg                                                                                                                                                                                                                <  1:92000‑M1/251‑43 (MQ=255)
                                   aaaCCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTCTCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGCCActgactagctgtgcagttcctgcacctttgctttaggcagcgc                                                                                                                                                                                                               >  2:18936‑M1/1‑208 (MQ=255)
                                   aaaCCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAACGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcattgc                                                                                                                                                                                                               >  1:263370‑M1/1‑208 (MQ=255)
                                   aaaCCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctatttcggcagtgg                                                                                                                                                                                                               >  2:334777‑M1/1‑208 (MQ=255)
                                     aCCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgcca                                                                                                                                                                                                             <  2:21460‑M1/251‑46 (MQ=255)
                                           aGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgggcagttcctgcaccttctctttcggcagtgccagttctt                                                                                                                                                                                                       >  2:359342‑M1/1‑200 (MQ=255)
                                            gATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttattt                                                                                                                                                                                                      >  2:263974‑M1/1‑199 (MQ=255)
                                                  gATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgct                                                                                                                                                                                                <  1:297103‑M1/251‑59 (MQ=255)
                                                  gATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgct                                                                                                                                                                                                <  1:263974‑M1/251‑59 (MQ=255)
                                                  gATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAAACAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgacagttcttcgttgct                                                                                                                                                                                                >  2:92972‑M1/1‑193 (MQ=255)
                                                     tACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTCAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgtt                                                                                                                                                                                             >  1:52234‑M1/1‑190 (MQ=255)
                                                        cAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgat                                                                                                                                                                                          <  2:140456‑M1/251‑65 (MQ=255)
                                                        cAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaac                                                                                                                                                                                                                                <  2:145191‑M1/213‑27 (MQ=255)
                                                        cAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaac                                                                                                                                                                                                                                >  1:145191‑M1/1‑187 (MQ=255)
                                                               aTGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacg                                                                                                                                                                                   >  1:321860‑M1/1‑180 (MQ=255)
                                                                     gtgacggGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgt                                                                                                                                                                             <  1:217036‑M1/251‑78 (MQ=255)
                                                                       gacggGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcggtgctgttgatgcgtccgtcacgttc                                                                                                                                                                           >  2:131662‑M1/1‑172 (MQ=255)
                                                                        acggGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttcc                                                                                                                                                                          <  2:161646‑M1/251‑81 (MQ=255)
                                                                        acggGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCCTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttcc                                                                                                                                                                          >  2:338422‑M1/1‑171 (MQ=255)
                                                                           ggCTTACTCTGTTAATCCCGTACAGTGGCAGTAAGCTTTGGCAATTCCGTTAATCTCGGCCTTTGACCCAGCAGCGAAACAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgatgttgatgcgtacgtcccgttcc                                                                                                                                                                          >  2:136030‑M1/1‑168 (MQ=255)
                                                                               tACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaat                                                                                                                                                                     <  1:159225‑M1/249‑86 (MQ=255)
                                                                                  tctGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGAAGCGAACCAAACAGAGCTTCGGTGCCTATCATGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagctgtgcagttcctgccacttctctttcggcagtgccagttcttcgttgctgttgatgcgtccgtcacgttccagcatgctac                                                                                                                                                                >  1:416232‑M1/1‑161 (MQ=255)
                                                                                       tAATCAAGTACAGTGGCAATACGCTTTGGCAATTCCGTTACTATCGGCCCTTGACCAAGCACCGAACCCAACAGCGCTTCGGTGCCTAACCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGCATCTAACGTTTTATTGGCGCAAGACActgactaaatgtgcagttcctgcaccttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcaccttccagaatgctacgcgcc                                                                                                                                                           >  2:192418‑M1/1‑156 (MQ=255)
                                                                                             aGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcc                                                                                                                                                     >  1:130389‑M1/1‑150 (MQ=255)
                                                                                                ccAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcccgttccagaatgctacgcgcaatatcctgc                                                                                                                                                  >  2:200546‑M1/1‑147 (MQ=255)
                                                                                                     ggCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttc                                                                                                                                             >  1:314370‑M1/1‑142 (MQ=255)
                                                                                                     ggCAGTAAGCTTTGGCAACTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTTTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactatatgtgcagttcctgcaacttctctttcggcagtgccagttcgtcgttgctggagatgcgtacgtcacgttccagaagtctacgcgcaagatcctgcgcttc                                                                                                                                             <  2:26736‑M1/251‑110 (MQ=255)
                                                                                                       cAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcct                                                                                                                                           >  2:100710‑M1/1‑140 (MQ=255)
                                                                                                             gCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacg                                                                                                                                     <  1:166237‑M1/251‑118 (MQ=255)
                                                                                                                ttGGCAATTCCGTTACTATCGGCCTTTGGCCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCGGCCGTTTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttcgctttcggcagtgccatttctgcgttgctgttgatgcggacgtcacgttccagaatgctacgcgcaatatcctgcgctgccggcaacgag                                                                                                                                   <  2:418452‑M1/250‑120 (MQ=255)
                                                                                                                    cAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttccggcaacgagtgcat                                                                                                                              <  1:123477‑M1/251‑125 (MQ=255)
                                                                                                                           gTTACTATCGGCCGTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTAGCCTGCCGTCTCGCTTTCTGACGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACAgtgactagatgtgcagttcctgcaacttctctttcgggagttccagtgcttcgttgctgtggatgcggacgtcacgttccagaatgcgacgggcaatatcctgcgcttcctgcaacgagtgcatcgggtaa                                                                                                                       <  1:289605‑M1/251‑132 (MQ=255)
                                                                                                                                   cGGCCTGTGACCACGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACGCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgaccagatgtgcagtgcctgcaacttcactctcggcagggccagttcttcgttgctggtcatgcgtatgtcacgttccagaatgctacgcgcaatatccggcgcttcccgcaaccagtgcatctggtaagtgcca                                                                                                                 <  2:272995‑M1/249‑138 (MQ=255)
                                                                                                                                    ggCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacac                                                                                                              <  1:53290‑M1/251‑141 (MQ=255)
                                                                                                                                         ttGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggt                                                                                                          <  1:102593‑M1/250‑145 (MQ=255)
                                                                                                                                           gcccaaGCAGCGAACCAATCAGAGCGTCGGTGGCTAGCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActtactagatgtgcagttcctgcaacttcgctgtcggcggtgccagttcgtcgttgctgttgatgcgtacgtcacgttccagaatgcttcgcgcaatatcctgcgcggccggcaacgggtgcatcttggaagtgccacactggtag                                                                                                        <  2:219689‑M1/248‑147 (MQ=255)
                                                                                                                                            aCCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCGGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAAGCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagggccagttcttcgttgctttggatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcgtcctgcaacgagtgcatctggtaagtgccacactggtagac                                                                                                      <  1:221623‑M1/251‑149 (MQ=255)
                                                                                                                                               aaGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgtt                                                                                                   >  1:323590‑M1/1‑100 (MQ=255)
                                                                                                                                                 gcagcGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTATCGATTTCTGATGAACGTAAACTCAGAGCCGAATCTAACGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtcgacgttc                                                                                                  >  1:75132‑M1/1‑98 (MQ=255)
                                                                                                                                                   agcGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggcaagtgccacactggtagacgttcagt                                                                                               >  1:312277‑M1/1‑96 (MQ=255)
                                                                                                                                                            aaCAGAGATTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGCATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaccttctctttcggcagagccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtcagtgccacactggtagacgttcagttccgggatc                                                                                      >  2:374576‑M1/1‑87 (MQ=255)
                                                                                                                                                              cAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctg                                                                                    >  1:15530‑M1/1‑85 (MQ=255)
                                                                                                                                                               agagATTCGGTGCCCTTCCTGCCGTCTCGCATTCTGATGCACGTAAACTCAGAGCCGTATCTAAAGTTTTATTGGCGAAAGACActgtctagatgtgcagttcctgcaacttctctttcggcagtgccagtgcttcgttgctgttgatgctgacgtcacgttccagaatgcgacgcgcaatattctgcgctgcctgcaacgagtgcatcttgtaagtgccacacttgtaagcgttcaggtccgtgatctga                                                                                   <  2:38989‑M1/251‑168 (MQ=255)
                                                                                                                                                                gagCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgat                                                                                  <  1:310022‑M1/251‑169 (MQ=255)
                                                                                                                                                                 agCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgatt                                                                                 >  2:369081‑M1/1‑82 (MQ=255)
                                                                                                                                                                     tCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattct                                                                               <  1:359342‑M1/249‑172 (MQ=255)
                                                                                                                                                                              aTCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcact                                                                    <  1:18936‑M1/251‑183 (MQ=255)
                                                                                                                                                                                ccTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagttccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtggatgttcagtgccgggatctgattctgatcctgcacttt                                                                  <  1:164745‑M1/251‑185 (MQ=255)
                                                                                                                                                                                ccTGCCGTCGCGCGGTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggtagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggagctgattctgatcctgcac                                                                     <  2:93027‑M1/248‑182 (MQ=255)
                                                                                                                                                                                   gCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttc                                                                 <  1:83004‑M1/249‑186 (MQ=255)
                                                                                                                                                                                   gCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacg                                                                                                     >  2:285162‑M1/1‑64 (MQ=255)
                                                                                                                                                                                   gCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacg                                                                                                     <  1:285162‑M1/213‑150 (MQ=255)
                                                                                                                                                                                          cGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctgggaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtc                                                         <  1:49123‑M1/250‑194 (MQ=255)
                                                                                                                                                                                           gCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtctt                                                       <  2:321860‑M1/251‑196 (MQ=255)
                                                                                                                                                                                                 tGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGAAActgactagatgtgcagttcctgccacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattg                                                 >  1:26423‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                  gATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgc                                                <  1:205697‑M1/251‑203 (MQ=255)
                                                                                                                                                                                                   aTGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgcc                                               <  1:339261‑M1/251‑204 (MQ=255)
                                                                                                                                                                                                    ggcACGTAAACTCAGAGCCGAATCGTAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccaggtcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgctgcctgcaacgagtgcatctgttaagtgccagacgggtagacgttcatgtccttgatctgagtcgaatcctgcactttcagcatgtcttccattgccg                                              <  2:189041‑M1/250‑205 (MQ=255)
                                                                                                                                                                                                     gCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgc                                             >  1:100463‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                      cACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGACAGACActaactagatgtgcagttcctgcaacttgtctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgct                                            >  1:431075‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                                        aaTCTAAAGTCTTATTGGCGAAAGACActgacttgatgtgcagttcctgaaacttctctttcggaagtgccagttcttcgttgctggtgatgcgtacgtcacgttccagaatgctacgcggaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtcgacgttcagttccgggatctgattctgatccttcactttcagcacgtcttccattgccgcgttccaggcatcagcaaca                          <  2:52234‑M1/251‑225 (MQ=255)
                                                                                                                                                                                                                              cagTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagaagttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttacaggcatcagcaacacgctgc                    >  2:374535‑M1/2‑21 (MQ=255)
                                                                                                                                                                                                                                gTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttccaggcatcagca                             <  2:301312‑M1/240‑222 (MQ=255)
                                                                                                                                                                                                                                gTTTTATTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttccaggcatcagca                             >  1:301312‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                                     aTTGGCGAAAGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttccaggcatcagcaacacgctgctcatctg             >  2:112248‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                            aaaGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcattatcctgcgcttcatgcaacgagtgcatctggtcagtgcaacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcccgtcttccattgccgctttcaaggcatcagcaacacgctgctcatctggcgtccc      >  2:137923‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                            aaaGACActgactagatgtgcagttcctgcaacttctcttacggcagtgccagttcttcgttgcggttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccatggccgctttccaggcatcagcaacacgctggtcatctggcgtacc      <  2:191389‑M1/251‑245 (MQ=255)
                                                                                                                                                                                                                                             aaGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttccaggcatcagcaacacgctgctcatctggcgtacc      >  1:368955‑M1/1‑6 (MQ=255)
                                                                                                                                                                                                                                              aGACActgactagatgtgcagttcctgcaacttctctttcggcagtgccagttcttcgttgatgttgatgagtacgtaacgttacagaatgctacgcgcaatatcctgcgattcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgatcctgcactttcagcacgtcttccattgccgctttccaggcatcagcaacacgctgctcatctggcgtacca     >  1:383901‑M1/1‑5 (MQ=255)
                                                                                                                                                                                                                                                 cactgactagatgtgcagttcctgcaacttctctttaggcagtgccagttcttcgttgctgttgatgcgtacgtcacgttccagaatgctacgcgcaatatcctgcgcttcctgcaacgagtgcatctggtaagtgccacactggtagacgttcagttccgggatctgattctgttcctgcactttcagcacgtcttccattgccgctttccaggcatcagcaccacgctgctcctctggcgtaccaatc  >  2:321547‑M1/1‑2 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
CCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTTGGCAATTCCGTTACTATCGGCCTTTGACCAAGCAGCGAACCAAACAGAGCTTCGGTGCCTATCCTGCCGTCTCGCTTTCTGATGCACGTAAACTCAGAGCCGAATCTAAAGTTTTATTGGCGAAAGACATTGATCCTCAGGAACATCAGAAAGAACAGGTGAGGAATTCTCAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTA  >  NC_000913/2756172‑2756662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: