Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,211,652 C→T intergenic (+75/+28) icdC → / ← iraM pseudogene, isocitrate dehydrogenase C‑terminal gene fragment, idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,211,6520CT100.0% 84.7 / NA 28intergenic (+75/+28)icdC/iraMpseudogene, isocitrate dehydrogenase C‑terminal gene fragment, idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/16);  total (12/16)

AGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCT  >  NC_000913/1211520‑1211775
                                                                                                                                    |                                                                                                                           
aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGTCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCg                                                                                                                       <  1:708361/139‑1 (MQ=255)
  gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     <  2:1257231/139‑1 (MQ=255)
  gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     <  1:27642/139‑1 (MQ=255)
  gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     >  1:418479/1‑139 (MQ=255)
       aaGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCatat                                                                                                                <  1:1828379/139‑1 (MQ=255)
                           gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACtt                                                                                            <  2:418479/139‑1 (MQ=255)
                               ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCata                                                                                                                 >  1:748249/1‑114 (MQ=255)
                               ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCata                                                                                                                 <  2:748249/114‑1 (MQ=255)
                                        atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                                                                               <  1:592048/139‑1 (MQ=255)
                                             cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCatat                                                                          <  2:1605979/139‑1 (MQ=255)
                                                       tAATCTCTCCATGTGTTAAATCTTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     >  1:1347146/1‑86 (MQ=255)
                                                       tAATCTCTCCATGTGTTAAATATTGAACCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     <  2:1347146/86‑1 (MQ=255)
                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     >  2:373421/1‑76 (MQ=255)
                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                     <  1:373421/76‑1 (MQ=255)
                                                                              tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa                                                                             >  1:1491685/1‑103 (MQ=255)
                                                                              tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa                                                                             <  2:1491685/103‑1 (MQ=255)
                                                                                aaaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                                       >  2:373401/1‑139 (MQ=255)
                                                                                                  cccGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTg                     >  1:704855/1‑139 (MQ=255)
                                                                                                      ttgttTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                                                                               >  2:970928/1‑77 (MQ=255)
                                                                                                      ttgttTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                                                                               <  1:970928/77‑1 (MQ=255)
                                                                                                             atttatGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAg                                                                                         <  2:1823216/60‑1 (MQ=255)
                                                                                                             atttatGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAg                                                                                         >  1:1823216/1‑60 (MQ=255)
                                                                                                                        tattatTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCt  >  2:1702748/1‑136 (MQ=255)
                                                                                                                        tattatTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCt  <  1:1702748/136‑1 (MQ=255)
                                                                                                                            attaATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                           >  2:1557969/1‑75 (MQ=255)
                                                                                                                            attaATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                           <  1:1557969/75‑1 (MQ=255)
                                                                                                                                aTAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                                >  1:415896/1‑66 (MQ=255)
                                                                                                                                aTAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                                <  2:415896/66‑1 (MQ=255)
                                                                                                                                    |                                                                                                                           
AGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCT  >  NC_000913/1211520‑1211775

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