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,6520CT94.4% 48.3 / ‑6.0 18intergenic (+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 (1/0);  new base T (5/12);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTG  >  NC_000913/1211517‑1211756
                                                                                                                                       |                                                                                                        
ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGtata                                                                                                       <  2:128762/139‑1 (MQ=255)
     gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                  <  2:855394/139‑1 (MQ=255)
       gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCa                                                                                                <  1:507383/139‑1 (MQ=255)
                                ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGa                                                                       <  1:820483/139‑1 (MQ=255)
                                       aGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTATTTGAGCCCata                                                                <  1:31588/139‑1 (MQ=255)
                                             catcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa                                                          <  2:1159300/139‑1 (MQ=255)
                                                     aCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta                                                  <  2:722873/139‑1 (MQ=255)
                                                      cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat                                                 <  1:891691/139‑1 (MQ=255)
                                                        tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa                                               <  1:826907/139‑1 (MQ=255)
                                                         gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt                                              >  1:1167309/1‑139 (MQ=255)
                                                           aaTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGc                                            <  2:99629/139‑1 (MQ=255)
                                                                 tcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGaa                                                                                    <  2:60986/93‑1 (MQ=255)
                                                                 tcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGaa                                                                                    >  1:60986/1‑93 (MQ=255)
                                                                 tcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtat                                      >  2:960489/1‑139 (MQ=255)
                                                                           aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                 <  2:281364/70‑1 (MQ=255)
                                                                           aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                 >  1:281364/1‑70 (MQ=255)
                                                                                        ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGATTCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTTTTATGTAAAAATTTAGTTTTTAtt               >  2:804978/1‑139 (MQ=255)
                                                                                                   cGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTcc    >  2:461710/1‑139 (MQ=255)
                                                                                                     cccGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTg  <  2:486739/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                        
GGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTG  >  NC_000913/1211517‑1211756

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