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

GAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGC  >  NC_000913/4492352‑4492827
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
gAGCGTGCTGCAGATAACATCTTGTATAGAGTGCAGGACGGCCTGTGCATACATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATAATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTa                                                                                                                                                                                                                                     >  2:28743/1‑250 (MQ=255)
           cGATAACATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAt                                                                                                                                                                                                                         >  2:247125/1‑251 (MQ=255)
                                                                                        aCGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGc                                                                                                                                            >  2:88983/1‑251 (MQ=255)
                                                                                          gaggaaaaTAAAACCGTAAATTTGAAAGACGATCATAAGGGTAAAAAGCAGCGGGACGCCTGCCGCAAAAATGATTAATGGCATAGTGTTTATCTCTTATTAGTGCAAGAAAAGCAGCGCTTACGCGCTGCACGG‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAAGCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGGTCCTGCgg                                                                                                                                          <  2:70635/245‑1 (MQ=255)
                                                                                              gCAATAAAGCCGTTAACTTTAAAGACGAGCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAtct                                                                                                                                      <  1:168026/251‑1 (MQ=255)
                                                                                                       gcGTTAACTTTACAGAAGATTATCAGGATAAGAAGCAGCCCGACGCCTGCCGCAATAATGATTAATGGCATAGCCTTTATCTCTTATTAGTACCAGAAAAGCACCGCTTTCGCGCTGCACGTCGGGGGGG‑CTGTTAAACAGCCCATAACCTGCCGCCATCAACAACCAACAGGTGGCAGTTTAAGCAATCAGAGGCTTTTAACGAAAGGAACACAGCAGCAACAATCATTTcaagcggtgttccccagcg                                                                                                                              >  2:95364/2‑232 (MQ=255)
                                                                                                        cGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTg                                                                                                                            >  2:100353/1‑251 (MQ=255)
                                                                                                           tAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCAATAAACAAACACCAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAACAATCAGTTCCTGCGGATCTACCCAGCGTgcgg                                                                                                                         >  1:120297/1‑251 (MQ=255)
                                                                                                                            atcaGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        >  1:8125/1‑251 (MQ=255)
                                                                                                                               aGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAc                                                                                                     >  1:195162/1‑251 (MQ=255)
                                                                                                                                      gaagCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCgg                                                                                              >  1:187522/1‑251 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        >  1:94697/1‑237 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcac                                                                                                        <  2:94697/237‑1 (MQ=255)
                                                                                                                                          cagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGaa                                                                                          <  1:88983/251‑1 (MQ=255)
                                                                                                                                                              aaaggATTAAGGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTg                                                                      <  2:195162/247‑1 (MQ=255)
                                                                                                                                                                                                                  tttCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATCCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGAAAGGCTTCGTCCTCAACCAGTGCTTTAGTCCTTTCTGTTTTGAAATAGCCCGGCGCCATAACGTTTGCCTGAat                  >  1:230379/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcacttaggggggg‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAACTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGgc          <  2:200902/242‑1 (MQ=255)
                                                                                                                                                                                                                                  ggggggg‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGAGGGGGTCCGTTTGCACAACCAGGCGGAGAAGGCTTCGTCCTCAACCCGGGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAAAACCGTTGCCCTAAATATTGTGGCGCGcccgc  >  1:262619/2‑248 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
GAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGC  >  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: ‑