Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,583 (G)6→7 intergenic (‑58/+4) idnT ← / ← idnO L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,5831.G93.8% 43.9 / ‑3.8 16intergenic (‑58/+4)idnT/idnOL‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (10/5);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGC  >  NC_000913/4492352‑4492827
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
gAGCGTGCTGCAGATAACATCTTGTATAGAGTGCAGGACGGCCTGTGCATACATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATAATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTa                                                                                                                                                                                                                                     >  2:28743/1‑250 (MQ=255)
           cGATAACATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAt                                                                                                                                                                                                                         >  2:247125/1‑251 (MQ=255)
                                                                                        aCGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGc                                                                                                                                            >  2:88983/1‑251 (MQ=255)
                                                                                          gaggaaaaTAAAACCGTAAATTTGAAAGACGATCATAAGGGTAAAAAGCAGCGGGACGCCTGCCGCAAAAATGATTAATGGCATAGTGTTTATCTCTTATTAGTGCAAGAAAAGCAGCGCTTACGCGCTGCACGGAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAAGCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGGTCCTGCgg                                                                                                                                          <  2:70635/245‑1 (MQ=255)
                                                                                              gCAATAAAGCCGTTAACTTTAAAGACGAGCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAtct                                                                                                                                      <  1:168026/251‑1 (MQ=255)
                                                                                                       gcGTTAACTTTACAGAAGATTATCAGGATAAGAAGCAGCCCGACGCCTGCCGCAATAATGATTAATGGCATAGCCTTTATCTCTTATTAGTACCAGAAAAGCACCGCTTTCGCGCTGCACGTCGGGGGGGCTGTTAAACAGCCCATAACCTGCCGCCATCAACAACCAACAGGTGGCAGTTTAAGCAATCAGAGGCTTTTAACGAAAGGAACACAGCAGCAACAATCATTTcaagcggtgttccccagcg                                                                                                                              >  2:95364/2‑232 (MQ=255)
                                                                                                        cGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTg                                                                                                                            >  2:100353/1‑251 (MQ=255)
                                                                                                           tAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCAATAAACAAACACCAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAACAATCAGTTCCTGCGGATCTACCCAGCGTgcgg                                                                                                                         >  1:120297/1‑251 (MQ=255)
                                                                                                                            atcaGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        >  1:8125/1‑251 (MQ=255)
                                                                                                                               aGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAc                                                                                                     >  1:195162/1‑251 (MQ=255)
                                                                                                                                      gaagCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCgg                                                                                              >  1:187522/1‑251 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        >  1:94697/1‑237 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        <  2:94697/237‑1 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGaa                                                                                          <  1:88983/251‑1 (MQ=255)
                                                                                                                                                              aaaggATTAAGGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTg                                                                      <  2:195162/247‑1 (MQ=255)
                                                                                                                                                                                                                  tttCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATCCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGAAAGGCTTCGTCCTCAACCAGTGCTTTAGTCCTTTCTGTTTTGAAATAGCCCGGCGCCATAACGTTTGCCTGAat                  >  1:230379/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcacttagggggggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAACTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGgc          <  2:200902/242‑1 (MQ=255)
                                                                                                                                                                                                                                  gggggggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGAGGGGGTCCGTTTGCACAACCAGGCGGAGAAGGCTTCGTCCTCAACCCGGGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAAAACCGTTGCCCTAAATATTGTGGCGCGcccgc  >  1:262619/2‑248 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
GAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGC  >  NC_000913/4492352‑4492827

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