Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,299,522 Δ12,535 bp mdtPalsB 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 43 [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)43 (0.720) 40/494 0.3 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

GCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCG  >  NC_000913/4299275‑4299764
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gCCTCCGCCCGGACGATTTTGACGTTATCGACGACCAGGTTGCTTTTGCAACCTGATGCGCAATAGAGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCCCTCGACTACAGTATGCATCTTTTATGCCACATTTTCTGTGGGGTCCTTCCCTGATATTACGGGCACTATTTATTCAAAAATCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTTAggt                                                                                                                                                                                                                                                  >  2:444350‑M1/1‑247 (MQ=255)
     cGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaagggg                                                                                                                                                                                                                                            <  2:7402‑M1/251‑10 (MQ=255)
          ggACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagt                                                                                                                                                                                                                                       >  2:513537‑M1/1‑237 (MQ=255)
          ggACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagt                                                                                                                                                                                                                                       <  1:409001‑M1/251‑15 (MQ=255)
          ggACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagt                                                                                                                                                                                                                                       <  1:262767‑M1/251‑15 (MQ=255)
              gATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGGAATAGGGGCAAACCGTGACCACTTTTTTCATCGGTTGCGCTCCAGTTAATCAAATCACTCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACAATTTATGTGGGGTCATTCCCTGAGATGACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTGAGTTTTGACGTGAggtaaggggacagtaact                                                                                                                                                                                                                                   <  2:283478‑M1/251‑19 (MQ=255)
                          aTCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcg                                                                                                                                                                                                                        <  1:128193‑M1/250‑30 (MQ=255)
                                  ccAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacct                                                                                                                                                                                                               <  2:264651‑M1/251‑39 (MQ=255)
                                    aGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattaccttt                                                                                                                                                                                                             >  1:102034‑M1/1‑211 (MQ=255)
                                    aGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattaccttt                                                                                                                                                                                                             <  1:415155‑M1/251‑41 (MQ=255)
                                         gATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttctt                                                                                                                                                                                                        <  1:101044‑M1/251‑46 (MQ=255)
                                            tttGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttat                                                                                                                                                                                                      <  2:479212‑M1/250‑48 (MQ=255)
                                               gCAACCGGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttca                                                                                                                                                                                                  <  2:262357‑M1/251‑52 (MQ=255)
                                                           gcAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtg                                                                                                                                                                                      >  2:115988‑M1/1‑188 (MQ=255)
                                                             aaTAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgag                                                                                                                                                                                    <  1:23348‑M1/251‑66 (MQ=255)
                                                                  gggCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaa                                                                                                                                                                               <  2:337090‑M1/251‑71 (MQ=255)
                                                                     cAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatc                                                                                                                                                                            <  1:277914‑M1/251‑74 (MQ=255)
                                                                                      tttCAGCGGTCTCGCTCCAGTTAAGCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggtt                                                                                                                                                            <  2:30985‑M1/250‑90 (MQ=255)
                                                                                                  cGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatat                                                                                                                                                  <  2:421510‑M1/248‑100 (MQ=255)
                                                                                                         gTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagc                                                                                                                                        >  2:468013‑M1/1‑142 (MQ=255)
                                                                                                         gTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagc                                                                                                                                        <  1:390236‑M1/251‑110 (MQ=255)
                                                                                                            aatcaaatcaCGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagctt                                                                                                                                      <  1:468013‑M1/250‑112 (MQ=255)
                                                                                                                    caCGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgt                                                                                                                             >  2:115522‑M1/1‑131 (MQ=255)
                                                                                                                         aTACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagc                                                                                                                        >  1:366718‑M1/1‑126 (MQ=255)
                                                                                                                              cgcTCTCGACTACAGTATGCATCTTTTATGCCGCAGTTTATGTGTGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggc                                                                                                                   <  2:520418‑M1/251‑131 (MQ=255)
                                                                                                                                           aGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttcagctgggcgacagagacaaa                                                                                                       >  2:177121‑M1/1‑108 (MQ=255)
                                                                                                                                                              acaTTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttg                                                                                   >  1:270548‑M1/1‑89 (MQ=255)
                                                                                                                                                                              cATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttg                                                                    <  1:443744‑M1/250‑178 (MQ=255)
                                                                                                                                                                                 tCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccg                                                                 <  1:115988‑M1/250‑181 (MQ=255)
                                                                                                                                                                                  cccTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccg                                                                 >  1:92791‑M1/1‑69 (MQ=255)
                                                                                                                                                                                         atTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattac                                                        <  2:500461‑M1/251‑190 (MQ=255)
                                                                                                                                                                                          ttACGGTCACTATTTAGTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAgttaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggtttttcccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattac                                                        <  2:472866‑M1/250‑190 (MQ=255)
                                                                                                                                                                                              gggCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacct                                                   <  1:186023‑M1/251‑195 (MQ=255)
                                                                                                                                                                                                     atttattCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataac                                             <  2:211273‑M1/250‑201 (MQ=255)
                                                                                                                                                                                                        tattCAAAACTCTGACGCAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcg                                         >  2:261204‑M1/1‑47 (MQ=255)
                                                                                                                                                                                                        tattCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggc                                          <  1:81307‑M1/250‑204 (MQ=255)
                                                                                                                                                                                                         attCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgt                                        >  1:326385‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                         attCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgt                                        <  1:115522‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                            cAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttc                                     >  1:236612‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                                          aaaaCAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagc                       >  2:180139‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                             aCAGGCTGTCGTCAGTTTTGACGTGAgggtaggggacagaaaatccgcgatggcggacttacatttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggatggtagtaatattcttgccgaggaattacgacctaaaaacggcgtttccgggtctgggtcgcaa                     >  2:180821‑M1/1‑26 (MQ=255)
                                                                                                                                                                                                                                ggCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgctggcggcattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttctatattactattagcttgttcaatgttgagctgggcgacagagacacaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaacta                  >  2:207946‑M1/1‑23 (MQ=255)
                                                                                                                                                                                                                                               tGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaaataatatgatcggccagtt  >  1:150263‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GCCTCCGCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCG  >  NC_000913/4299275‑4299764

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