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% 99.4 / NA 32intergenic (+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/20);  total (12/20)

GGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTG  >  NC_000913/1211521‑1211773
                                                                                                                                   |                                                                                                                         
ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGa                                                                                                                    >  1:127351/1‑139 (MQ=255)
  cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                                  <  1:264341/139‑1 (MQ=255)
  cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                                  <  1:289399/139‑1 (MQ=255)
   gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCa                                                                                                                 <  2:372088/139‑1 (MQ=255)
             tgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATg                                                                                                       <  2:286141/139‑1 (MQ=255)
                aaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAg                                                                                                    <  1:145539/139‑1 (MQ=255)
                    tttaaagTTTGGTCAAGCGATCATCGAAAACATGTAGTTTGTCCATATGTTAAATATTGAAAAGGCCGAATAACCAGCCTGTTGTTTTATTTATGTGAATATTATTAATAGTATACCGAGCATATTTATAGGAAGCCCa                                                                                                <  1:223272/134‑1 (MQ=255)
                    ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCa                                                                                                <  2:742119/139‑1 (MQ=255)
                    ttCAGAGTTAGGTGAAGCGATCATCAAATACATGTAATCTCTCCATGTGTTAAACAGTGCAACGGCCGTATACCACGCCCGGTGTTTTATGTATGTGTATATTATTAATAGTATAACGAGCATATTTATATGAAGCCCa                                                                                                <  1:458056/139‑1 (MQ=255)
                          gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACtt                                                                                          <  1:412834/139‑1 (MQ=255)
                             tGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAg                                                                                       <  2:127351/139‑1 (MQ=255)
                                         catcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa                                                                           <  2:610945/139‑1 (MQ=255)
                                                      tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                              <  2:479105/139‑1 (MQ=255)
                                                      tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                              >  2:311586/1‑139 (MQ=255)
                                                                     ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTaaa                                               >  1:252612/1‑139 (MQ=255)
                                                                       aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa                                                                >  1:216860/1‑120 (MQ=255)
                                                                       aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa                                                                <  2:216860/120‑1 (MQ=255)
                                                                       aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                         >  2:656300/1‑127 (MQ=255)
                                                                       aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                         <  1:656300/127‑1 (MQ=255)
                                                                             tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTTGATATTATTAATAGtata                                                                                                                        >  1:335299/1‑58 (MQ=255)
                                                                             tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGtata                                                                                                                        <  2:335299/58‑1 (MQ=255)
                                                                                         ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt                           <  1:556439/139‑1 (MQ=255)
                                                                                               cGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTcc                     >  2:801566/1‑139 (MQ=255)
                                                                                                  ccGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                         <  2:308151/100‑1 (MQ=255)
                                                                                                  ccGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACt                                                         >  1:308151/1‑100 (MQ=255)
                                                                                                     ttgttTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACta                                                        <  2:229909/98‑1 (MQ=255)
                                                                                                     ttgttTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACta                                                        >  1:229909/1‑98 (MQ=255)
                                                                                                                  gtgGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt                           <  2:102795/114‑1 (MQ=255)
                                                                                                                  gtgGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt                           >  1:102795/1‑114 (MQ=255)
                                                                                                                  gtgGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagg               >  2:372625/1‑126 (MQ=255)
                                                                                                                  gtgGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagg               <  1:372625/126‑1 (MQ=255)
                                                                                                                  gtgGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTg  >  1:563970/1‑139 (MQ=255)
                                                                                                                               aTAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                                       >  2:84457/1‑89 (MQ=255)
                                                                                                                               aTAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                                       <  1:84457/89‑1 (MQ=255)
                                                                                                                                   |                                                                                                                         
GGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTG  >  NC_000913/1211521‑1211773

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