Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,299,522 Δ12,535 bp mdtP–alsB mdtP, mdtO, mdtN, ytcA, yjcS, alsK, alsE, alsC, alsA, alsB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4299522 4312056 12535 41 [0] [0] 41 mdtP–alsB mdtP,mdtO,mdtN,ytcA,yjcS,alsK,alsE,alsC,alsA,alsB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42995210 (0.000)41 (0.630) 35/490 0.5 100% intergenic (‑155/+43) fdhF/mdtP formate dehydrogenase H/putative multidrug efflux pump outer membrane channel
?NC_000913 4312057 = 0 (0.000)intergenic (‑15/+44) alsB/alsR D‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR

ATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGT  >  NC_000913/4311835‑4312299
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
gatcgacgaccaggttgattttgcaacctgatgcgcaatagggtcaaaccgtgacgtcttttttcatcgctctcgctccagttaatcaaatcacgcatacgcgctctcgactacaatatgcatcttttattccacattttatgtggggtcattccctgatattacgggctctatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGc                                                                                                                                                                                                                        <  2:169014‑M2/28‑1 (MQ=255)
  tcgacgaccaggttgattttgcaaactgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGcggg                                                                                                                                                                                                                     >  2:204527‑M2/221‑250 (MQ=255)
     ccgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAAtt                                                                                                                                                                                                                  <  1:417750‑M2/34‑1 (MQ=255)
      cgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAAtt                                                                                                                                                                                                                  >  1:177416‑M2/217‑250 (MQ=255)
       gaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTAc                                                                                                                                                                                                                >  2:26511‑M2/216‑251 (MQ=255)
           aggttgcttttgcaacctgatgcgcaatcggggcaaaccgtgacgacttttttcatcgggctcgctccatttaatcaaatcacgcatacgcgctctcgattacagtatgcatcttttatgccacattttatgttgggtcattccctgatattacgggcactatttattcaaaactcggacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGTACA‑GTAACTCCGCGATGGCGGAAGTACCttt                                                                                                                                                                                                            <  2:466873‑M2/40‑1 (MQ=255)
               tgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                                                                                                         <  2:330775‑M2/43‑1 (MQ=255)
               tgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                                                                                                         <  1:392364‑M2/43‑1 (MQ=255)
                gattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                                                                                                                                       <  2:218673‑M2/45‑1 (MQ=255)
                    ttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGCAATTACCTTTTTCTTAtt                                                                                                                                                                                                    <  2:446103‑M2/48‑1 (MQ=255)
                      gcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                                                                                                                                                  >  2:215848‑M2/201‑250 (MQ=255)
                        aacctgatgcgcaatcgcggcaaaccgtgacgacttttttcatcgctctcgctccagctaatcaaatcacgcatacgcgatctcgactacagtatgcatcttttatgccacattttatgtggggtaataccctgatattacgggccctatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGCAAGGGGACA‑GTAACTCCGCGATGGGGGAATTACCTTTTTCTTAtttcggc                                                                                                                                                                                               >  2:35807‑M2/199‑248 (MQ=255)
                          cctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCagcgt                                                                                                                                                                                             <  2:197353‑M2/55‑2 (MQ=255)
                               ggcgcaataggtgcaaaccgttgcgacttttttcagcggtctagatgcagttaatcaaatcacgcatacgagctctcgactacagtatgcatcttttatgccacattttaggtgtggttattcccttgtattacgggcactatttattcaaaactctgactaaaaacaggctgtcgtcagttttgtcgtgcGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTtcagcgcgtcg                                                                                                                                                                                         <  2:282981‑M2/59‑4 (MQ=255)
                                         gggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgcGGGAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTCTCTTATTTCAGCGCGACTGGTGagaaga                                                                                                                                                                               >  1:316487‑M2/182‑250 (MQ=255)
                                                 cgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                                                                                                                                      <  2:114699‑M2/78‑1 (MQ=255)
                                                              ttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                                                                                                                                          <  1:269929‑M2/90‑1 (MQ=255)
                                                                catcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                                                                                                       <  1:98021‑M2/93‑1 (MQ=255)
                                                                           ctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                                                                                >  1:257019‑M2/148‑215 (MQ=255)
                                                                           ctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                                                                                <  2:257019‑M2/68‑1 (MQ=255)
                                                                                gttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacagtctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAGTGCTCCGGTTTTTTCCATATTAATATTAGc                                                                                                                                       <  1:161804‑M2/109‑1 (MQ=255)
                                                                                 ttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCt                                                                                                                                      >  2:214227‑M2/142‑251 (MQ=255)
                                                                                                  acgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTg                                                                                                                     <  1:49764‑M2/127‑1 (MQ=255)
                                                                                                      gctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctttcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCg                                                                                                                 <  1:420673‑M2/131‑1 (MQ=255)
                                                                                                                         atcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGt                                                                                              >  1:182893‑M2/102‑251 (MQ=255)
                                                                                                                            ttttatgcaaaattttatggggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGcgttc                                                                                            >  1:160039‑M2/99‑248 (MQ=255)
                                                                                                                              ttatgccacattttatgtggggtcattccctgatattacgggcactatttagtcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATAACAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAg                                                                                          <  2:7889‑M2/154‑1 (MQ=255)
                                                                                                                               tatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGc                                                                                         <  1:471827‑M2/155‑1 (MQ=255)
                                                                                                                                      cattttatgtggggtcattccctgatattacgggaactatttcttcaaaactctgacgaaaaacaggctgtcgtcagttttgccgtgtGGTACGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGTCAGAGACAAAAAAAGCgtccagcatggttgg                                                                                 >  2:463859‑M2/89‑244 (MQ=255)
                                                                                                                                       attttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTGCCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGTGCGACAGAGACAAAAAAAGCGTCCAGcaaggttggt                                                                                <  1:434839‑M2/164‑4 (MQ=255)
                                                                                                                                                ggggtaattccctgatattacgggctctagtcattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTGTCTTAGTGCAGCGCGACTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGTGACAGATACaaaaagagcgtccagcaagtgtggttggaatatt                                                                       <  2:60809‑M2/173‑17 (MQ=255)
                                                                                                                                                  ggtcattccctgatattacgggcactatttattcaaatctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATGTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCt                                                                     <  2:412948‑M2/175‑1 (MQ=255)
                                                                                                                                                              atattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCAtt                                                         <  1:114087‑M2/187‑1 (MQ=255)
                                                                                                                                                                attacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTAc                                                       >  1:242439‑M2/63‑251 (MQ=255)
                                                                                                                                                                               tattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGTGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGc                                         <  1:52288‑M2/203‑1 (MQ=255)
                                                                                                                                                                                  tcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGttt                                     >  1:35012‑M2/45‑251 (MQ=255)
                                                                                                                                                                                               gaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTg                         >  1:1493‑M2/32‑250 (MQ=255)
                                                                                                                                                                                                   aacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCa                     >  1:365666‑M2/28‑250 (MQ=255)
                                                                                                                                                                                                     caggctgtcgtcagttttgacgtgaGGTAAAGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaat                  >  2:74434‑M2/26‑251 (MQ=255)
                                                                                                                                                                                                             cgtcagttttgactggaGGTAAGGGGACA‑TTTACTCCGCGATGGATGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTGCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCaaataatatga           <  2:316487‑M2/233‑3 (MQ=255)
                                                                                                                                                                                                                      tgacgggaGGTAAGGGTACA‑GTAACTTCGCGATGGCAGAATTACCTTTTTCTGATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTGCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATAGGATCGGCCAGt  <  2:542805‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
ATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGT  >  NC_000913/4311835‑4312299

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