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% 80.6 / NA 26intergenic (+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 (10/16);  total (10/16)

AGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTTT  >  NC_000913/1211520‑1211777
                                                                                                                                    |                                                                                                                             
aGGCGCTAAGCTGCTGCAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCg                                                                                                                         <  1:791805/139‑1 (MQ=255)
  gggcTAAGCTGCTGAAATGTTCAGAGGTTGGTGAAGCGATCATCGAAAACATGTAATATCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCTCGTTGTTTTATTTATGTTGATATTATTAATAGTATATCGAg                                                                                                                       <  2:658710/137‑1 (MQ=38)
   cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                       <  1:28592/138‑1 (MQ=255)
   cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                       >  2:28592/1‑138 (MQ=255)
   cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                                      <  2:640258/139‑1 (MQ=255)
     cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                       <  1:1055931/136‑1 (MQ=255)
     cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                       >  2:1055931/1‑136 (MQ=255)
     cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCat                                                                                                                    >  1:46261/1‑139 (MQ=255)
               gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGa                                                                                                          <  2:253631/139‑1 (MQ=255)
                                      cGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                                                                                   >  2:335971/1‑139 (MQ=255)
                                        atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAg                                                                                                        >  1:848078/1‑116 (MQ=255)
                                        atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAg                                                                                                        <  2:848078/116‑1 (MQ=255)
                                                     tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATTTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa                                                                    <  2:794925/139‑1 (MQ=255)
                                                       tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                                  <  1:335971/139‑1 (MQ=255)
                                                             ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCTCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGTCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACta                                                            <  1:472997/139‑1 (MQ=255)
                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCa                                                                                                    >  2:716503/1‑95 (MQ=255)
                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCa                                                                                                    <  1:716503/95‑1 (MQ=255)
                                                                              tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                                           >  1:896442/1‑139 (MQ=255)
                                                                                    gggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCa                                                                                                                     >  1:1062205/1‑59 (MQ=255)
                                                                                    gggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCa                                                                                                                     <  2:1062205/59‑1 (MQ=255)
                                                                                       cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATg                                                                                                           <  2:291034/66‑1 (MQ=255)
                                                                                       cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATg                                                                                                           >  1:291034/1‑66 (MQ=255)
                                                                                                                       atatTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCttt  <  2:527703/139‑1 (MQ=255)
                                                                                                                            attaATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat                                                                      >  1:959483/1‑66 (MQ=255)
                                                                                                                            attaATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat                                                                      <  2:959483/66‑1 (MQ=255)
                                                                                                                                 tAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa                                                                    <  2:1049812/63‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
AGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTTT  >  NC_000913/1211520‑1211777

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