Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,518,434 Δ81,761 bp [ykiE]yjjA 74 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4518434 4600194 81761 84 [0] [0] 84 [ykiE]–yjjA [ykiE],insA7,yjhU,yjhF,yjhG,yjhH,yjhI,sgcR,sgcE,sgcA,ryjB,sgcQ,sgcC,sgcB,sgcX,yjhP,yjhQ,topAI,ythB,nanS,nanM,nanC,fimB,fimE,fimA,fimI,fimC,fimD,fimF,fimG,fimH,gntP,uxuA,uxuB,uxuR,yjiC,ytiC,ytiD,idlP,iraD,hypT,iadA,yjiG,yjiH,kptA,yjiJ,yjiK,ytiA,yjiL,yjiM,yjiN,mdtM,ytiE,yjiR,yjiS,yjiT,mcrC,mcrB,mcrB,symE,symR,hsdS,hsdM,hsdR,mrr,yjiA,yjiX,btsT,tsr,lgoT,lgoR,lgoD,opgB,yjjA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45184330 (0.000)84 (0.870) 73/488 0.1 100% coding (62/69 nt) ykiE protein YkiE
?NC_000913 4600195 = 0 (0.000)intergenic (‑6/+43) yjjA/dnaC DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC

ATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCG  >  NC_000913/4518185‑4518677
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aTAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAct                                                                                                                                                                                                                                                    <  1:469125‑M1/251‑3 (MQ=255)
                            gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccg                                                                                                                                                                                                                        >  1:471207‑M1/1‑221 (MQ=255)
                            gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatacc                                                                                                                                                                                                                         <  1:438043‑M1/250‑30 (MQ=255)
                                     gcgcGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacg                                                                                                                                                                                                               <  1:328555‑M1/251‑40 (MQ=255)
                                       gcgGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggg                                                                                                                                                                                                             >  1:147061‑M1/1‑210 (MQ=255)
                                             aGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggatta                                                                                                                                                                                                        >  1:535454‑M1/1‑204 (MQ=255)
                                               aCATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatt                                                                                                                                                                                                     >  1:23361‑M1/1‑202 (MQ=255)
                                                     ggAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTTTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactcttt                                                                                                                                                                                               <  2:734818‑M1/251‑56 (MQ=255)
                                                       aGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttac                                                                                                                                                                                             <  1:348656‑M1/251‑58 (MQ=255)
                                                        gTGCATCAAAAGTTAATTATCGCGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCGActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacc                                                                                                                                                                                            <  1:252216‑M1/251‑59 (MQ=255)
                                                               aaaaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacc                                                                                                                                                                                     <  1:105401‑M1/251‑66 (MQ=255)
                                                               aaaaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacc                                                                                                                                                                                     >  1:567374‑M1/1‑186 (MQ=255)
                                                                 aaGTTAATTATCACGTGTTCATATTCATATGAGACTGGTTATCATTACAATTGGAAATAAAATTGTTTCCAAAAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGATAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGAGGTGGTCTGCTAcgcaaggtggtattgtgtcgcagtataccgccagaaacggggattaatactctttcactgttacccg                                                                                                                                                                                   >  2:46052‑M1/1‑184 (MQ=255)
                                                                 aaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccg                                                                                                                                                                                   <  1:414300‑M1/251‑68 (MQ=255)
                                                                   gTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaata                                                                                                                                                                                                     >  2:568472‑M1/1‑182 (MQ=255)
                                                                   gTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaata                                                                                                                                                                                                     <  1:568472‑M1/231‑50 (MQ=255)
                                                                      aaTTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggcta                                                                                                                                                                               <  1:676193‑M1/250‑72 (MQ=255)
                                                                      aaTTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggcta                                                                                                                                                                               <  2:567374‑M1/250‑72 (MQ=255)
                                                                          aTCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggta                                                                                                                                                                          <  1:310104‑M1/251‑77 (MQ=255)
                                                                          aTCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGGTTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatcctctttacctgttacccggctacggta                                                                                                                                                                          <  2:111647‑M1/251‑77 (MQ=255)
                                                                              cGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTATGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagct                                                                                                                                                                       >  1:217740‑M1/1‑171 (MQ=255)
                                                                                tAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGAATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGACGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatc                                                                                                                                                                    >  2:480863‑M1/1‑169 (MQ=255)
                                                                                     atatTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctattccagt                                                                                                                                                               <  1:228759‑M1/251‑88 (MQ=255)
                                                                                       atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttg                                                                                                                                                             >  1:188127‑M1/1‑162 (MQ=255)
                                                                                           atatGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttga                                                                                                                                                            >  1:755912‑M1/1‑158 (MQ=255)
                                                                                           atatGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCAGGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaaga                                                                                                                                                         <  1:577419‑M1/251‑94 (MQ=255)
                                                                                                agaAGGGTTATCATTACAATTGGAAATAAAAGTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCGAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGGGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatgccgccagaaacggggattaatactctttacctgttagccggctacggtagctatcccagttgaagatcacc                                                                                                                                                    <  2:23361‑M1/251‑99 (MQ=255)
                                                                                                 gaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaccc                                                                                                                                                   >  1:21204‑M1/1‑152 (MQ=255)
                                                                                                       ttATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAAAATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGAATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaac                                                                                                                                             >  2:195291‑M1/1‑146 (MQ=255)
                                                                                                            aTTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaatgtggaattgtgtcgcagtataccgccagaaacgtggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgtta                                                                                                                                        <  2:609380‑M1/251‑111 (MQ=255)
                                                                                                            aTTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttc                                                                                                                                        >  1:220577‑M1/1‑141 (MQ=255)
                                                                                                             ttACAATTGGAAATAAAATTGTTTCAAATAGCCATTTTTAACATGTTGATTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGCGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcacggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttac                                                                                                                                       >  1:263551‑M1/1‑140 (MQ=255)
                                                                                                                 aattgtgaataaaattGTTTCCAATAGACATNTTTAACATGTTGTTTTTCTAAGTGTTGTAATGTAGGTATAAAATGGGATGGGGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTTGTCTGCTActcaaggtggaattttgtcgcagtataccgccagaaacggggattaatactctttccctgttacccggctacgggagctatcccaggtgaaggtctcccacaaactgttgccctg                                                                                                                                   <  2:338565‑M1/244‑116 (MQ=255)
                                                                                                                 aaTTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttaccc                                                                                                                                     >  1:642602‑M1/1‑136 (MQ=255)
                                                                                                                         aTAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatac                                                                                                                           <  1:239179‑M1/251‑124 (MQ=255)
                                                                                                                               ttGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacncacaaactgttacccaggcgcatacggtcc                                                                                                                      >  2:338952‑M1/1‑122 (MQ=255)
                                                                                                                                tatttccACTAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGGCAGAATGACTCCTGATGTGGTCTGCTActcaaggtggaattgtgacgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaaccacaaactgttacccaggcgcatactgtccat                                                                                                                    >  2:559213‑M1/3‑121 (MQ=255)
                                                                                                                                tGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATTATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcccccacaaactgttacccaggcgcatacggtcca                                                                                                                     <  2:193563‑M1/250‑130 (MQ=255)
                                                                                                                                     ccAATAGACATTTTTATTATGTTGTTTTTCTAAGTGTTGTAAGGGAGGGATAAAATGGGATGGAGACTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTCctcaaggtggaattgcgtcgcagtataccgccagaatcggggattaatactctatccctcttacccggctacggtagctatcccagttcaagatcacccacaatctgttacccaggctcatacggtccatcactc                                                                                                               <  2:217740‑M1/251‑136 (MQ=255)
                                                                                                                                     ccAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgc                                                                                                               >  1:668268‑M1/1‑116 (MQ=255)
                                                                                                                                          aGACATTTTTATCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttc                                                                                                           <  1:559213‑M1/250‑140 (MQ=255)
                                                                                                                                          aGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATAGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACCCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatctcagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcg                                                                                                          >  1:261646‑M1/1‑111 (MQ=255)
                                                                                                                                                      cATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttg                                                                                              >  1:546516‑M1/1‑99 (MQ=255)
                                                                                                                                                       aTGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGAGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttgg                                                                                             >  1:283338‑M1/1‑98 (MQ=255)
                                                                                                                                                            gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggcttaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcccccagcaacttggtcctc                                                                                        >  2:680720‑M1/1‑93 (MQ=255)
                                                                                                                                                            gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatc                                                                                        <  1:185050‑M1/251‑159 (MQ=255)
                                                                                                                                                              ttttCTAATTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCGTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgtttcccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatct                                                                                       <  2:618948‑M1/250‑160 (MQ=255)
                                                                                                                                                                     aGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccat                                                                                >  1:110961‑M1/1‑84 (MQ=255)
                                                                                                                                                                               ggtaggtaTAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttgg                                                                     <  2:147061‑M1/251‑178 (MQ=255)
                                                                                                                                                                                        aaaaTGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcc                                                            <  1:480863‑M1/251‑187 (MQ=255)
                                                                                                                                                                                         aactgggATGGAGCCTCTGCTTCTGGCATGTGTAGGTCAGAATGACTCATGATGTGGTCTGCTActcaagggggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggcccatcacgcgttcgcccagcaacttggtcatcccttccctattactgttggtcagaatccc                                                           >  2:619050‑M1/4‑64 (MQ=255)
                                                                                                                                                                                         aaaggggatggAGCCTCTGCTTCTGGCATGTTTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgtttcccggctacggtagctatcccacttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatatcactgttggtcagcatccc                                                           <  2:631837‑M1/247‑188 (MQ=255)
                                                                                                                                                                                         aaaTGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatccc                                                           >  1:584395‑M1/1‑64 (MQ=255)
                                                                                                                                                                                         aaaTGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatccc                                                           >  1:635872‑M1/1‑64 (MQ=255)
                                                                                                                                                                                          aaTGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaagggggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccg                                                          >  1:508222‑M1/1‑63 (MQ=255)
                                                                                                                                                                                            tGGGATGGAGCCTATGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggt                                                        >  2:761691‑M1/1‑61 (MQ=255)
                                                                                                                                                                                                  ggAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcg                                                  >  1:500417‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                   gAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGCCTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgt                                                 >  2:311425‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                   gAGCCTCTGCTTCGGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggaagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgt                                                 >  2:523295‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                    ggcCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTTGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtt                                                <  1:457452‑M1/250‑199 (MQ=255)
                                                                                                                                                                                                       ctctGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttgg                                             >  1:710848‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                           gCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaag                                          <  1:338952‑M1/250‑205 (MQ=255)
                                                                                                                                                                                                           gCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttgg                                             <  2:716413‑M1/247‑202 (MQ=255)
                                                                                                                                                                                                           gCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttgg                                             >  1:716413‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                               cTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaag                                       >  1:676997‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                                 ggCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagc                                    <  1:231929‑M1/250‑211 (MQ=255)
                                                                                                                                                                                                                  gCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcag                                                                 >  2:688678‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                  gCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcag                                                                 <  1:688678‑M1/220‑182 (MQ=255)
                                                                                                                                                                                                                      gagtcgGTCAGAATGACTCATGATGTGGTCTGCTCctcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctctcccagttgaagatcacccacaaactgttaccccggcgcatacggtccatcacgcgttcgcccaccaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgat                              >  2:220316‑M1/3‑35 (MQ=255)
                                                                                                                                                                                                                           gtgcagaATGACTCATGATGTGGGCTGCTActcaaggtggacttttttcgcaggataccgccagaaacggggattaatcctctttacctgttacccggttacggtagctatcccagttgaagtgcacccccaaactgtttcccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccagattcctgttggtcatcatcccgggtgggcgtttggaataagagcggcgatctac                          <  2:471207‑M1/247‑221 (MQ=255)
                                                                                                                                                                                                                              cAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccatttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatc                      <  2:427523‑M1/251‑225 (MQ=255)
                                                                                                                                                                                                                                gAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgat                       >  1:349707‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                gAATGACTCATGAGGTGGTCTGCTActcnaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgat                       <  2:349707‑M1/248‑224 (MQ=255)
                                                                                                                                                                                                                                    gACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggtt                 <  1:761691‑M1/250‑230 (MQ=255)
                                                                                                                                                                                                                                       tCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatatta                                                                            >  1:99931‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                       tCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatatta                                                                            <  2:99931‑M1/188‑171 (MQ=255)
                                                                                                                                                                                                                                        cATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatgctctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggttgat              <  1:366055‑M1/249‑233 (MQ=255)
                                                                                                                                                                                                                                           gatgTGGTCTGCTActcaaggtggaattgggttgcagtataccggcagaaacggggattaatactctttacctggtacccggctacggtagctatcccagttgaagatcccccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcctcccggttgggcgtttggaagaagatcggcgatcgacgagctggttgatgatc          <  2:323187‑M1/250‑237 (MQ=255)
                                                                                                                                                                                                                                             tgtgGTCTGCTActaaaggtggaattgtgtcgcagtataccgcccgaaacggggattaatactctttacctgttacccggctacggtagctatcccagcttaagataacccacaaactgttacccaggcgcattcggtccatcacgcgttagcccagcaacttggtcatctcttacatattactgttggtcagcatcccggttgggcgtttggaagaagngcggcgatcgccgctctggttgatgctcact       >  2:351053‑M1/1‑12 (MQ=255)
                                                                                                                                                                                                                                              gtgGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggttgatgatcactt      >  1:731829‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                                 gTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggttgatgatcacttttt   >  2:100268‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                                 gTCTGCTActcaaggtggaattgtgtcgcagaataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccgtcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctgtttgatgatcacttttt   >  1:556613‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                                  tCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggttgatgatcactttttc  >  1:49664‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                                     gCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatattactgttggtcagcatcccggttgggcgtttggaagaagagcggcgatcgacgatctggttgatgatcactttttc  <  1:680720‑M1/248‑245 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCG  >  NC_000913/4518185‑4518677

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