Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 417,139 2043 bp→CATAGGTGGCGGCCGCTCTAGAACTAGTGGATCCCCC phoB[phoR] phoB, [phoR]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 417139 419181 2043 18 [0] [0] 17 phoB–[phoR] phoB,[phoR]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4171380 (0.000)13 (0.380)
+CATAGGTGGCGGCCGCTCTAGAACTAGTGGATCCCCC
13/406 1.1 100% intergenic (‑186/‑4) sbcD/phoB ATP‑dependent structure‑specific DNA nuclease ‑ SbcD subunit/DNA‑binding transcriptional dual regulator PhoB
?NC_000913 419182 = 0 (0.000)coding (1294/1296 nt) phoR sensor histidine kinase PhoR

CGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTAT  >  NC_000913/418947‑419400
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                          
aaggatgcgcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCa                                                                                                                                                                                                             <  1:411624‑M2/16‑1 (MQ=255)
                                            ttatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTg                                                                                                                                                                 >  2:118826‑M2/192‑251 (MQ=255)
                                                                                            caataaccggaagatatgtgcgacgagcttttcataaagctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctattgattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGACATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGAttgtttgt                                                                                                                 <  2:339360‑M2/108‑1 (MQ=255)
                                                                                                                cgacgagcttttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCAGAAGCTAGTCGATGCTGGCTTATTTTCTGTGCAGTCAAAATACGGGCGTTAGATTTTACAACTATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCt                                                                                             <  1:187243‑M2/128‑1 (MQ=255)
                                                                                                                       cttttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTg                                                                                                                                                                                                                >  2:299083‑M2/117‑129 (MQ=255)
                                                                                                                       cttttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTg                                                                                                                                                                                                                <  1:299083‑M2/13‑1 (MQ=255)
                                                                                                                          ttcataaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACgg                                                                                    >  1:129377‑M2/114‑250 (MQ=255)
                                                                                                                               aaatctgtcataaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTGGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAg                                                                              <  1:241372‑M2/143‑1 (MQ=255)
                                                                                                                                          aaatctgacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACt                                                                   >  2:5473‑M2/98‑251 (MQ=255)
                                                                                                                                                                 cattaatgatcgcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCtat                                            >  1:274827‑M2/75‑251 (MQ=255)
                                                                                                                                                                            gcaacctatttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGt                                 >  1:303582‑M2/64‑251 (MQ=255)
                                                                                                                                                                                     ttattacaacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTATTAGAt                         <  1:118826‑M2/196‑1 (MQ=255)
                                                                                                                                                                                            aacagggcaacataggtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCAGGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCTCTATATACCTCTGGTTTTTATATCCCTCCtt                 <  2:255353‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                      cataggtggcggctgctctagaactagtggatcccccTAATCCACCTTTGTCATCTTTTATTGCCATAAGCCTGTCGATGCTGGCTAATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGAGTTTTTTACTGGATGGTGATTTCTTATCGATATATACCTCTGGTTTTTAGATCCCTCCTGGCTTTAAAAc       >  2:106182‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                           gtggcggccgctctagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTtat  <  1:155341‑M2/219‑1 (MQ=255)
                                                                                                                                                                                                                        tagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAg                                                                              >  2:140644‑M2/20‑162 (MQ=255)
                                                                                                                                                                                                                        tagaactagtggatcccccTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAg                                                                              <  1:140644‑M2/143‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                          
CGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTAT  >  NC_000913/418947‑419400

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