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 19 [0] [0] 20 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)19 (0.390) 19/488 1.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)

AAACAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTA  >  NC_000913/4100297‑4100776
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
gagcgctgtaatctgacactggtaggtttttgtaaaccgggtcgggcaacggtttatacccatccgcagcgtttgagcaattaattaaaataagtcctgcgaaatatattaaattgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTtc                                                                                                                                                                                                                                       <  1:217275‑M2/21‑1 (MQ=255)
           tctgacactggtaggtttttgtaaaccgggtctggcaacggtttatacccatccgcagcgtttgagcaattaattaaaataagtcctgctaaatatattaaattgcaggacttattcatttcgtgcattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataagtcTTCTTATCCTCATCATTTTTCGTCTCGtcacc                                                                                                                                                                                                                            >  2:151701‑M2/220‑250 (MQ=255)
                   tggtaggtttttgtaaaccgggtcgggcaacggtttatacccatccgcagcgtttgagcaattaattaaaataagtcctgcgaaatatattaaattgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCcgac                                                                                                                                                                                                                    >  1:49296‑M2/212‑251 (MQ=255)
                            tttgtaaaccgggtcgggcaacggtttatacccatccgcagcgtttgagcaattaattaaaataagtcctgcgaaatatattaaattgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAg                                                                                                                                                                                                            <  1:64112‑M2/48‑1 (MQ=255)
                                                               ccgcagcgtttgagcaattaattaaaataagtcctgcgaaatatattaaattgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGt                                                                                                                                                                        >  2:220559‑M2/168‑251 (MQ=255)
                                                                               attaattaaaataagtcctgcgaaatatattaaattgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGAt                                                                                                                                                        >  2:154933‑M2/152‑251 (MQ=255)
                                                                                                                attgcaggacttattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGAttt                                                                                                                       >  1:42135‑M2/119‑251 (MQ=255)
                                                                                                                           tattcatttcgtgaattttattattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACAtct                                                                                                            >  1:223776‑M2/108‑251 (MQ=255)
                                                                                                                                                tattttatttataaacaaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAAc                                                                                        <  2:73361‑M2/164‑1 (MQ=255)
                                                                                                                                                                aaaatgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGCCGATTTACTGGCCATCCCGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAACCCATACCCTTACGaaa                                                                         >  2:414448‑M2/71‑249 (MQ=255)
                                                                                                                                                                   atgcgtgtagtggcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTAc                                                                    <  2:181377‑M2/184‑1 (MQ=255)
                                                                                                                                                                               gcgttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAATAGTACGGCATTGATAAt                                                        <  1:80530‑M2/196‑1 (MQ=255)
                                                                                                                                                                                  ttttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCAt                                                     >  1:139199‑M2/53‑251 (MQ=255)
                                                                                                                                                                                    ttcttttttcatttattatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCAttt                                                   >  1:67718‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                  attatgttcactctcttatctatatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATtta                                      >  1:387357‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                                       atatatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTg                <  1:281902‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                          tatttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCtt             >  1:185062‑M2/13‑251 (MQ=255)
                                                                                                                                                                                                                           atttataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTa            <  2:106018‑M2/240‑1 (MQ=255)
                                                                                                                                                                                                                              tataaggcTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTAcgc         >  2:414512‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                      ttcttATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTa  >  1:63924/1‑250 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
AAACAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTA  >  NC_000913/4100297‑4100776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.