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 24 [0] [0] 24 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)24 (0.520) 23/482 0.8 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)

GAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTT  >  NC_000913/4086757‑4087220
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
gAAGTGGCTGAGCGCTGTAATCTGGCACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCGGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTGCTTTTTTCATTTATTATTTTCGCGCTCTTATCTATATATTTGTAAGGCttctttgcctca                                                                                                                                                                                                                       <  2:443897‑M1/251‑13 (MQ=255)
                         cacTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctcc                                                                                                                                                                                              <  2:80861‑M1/251‑38 (MQ=255)
                                             ccGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTAATTATAAACAAAATGCGTGTAGTGGCGTTTTCTTGTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttctcgtcgcgacacatctccgacgagatgagtgtaaaaat                                                                                                                                                                          <  1:420039‑M1/251‑58 (MQ=255)
                                                          cGGTTTATACCCATCCGCAGCGGTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgatt                                                                                                                                                             <  2:17400‑M1/251‑71 (MQ=255)
                                                                                    aGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgt                                                                                                                                   <  1:130605‑M1/251‑97 (MQ=255)
                                                                                        attaattaAAATAAGTCCTTCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcttcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaactcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgat                                                                                                                               >  2:324838‑M1/1‑151 (MQ=255)
                                                                                              taaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTATCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcccatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactg                                                                                                                         >  2:434552‑M1/1‑145 (MQ=255)
                                                                                                  aTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATCCATTTCGCGACTCTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTTTATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggca                                                                                                                      >  1:441419‑M1/1‑141 (MQ=255)
                                                                                                                    atTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggt                                                                                                   >  1:12056‑M1/1‑123 (MQ=255)
                                                                                                                                               tGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttc                                                                        >  2:94792‑M1/1‑96 (MQ=255)
                                                                                                                                                  attttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgattt                                                                                              >  1:101554‑M1/1‑93 (MQ=255)
                                                                                                                                                  attttattagtttatgTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgttgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgattt                                                                                              <  2:101554‑M1/226‑134 (MQ=255)
                                                                                                                                                               atttatAAACAAAATGGGGGTAGTGGCGTTTTCTTTTGTCATGTATTATGTTCACTCTCTTATCGATATATTTATAAGGCttcttagcctcatcattgttcgtcgggtgacatctccgacgagatgagtgtaaaaagcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgttttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccataccc                                                        <  1:447911‑M1/251‑172 (MQ=255)
                                                                                                                                                                 ttatAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATCTTCACCCTCTTATCTATATATTTATAACGCttcttatcctcatcctttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttagcctgttgccgccgttgtcgctttactggcaatcacggcattaagtgggtgatttgcttcccatctcgggcattttcctgcaaaaccataccctt                                                      >  1:232574‑M1/1‑78 (MQ=255)
                                                                                                                                                                   ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttac                                                    >  2:365468‑M1/1‑76 (MQ=255)
                                                                                                                                                                       aCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaa                                                 >  2:134290‑M1/1‑72 (MQ=255)
                                                                                                                                                                              gCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacgg                                         >  1:6337‑M1/1‑65 (MQ=255)
                                                                                                                                                                                gtgtAGTGGCGTTTTCTTTTTTCATTTATTATCTTCAATCTCTTATCTATATATTTATCACGCttcttatcctaattctttttcgtcgagtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgactgttgccgccgttgtcgctttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggc                                        >  2:212783‑M1/1‑63 (MQ=255)
                                                                                                                                                                                                 tttttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcatttt                         <  2:384611‑M1/248‑203 (MQ=255)
                                                                                                                                                                                                    ttCATTTCTTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatat                   >  1:158595‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                              atGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaatt          >  1:322358‑M1/1‑33 (MQ=255)
                                                                                                                                                                                                              aggttCACTCTCTTATCTTAATATTTATAAGTCttcgtatcctcatcatttttcgtcgcttcacatctccgacgagatgagtgtaaaaatcgggctgtcgattaaccgttcgcctgttgctgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaatt          <  2:359832‑M1/248‑218 (MQ=255)
                                                                                                                                                                                                                    cctctcttATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactata   <  2:6337‑M1/250‑225 (MQ=255)
                                                                                                                                                                                                                       ctctTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtggaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggtattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataa  <  1:324838‑M1/249‑226 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
GAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTT  >  NC_000913/4086757‑4087220

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