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 42 [0] [0] 43 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)42 (0.580) 39/492 0.4 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

CTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGC  >  NC_000913/4311814‑4312304
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
ccgcccggacgattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAg                                                                                                                                                                                                                                                    <  1:74735‑M2/6‑1 (MQ=255)
              ttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaatccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcttcagttttgacgtgaGGTAAGGGGACA‑GTAACTCcg                                                                                                                                                                                                                                    <  1:463562‑M2/21‑1 (MQ=255)
                  cgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATg                                                                                                                                                                                                                               <  1:249922‑M2/26‑1 (MQ=255)
                     tatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGc                                                                                                                                                                                                                             <  1:133032‑M2/28‑1 (MQ=255)
                      atcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑gtaactccgcgcgggcgg                                                                                                                                                                                                                           >  2:299621‑M2/222‑244 (MQ=255)
                        cgacgaccaggttgattttgcaccctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaaccaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGa                                                                                                                                                                                                                          <  1:282618‑M2/31‑1 (MQ=255)
                             accaggttgattttgcaacctgatgcgcactaggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTAcc                                                                                                                                                                                                                    >  1:90389‑M2/215‑251 (MQ=255)
                                  gttgattttgcaacctgatgcgcaataggggcaaaccgagacgacatttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GCAACTCCGCGCTGGCGGAATTACCttttt                                                                                                                                                                                                               >  2:569049‑M2/210‑251 (MQ=255)
                                     gattttgcaacctgatgcgtaataggggcaaaccgtgacggcttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                                                                                                                                             <  2:528334‑M2/44‑1 (MQ=255)
                                        tttgcaagctgatgcgcaatagggacaaaccgtgacgacttctttcatcggtctcgctacagttaatcaaatcacgcatacgagctctcgacaacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaacctctgacgaaaaacaggctgtcgtcagttttgacgttaGGGAAGGGG‑CAAGTAACTCCGCGAGGGCGGAAATACCTTTttcttctt                                                                                                                                                                                                         >  2:203574‑M2/204‑248 (MQ=255)
                                               cctgatgcgcaatagggtcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctcgcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaacgctgacgaaaaacaggctgtcgtcagttttggcgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGGGTAATTACCTTTTTCTTGTTTgagctc                                                                                                                                                                                                  <  2:195345‑M2/55‑3 (MQ=255)
                                                          ataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                                                                                                                                          >  1:545307‑M2/186‑248 (MQ=255)
                                                             gggtcaaaccgagacgacttttttcatcggtctcgctccagttaataaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactctttattcaaaactctgccgaaaaacaggctgtcgtcatttttgtcgcgaGGTAAGGGGACA‑GCAACTCCGCGATGGTGGAATTACCTTTTTCTTATTTCAGCGCGCATGGTgaga                                                                                                                                                                                       >  2:314542‑M2/183‑248 (MQ=255)
                                                               ggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcaggtttgacgagaGGTAAGGGGA‑AAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgcagaa                                                                                                                                                                                   >  1:321837‑M2/181‑250 (MQ=255)
                                                                 caaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                                                                                                                                                <  1:399095‑M2/73‑1 (MQ=255)
                                                                      cgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGGAGAAATcaaca                                                                                                                                                                            >  2:234586‑M2/174‑250 (MQ=255)
                                                                             acttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                                                                                                                                                    <  1:299621‑M2/85‑1 (MQ=255)
                                                                                tttttcatcggtctcgctccagttaatcaaatcacgcatacgcgcgctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                                                                                                                                                  <  1:73770‑M2/87‑1 (MQ=255)
                                                                                   ttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                                                                                                                                               <  2:135192‑M2/90‑1 (MQ=255)
                                                                                      atcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGgttttttc                                                                                                                                                           >  1:388653‑M2/158‑250 (MQ=255)
                                                                                              cgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgcctcattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggcggtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTGTTTTCCatattaa                                                                                                                                                   <  2:491251‑M2/102‑1 (MQ=255)
                                                                                                 tccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaatat                                                                                                                                                >  1:229655‑M2/147‑251 (MQ=255)
                                                                                                       taatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCt                                                                                                                                           <  1:311841‑M2/110‑1 (MQ=255)
                                                                                                            aaaacacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCa                                                                                                                                     <  1:477619‑M2/116‑1 (MQ=255)
                                                                                                                                cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAc                                                                                                                    <  1:509218‑M2/133‑1 (MQ=255)
                                                                                                                                cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACag                                                                                                                  <  1:225509‑M2/135‑1 (MQ=255)
                                                                                                                                  actacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAGTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                                                                                               <  2:233833‑M2/138‑1 (MQ=255)
                                                                                                                                                    tatgccacatgttcagtggggtcattccctgatattacgggcactatttattcaaaactctgacgaacaacaggatgtcggcagttttgacgtgaTGTTAGGGGACA‑GTAACTCCGCGATGGCGGAATTACGTTTGTCTTATTTCAGCGCTCCTTGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATAGTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGc                                                                                              <  2:313268‑M2/155‑1 (MQ=255)
                                                                                                                                                        ccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCagcaaggt                                                                                         <  1:553254‑M2/160‑4 (MQ=255)
                                                                                                                                                         cacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGGTCAATGTTTAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAggtt                                                                                        <  2:40502‑M2/161‑1 (MQ=255)
                                                                                                                                                           cattttatgtggggttattccctgatatcacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGACGGCGGAATTACCTTTTTCTTATTTCAGCGCGCATGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCgtccagcacggttgg                                                                                      >  2:544977‑M2/89‑244 (MQ=255)
                                                                                                                                                                 atgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTa                                                                                 >  1:141842‑M2/83‑250 (MQ=255)
                                                                                                                                                                    tggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAatat                                                                             >  1:96955‑M2/80‑251 (MQ=255)
                                                                                                                                                                                    tattacgggcactatttattcaaaactctgacgaaaaacaggatgtcgtcagtcttgacgtgaGGTAAGGGGACA‑GTAACTACGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCCTATTAATATTAGCTTGATCAATGTTGAGCTGGGCGACAGCGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAAtt                                                              >  2:417339‑M2/64‑250 (MQ=255)
                                                                                                                                                                                        acgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGAc                                                         >  1:309764‑M2/60‑251 (MQ=255)
                                                                                                                                                                                                              tctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGc                                   >  1:304933‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                                         acaggctgtcgtctgttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGGATTACCTTTTTCTTATTACAGCGCGCCTGGTGAGAAGAAATAAACAATTGCGCCGGTTTTTTCCATATTATTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTATTGCCGAGGGATTAAGACCTAATAACGGCGTTTCCGGGGctggtgggaaa                         >  2:129‑M2/27‑245 (MQ=255)
                                                                                                                                                                                                                             gctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaataat                    <  1:499150‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                                      agttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc           <  1:125738‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                         tttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGcc          >  2:442872‑M2/11‑249 (MQ=255)
                                                                                                                                                                                                                                           tgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGtt      >  1:20637‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                             acgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCGTGGTGATAATAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCTTCCAGCAAGGTTGGTTGTAATATTCTTGGCGTTGCATTGCGTCCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc    <  1:32064‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                               gtgaggTAAGGGGACAG‑TAACTCGGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGGGCGAAGAAATCAACAATTGCGACGGTTTTTTCCGTATTAATATTAGATTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCAAGCAAGGTTGTTTGTAATATTCCTGCCGAGGCAATACGACCTAATAACTGCGTTTGCGGGGCTGGTGAGCAAATAAAATGATCGGGCCGTTTCGc  >  2:571413/5‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
CTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGC  >  NC_000913/4311814‑4312304

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: ‑