Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,086,996 Δ13,531 bp yiiGrhaT yiiG, frvR, frvX, frvB, frvA, rhaM, rhaD, rhaA, rhaB, rhaS, rhaR, rhaT

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4086996 4100526 13531 44 [0] [1] 46 yiiG–rhaT yiiG,frvR,frvX,frvB,frvA,rhaM,rhaD,rhaA,rhaB,rhaS,rhaR,rhaT

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40869950 (0.000)44 (0.680) 40/476 0.2 100% intergenic (+146/‑7) fdhD/yiiG sulfurtransferase for molybdenum cofactor sulfuration/DUF3829 domain‑containing lipoprotein YiiG
?NC_000913 4100527 = 0 (0.000)intergenic (‑2/‑283) rhaT/sodA rhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)

CGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTTTTCGACGAAAAAATGGATG  >  NC_000913/4086746‑4087239
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cGCTTGCTGTAGAAGGGGGTGGGCGCTGTAATCTGACACTGGTAGGGTTTTGTAAACCGGTTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTGTTTCATTTATTAGGTTCACTCTCTTATCTATATATTTATAAGGCt                                                                                                                                                                                                                                                     <  2:311210‑M1/251‑2 (MQ=255)
                   tGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcattttt                                                                                                                                                                                                                                  >  1:332852‑M1/1‑231 (MQ=255)
                       cgcTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtc                                                                                                                                                                                                                              >  2:213033‑M1/1‑227 (MQ=255)
                             aaTCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtc                                                                                                                                                                                                                         >  2:202465‑M1/1‑221 (MQ=255)
                               tCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtctcgtcac                                                                                                                                                                                                                       >  2:51824‑M1/1‑219 (MQ=255)
                                                           ggTCGGGCAACGCTTTATACCCATCCGCAGCGTGTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACATATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTATTGGCGTTTTCTTTTTTCATTTATTATGTTCTCTCGCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcttcgcgtaccatcgccgacgagatgagtgtaaaaatcgt                                                                                                                                                                                          <  1:279889‑M1/251‑61 (MQ=255)
                                                            gTCGCGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtg                                                                                                                                                                                         >  1:234025‑M1/1‑190 (MQ=255)
                                                                  cAACGGTTTATACCCATCCGCAGCGTTTGAGCATTTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcatcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcg                                                                                                                                                                                   >  1:230889‑M1/1‑184 (MQ=255)
                                                                   aaCGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcg                                                                                                                                                                                   >  2:69532‑M1/1‑183 (MQ=255)
                                                                                 aTCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcct                                                                                                                                                                    <  1:453448‑M1/251‑83 (MQ=255)
                                                                                       aGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgcc                                                                                                                                                              <  2:98629‑M1/251‑89 (MQ=255)
                                                                                                gCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtc                                                                                                                                                     <  1:184422‑M1/251‑98 (MQ=255)
                                                                                                          aaaaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTTTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgttgtcgattaacctttcgcctgttgccgccgttgtcgatttactgg                                                                                                                                           <  2:465438‑M1/251‑108 (MQ=255)
                                                                                                            aaTAAGTCCAGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttngcctgttgccgccgttgtcgatttactggc                                                                                                                                          >  2:250985‑M1/1‑142 (MQ=255)
                                                                                                                     tGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggc                                                                                                                                 >  1:155368‑M1/1‑133 (MQ=255)
                                                                                                                                 tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttg                                                                                                                     <  2:337251‑M1/250‑130 (MQ=255)
                                                                                                                                  aaaTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCTTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaaccgtgcgccggttgccgccgttgtcgatttactggcaatcacggcattaagtggttgat                                                                                                                   <  2:203360‑M1/251‑132 (MQ=255)
                                                                                                                                      tGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgcttt                                                                                                                                                >  2:19271‑M1/1‑116 (MQ=255)
                                                                                                                                      tGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgattt                                                                                                                                                <  1:19271‑M1/218‑103 (MQ=255)
                                                                                                                                           gACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcac                                                                                                          >  1:459822‑M1/1‑111 (MQ=255)
                                                                                                                                           gACTTAATCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcac                                                                                                          <  1:497969‑M1/251‑141 (MQ=255)
                                                                                                                                                   catgcgttgaattttATTATTTTATTTATAANAAAAATGCGTGTAGTGGCGTTTTTTTTTTGCATTTATTATGTTCCCTCTCTTATTGATATATTTATAAGGCgtcttagcctcatcatttttcgtcgcgtcacatctccggcgagatgcgtgtaaaaatcgtgctgtcgattaacctttctcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcagg                                                                                                  <  2:240284‑M1/244‑149 (MQ=255)
                                                                                                                                                        cGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGTCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatcccgggcattt                                                                                             >  2:13581‑M1/1‑98 (MQ=255)
                                                                                                                                                              ttttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCTCTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctcc                                                                                                                                                                                                                 <  2:259533‑M1/129‑38 (MQ=255)
                                                                                                                                                              ttttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgc                                                                                       >  1:459191‑M1/1‑92 (MQ=255)
                                                                                                                                                              ttttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTCAATCTCTTATCTATATATTTATAAGGCttattatcctcatcatttttcgtcgcgtcacatctcc                                                                                                                                                                                                                 >  1:259533‑M1/1‑92 (MQ=255)
                                                                                                                                                                 tattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaa                                                                                     <  1:450899‑M1/250‑162 (MQ=255)
                                                                                                                                                                 tattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgc                                                                                                                                                               >  2:26786‑M1/1‑89 (MQ=255)
                                                                                                                                                                    tattttatttATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaacc                                                                                 >  1:259529‑M1/1‑86 (MQ=255)
                                                                                                                                                                     attttatttATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaacc                                                                                 >  2:430379‑M1/1‑85 (MQ=255)
                                                                                                                                                                             tataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccataccctt                                                                         <  2:497811‑M1/250‑174 (MQ=255)
                                                                                                                                                                              ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttaatggcaatcacggcattaagtggttgatttgcttcccatctcgggcattttcctgcaaaaccatacccttac                                                                       >  2:488205‑M1/1‑76 (MQ=255)
                                                                                                                                                                                     aaaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgcaaagt                                                                >  2:203590‑M1/1‑69 (MQ=255)
                                                                                                                                                                                      aaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagt                                                                >  1:427920‑M1/1‑68 (MQ=255)
                                                                                                                                                                                        tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacg                                                             >  1:13548‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                         tCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattt                                             <  2:155368‑M1/250‑202 (MQ=255)
                                                                                                                                                                                                               ttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcggcgcgtcacatctccgacgagaggagtgtaaaaagcgtgctgtcgattaacctttcgccttttgccgccgttgtcgatttactggctatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaata                                       <  1:213033‑M1/250‑208 (MQ=255)
                                                                                                                                                                                                                               aCTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactat                       >  1:335060‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                           atatatTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataatgaaccaact           <  1:113001‑M1/250‑236 (MQ=255)
                                                                                                                                                                                                                                             atatTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcagttaatcaactataatgaaccaactgc         <  1:94669‑M1/250‑238 (MQ=255)
                                                                                                                                                                                                                                               atTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaattacggcattgataatcattttcaatatcatttaattaactataatgaaccaactgctt       >  2:471710‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                               atTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaaccttgcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattgtcaatatcatttaattaactataattaaccaactgctt       <  1:149003‑M1/250‑240 (MQ=255)
                                                                                                                                                                                                                                                  tatccggcttcttatcctcatcatttttcgtcgcgtcacatctccgccgagatgagtgtaaaaatcgtgctgtcgcttaacctttcgcctgttgccgccgttgtagatttactggcaatcacggcattaagtcgttgatttgcttcacatctcgggcattttcctgcacaaccatacccttacgcaaagtacggcattgataatcattttcaatatcatatatttaactataatgaaccagctgattacgc   >  2:186672‑M1/6‑8 (MQ=255)
                                                                                                                                                                                                                                                   ataAGGCttcttagcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtggttgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataatgaaccaactgcttacgcg  <  2:283453‑M1/251‑245 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
CGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTTTTCGACGAAAAAATGGATG  >  NC_000913/4086746‑4087239

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