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 21 [0] [0] 35 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)21 (0.580)
+AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC
21/384 0.6 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

TTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAA  >  NC_000913/4235708‑4236152
                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
ttGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATCTTCATTCCCATTATCCTAACACGCCATAGCTAATTCCCCGCGCCGTACTTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgct                                                                                                                                                                                                     >  1:307759‑M1/1‑197 (MQ=255)
           ttaAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaa                                                                                                                                                                                          >  1:91652‑M1/1‑186 (MQ=255)
               gTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacat                                                                                                                                                                                     >  1:326247‑M1/1‑182 (MQ=255)
                            aTTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggcagctaagtaaaggaaaacatcaccgttcctggc                                                                                                                                                                        >  2:36291‑M1/1‑169 (MQ=255)
                               gTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatc                                                                                                                                                                     >  1:350403‑M1/1‑166 (MQ=255)
                                                tGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccc                                                                                                                                                    >  1:364705‑M1/1‑149 (MQ=255)
                                                    aCTTATATTTGATATTCATTCCAATAATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacg                                                                                                                                                >  1:243400‑M1/1‑145 (MQ=255)
                                                                                 cTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacg                                                                                                                                                <  1:61667‑M1/222‑107 (MQ=255)
                                                                                 cTAAAACGCCATCGCTAATTACCCGCGCCGTACTTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcatggacggtgatgcccctacg                                                                                                                                                >  2:61667‑M1/1‑116 (MQ=255)
                                                                                   aaaaCGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaa                                                                                                                   <  1:21508‑M1/249‑136 (MQ=255)
                                                                                        gCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaa                                                                                                            <  1:206450‑M1/251‑143 (MQ=255)
                                                                                        gCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgta                                                                                                             >  1:159259‑M1/1‑109 (MQ=255)
                                                                                                 aaTTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaac                                                                                                    >  1:234655‑M1/1‑100 (MQ=255)
                                                                                                                     gCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcg                                                                               >  1:116695‑M1/1‑80 (MQ=255)
                                                                                                                         gCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgagttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccttacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattg                                                                           <  2:326247‑M1/251‑176 (MQ=255)
                                                                                                                                            gATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGACGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtat                                                        >  1:321512‑M1/1‑57 (MQ=255)
                                                                                                                                                        aaTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcatggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcactta                                             >  2:406641‑M1/1‑45 (MQ=255)
                                                                                                                                                                     aaaTCTCCAGTCATTCAAGGAAAGGAAATGGTaaggaaaggaaaacagcaacgttcgcaaaagaagatgccgatttgcgggcgctaagtaaaggaacacatcaccgtccctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagagcacttaatttgctgtacat                                <  2:256357‑M1/250‑219 (MQ=255)
                                                                                                                                                                          tcCATTCATTAAATGAAGGGACATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtacaggaaaacatcaacgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccg                          >  2:319168‑M1/1‑27 (MQ=255)
                                                                                                                                                                                    caATGAAGGGCAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacacttt                >  1:290049‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                   gtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgctgccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacacttttctgcgatacaatg  <  2:109528‑M1/250‑249 (MQ=255)
                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
TTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAA  >  NC_000913/4235708‑4236152

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