Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,235,905 16376 bp→AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC yjbEmalM yjbE, yjbF, yjbG, yjbH, yjbT, psiE, xylE, malH, malG, malF, malE, malK, lamB, malM

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4235905 4252280 16376 20 [0] [0] 19 yjbE–malM yjbE,yjbF,yjbG,yjbH,yjbT,psiE,xylE,malH,malG,malF,malE,malK,lamB,malM

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42359040 (0.000)13 (0.540)
+AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC
12/382 0.8 100% intergenic (+497/‑2) pgi/yjbE glucose‑6‑phosphate isomerase/uncharacterized protein YjbE
?NC_000913 4252281 = 0 (0.000)intergenic (+1807/‑225) malM/ubiC maltose regulon periplasmic protein/chorismate lyase

GGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACC  >  NC_000913/4252070‑4252515
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
aaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgaacccgcccgaaatgtttttgggtaaatctccatgcattcaatgaagggaaattgtaagtaaaggaaaacatcagcgttcgcaaaagaaaatgccgatatgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGAc                                                                                                                                                                                                      <  1:177128‑M2/39‑1 (MQ=255)
            aaacttatatttgatattcattccaattatcctaaaacgccatcgctaattcccagcgccgtaattcgcatgctttagctgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgcagatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGCCATCATGGACGGTGAGGCCCCt                                                                                                                                                                                          >  2:56551‑M2/200‑250 (MQ=255)
                            tgcattccaattatcctctatcgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgacatgattttgggtaaatctccattcattcaatgacgggaaattgtaagtacaggaaaacatccacgttcgcaaaggaaaatgacgatttgcggccgcTCAGTAAAGGAAAACCTCACCGTTCCTGGCATCCTGGACGGTGATGGGCCTACGGATACCCTCGcc                                                                                                                                                                           >  2:139253‑M2/184‑249 (MQ=255)
                                       tatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATc                                                                                                                                                               <  2:92597‑M2/78‑1 (MQ=255)
                                                                      ccgtaattcgcatgctttagttgtgtatactcgatcccgaccgaactgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTaccc                                                                                                                               >  1:177005‑M2/142‑250 (MQ=255)
                                                                        gtaattcgcatgctttagttgtgtatactcgatcccgcacgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATAGGTAAAGCGTAAGGTTCAACATCGTTTTACCACt                                                                                                                             >  1:278823‑M2/140‑251 (MQ=255)
                                                                                                           cgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTAtt                                                                                              >  1:162704‑M2/105‑247 (MQ=255)
                                                                                                                                       attcattcactgaagggaaattgtaagtaaaggaaaacatcaaagttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTa                                                               >  1:50348‑M2/77‑250 (MQ=255)
                                                                                                                                         tcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgacttccggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAac                                                            >  1:53585‑M2/75‑251 (MQ=255)
                                                                                                                                            ttcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACt                                                         >  1:22280‑M2/72‑251 (MQ=255)
                                                                                                                                                       gaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtct                                                                     >  2:181856‑M2/61‑228 (MQ=255)
                                                                                                                                                       gaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtct                                                                     <  1:181856‑M2/168‑1 (MQ=255)
                                                                                                                                                          attgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACaa                                           <  2:76860‑M2/194‑1 (MQ=255)
                                                                                                                                                                                           aagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCAt          >  1:133372‑M2/25‑251 (MQ=255)
                                                                                                                                                                                                   gccgatttgcggacgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGCTGCACCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCCCTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGACtgtcacccc  >  2:107560‑M2/17‑248 (MQ=255)
                                                                                                                                                                                                       atttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCGCCCGTAAACACGTGTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTc               <  2:289634‑M2/222‑1 (MQ=255)
                                                                                                                                                                                                       atttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGATGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTCAAGCGTCAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGCTACAATGCCTTTACGTTATGTAACGGAGAGTTc               >  1:289634‑M2/13‑234 (MQ=255)
                                                                                                                                                                                                             ggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCGCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATGGTTGCGTTTttgt                                                                                                      <  2:123141‑M2/135‑1 (MQ=255)
                                                                                                                                                                                                             ggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCACCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGCAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTttgt                                                                                                      >  1:123141‑M2/7‑141 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACC  >  NC_000913/4252070‑4252515

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