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 76 [0] [2] 77 [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)75 (0.980) 68/490 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

TATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATT  >  NC_000913/4599948‑4600443
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
aaagcgactcgaacgttaaggaagctgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTc                                                                                                                                                                                                                                                        >  2:283924‑M2/248‑250 (MQ=255)
 aagcgactcgaacgttaaggaagctgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacttagtcatattaatatgagaatggttatcattacaattggaaataaaattgttgccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCaa                                                                                                                                                                                                                                                      <  2:593458‑M2/5‑1 (MQ=255)
    cgactcgaacgttaaggaagctgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGt                                                                                                                                                                                                                                                   <  2:65130‑M2/8‑1 (MQ=255)
                 aaggaagctgtggtcataccgcggacagacatgcggcgagcatcaaaagttaattatcacgtagtcctattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggaataaaatgggctggagcctctgcttctggcatgtgtcgggcagaatgacccatgatgtggtttaataCTAAAGGAGGAATTGTGTCGc                                                                                                                                                                                                                                      >  2:511061‑M2/231‑251 (MQ=255)
                 aaggaagctgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCg                                                                                                                                                                                                                                       <  2:527553‑M2/20‑1 (MQ=255)
                       gctgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGtata                                                                                                                                                                                                                                <  1:209007‑M2/27‑1 (MQ=255)
                         tgtggtagtggcgcggtcagacatgcggagtgcatcaaaagttcattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcctgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATAcc                                                                                                                                                                                                                              >  2:349193‑M2/223‑251 (MQ=255)
                          gtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCg                                                                                                                                                                                                                             <  1:183439‑M2/30‑1 (MQ=255)
                          gtggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCg                                                                                                                                                                                                                             >  2:231475‑M2/222‑251 (MQ=255)
                            ggtagtggcgcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGcc                                                                                                                                                                                                                           <  1:238054‑M2/32‑1 (MQ=255)
                                     gcggtcagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACggg                                                                                                                                                                                                                  >  2:407805‑M2/211‑251 (MQ=255)
                                          cagacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtctgtcagaatgactcatgatgtggtctgctaCTCAAGGAGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGAtt                                                                                                                                                                                                              >  1:406902‑M2/206‑250 (MQ=255)
                                            gacatgcggagtgcatcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                           >  1:217754‑M2/204‑251 (MQ=255)
                                                          atcaaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGCTTAATACTCTTTACCTGtt                                                                                                                                                                                             >  2:306775‑M2/190‑251 (MQ=255)
                                                            caaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTa                                                                                                                                                                                            <  1:314990‑M2/63‑1 (MQ=255)
                                                             aaaagttaattatcacgtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAcc                                                                                                                                                                                          >  1:516272‑M2/187‑251 (MQ=255)
                                                                             gtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacttgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgacgcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAt                                                                                                                                                                          <  1:569132‑M2/81‑1 (MQ=255)
                                                                             gtagtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAt                                                                                                                                                                          >  1:525570‑M2/171‑251 (MQ=255)
                                                                               agtcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATcc                                                                                                                                                                        >  1:173305‑M2/169‑251 (MQ=255)
                                                                               agtcatattaatatgagaatggttatcattacaagtggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcaggtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATAACGCCAGAAACGGGGATTAATACTCTTTACCTGGTACCCGGCTACGGTAGctaac                                                                                                                                                                         <  2:18340‑M2/82‑3 (MQ=255)
                                                                                 tcatattaatatgagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                                                       >  1:630003‑M2/167‑250 (MQ=255)
                                                                                  catattaatatgagaatggttatcattacaattggaaataaaagtgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACTCGGCTACGGTAGCTATCCCAg                                                                                                                                                                     <  2:251475‑M2/86‑1 (MQ=255)
                                                                                             gagaatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtgtttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACGCTTTACCTGTTGCCCGGCTACGGTAGCTATCCCAGTTGAAGATCAc                                                                                                                                                          <  1:139633‑M2/97‑1 (MQ=255)
                                                                                             gagaatgggtatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcttctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCAAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGCAGATCAc                                                                                                                                                          >  2:270962‑M2/155‑251 (MQ=255)
                                                                                               gaatggatatcattccaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAccc                                                                                                                                                        >  2:463410‑M2/153‑251 (MQ=255)
                                                                                                aatggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgattttataagtgttataaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTAAAGATCAccc                                                                                                                                                        >  1:540589‑M2/152‑250 (MQ=255)
                                                                                                  gggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggctgtagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGACTTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAtcaccaac                                                                                                                                                      >  2:2267‑M2/150‑247 (MQ=255)
                                                                                                   ggttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcaca                                                                                                                                                     >  1:88499‑M2/149‑250 (MQ=255)
                                                                                                    gttatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACaa                                                                                                                                                    >  1:196991‑M2/148‑250 (MQ=255)
                                                                                                      tatcattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACt                                                                                                                                                 >  1:21235‑M2/146‑251 (MQ=255)
                                                                                                         cattacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGtt                                                                                                                                              >  2:334188‑M2/143‑251 (MQ=255)
                                                                                                           ttacaattggaaataaaattgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAc                                                                                                                                            >  1:546789‑M2/141‑251 (MQ=255)
                                                                                                                            attgtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTcc                                                                                                                           >  1:266194‑M2/124‑251 (MQ=255)
                                                                                                                               gtttccaatagacatttttaacatgttgtttttctaagtgttgtaatgtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTTCCTGTTTCCCGTCTACGGTAGCTATCCCAGTTGAAGATCACTCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                                                         <  2:502557‑M2/130‑1 (MQ=255)
                                                                                                                               gtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATc                                                                                                                        >  1:300481‑M2/121‑251 (MQ=255)
                                                                                                                               gtttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATAACGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATc                                                                                                                        >  2:241607‑M2/121‑251 (MQ=255)
                                                                                                                                tttccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATc                                                                                                                        >  1:590075‑M2/120‑250 (MQ=255)
                                                                                                                                  tccaatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAc                                                                                                                      <  1:407805‑M2/133‑1 (MQ=255)
                                                                                                                                    caatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCCTCAcgcg                                                                                                                   >  2:311687‑M2/116‑251 (MQ=255)
                                                                                                                                     aatagacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGt                                                                                                                  >  1:400377‑M2/115‑251 (MQ=255)
                                                                                                                                        agacatttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCg                                                                                                               <  1:311687‑M2/140‑1 (MQ=255)
                                                                                                                                           catttttaacatgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGccc                                                                                                            >  1:591759‑M2/109‑251 (MQ=255)
                                                                                                                                                  aacatgttgtttttctaagagttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGCATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTACATACCGCGTTCGCCCagcacat                                                                                                     >  2:451228‑M2/102‑248 (MQ=255)
                                                                                                                                                     atgttgtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGGGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACATgg                                                                                                  >  2:20365‑M2/99‑251 (MQ=255)
                                                                                                                                                         tgtttttctaagtgttgtaaggtatgtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCAAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAt                                                                                              >  1:227389‑M2/95‑251 (MQ=255)
                                                                                                                                                          gtttttctaagtgttgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtc                                                                                             <  2:196991‑M2/158‑1 (MQ=255)
                                                                                                                                                                      ggtgtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATAtt                                                                                  <  2:516272‑M2/169‑1 (MQ=255)
                                                                                                                                                                        tgtaaggtcggtataaaatgggatggcgcctctgcttctggcatgtgtaggtcagaatgactgatgctgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGATACGGGGATTAATAATCTTTACCTGTTACACGGCTCCGGTAGCTATCCCAGTTGAAGATCAACCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATCTTACt                                                                               >  2:546599‑M2/80‑251 (MQ=255)
                                                                                                                                                                         gtaaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACt                                                                               <  1:2267‑M2/172‑1 (MQ=255)
                                                                                                                                                                           aaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGtt                                                                            <  2:217754‑M2/175‑1 (MQ=255)
                                                                                                                                                                           aaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGtt                                                                            >  2:389779‑M2/77‑251 (MQ=255)
                                                                                                                                                                           aaggtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGtt                                                                            <  2:525570‑M2/175‑1 (MQ=255)
                                                                                                                                                                            aggcaggtatcaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTg                                                                           >  2:278795‑M2/76‑251 (MQ=255)
                                                                                                                                                                              gtaggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGt                                                                         >  1:19920‑M2/74‑251 (MQ=255)
                                                                                                                                                                                aggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCa                                                                       <  2:106607‑M2/180‑1 (MQ=255)
                                                                                                                                                                                 ggtataaaatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAg                                                                      <  2:458583‑M2/181‑1 (MQ=255)
                                                                                                                                                                                        aatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAgg                                                                                                                                       <  2:286069‑M2/116‑1 (MQ=255)
                                                                                                                                                                                        aatgggatggagcctctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAgg                                                                                                                                       >  1:286069‑M2/64‑179 (MQ=255)
                                                                                                                                                                                                      tctgcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTgg                                                  >  1:277279‑M2/50‑250 (MQ=255)
                                                                                                                                                                                                         gcttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCCCGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGga                                                 >  1:205041‑M2/47‑248 (MQ=255)
                                                                                                                                                                                                           ttctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                              >  2:526672‑M2/45‑249 (MQ=255)
                                                                                                                                                                                                             ctggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGaagaag                                            >  1:488534‑M2/43‑249 (MQ=255)
                                                                                                                                                                                                              tggcatgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCCTCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGaagagg                                         >  2:633565‑M2/42‑250 (MQ=255)
                                                                                                                                                                                                                  atgtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                     >  2:518717‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                                  atgtgtcggtaagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                     >  1:589900‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                                    gtgtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGAt                                   >  1:255519‑M2/36‑251 (MQ=255)
                                                                                                                                                                                                                      gtcggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATc                                  <  1:283924‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                                         ggtcagaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcgacg                              >  2:518565‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                         ggtcagaatgaatcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATTCCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcgacg                              <  2:432220‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                                              gaatgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTg                         <  1:334188‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                                atgactcatgatgtggtcggctaCTCAAGGTGGAATTTTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTGCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACACGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGAGATCGACGATCTGGt                       <  2:406902‑M2/228‑1 (MQ=255)
                                                                                                                                                                                                                                 tgactcatgatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGtt                      >  1:473048‑M2/23‑251 (MQ=255)
                                                                                                                                                                                                                                         gatgtggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATc               >  1:278116‑M2/15‑250 (MQ=255)
                                                                                                                                                                                                                                            gtggtctgctACTCCAGGTGTAATTGTGTCGCATTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTGACCTGGCTACCGTAGCTGTCCAAGTTGAATATCACCCACAAACTGTTCCCCAGGCTCATACGTTCCATCACGCGTTCGCCCAGCAACTTTTTCATCTCTTCCATATTACTGTTGCTCATCATCCCGGTTGGACGTTTGGAATAAGAGCGGCTATCTACGATCTATTTGTTGATCACtt           >  2:68962/10‑251 (MQ=255)
                                                                                                                                                                                                                                              ggtctgctaCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACtttt         >  1:611130‑M2/10‑251 (MQ=255)
                                                                                                                                                                                                                                                   gctaCTCAAGGTGGAATTGGGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAATCTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTCCGTa    <  2:18936‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                                                     tACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAtt  >  1:180580/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
TATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATT  >  NC_000913/4599948‑4600443

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