Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,804,405 T→C S302P (TCC→CCC)  nrdF → ribonucleoside‑diphosphate reductase 2 subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,804,4050TC94.4% 55.9 / NA 18S302P (TCC→CCC) nrdFribonucleoside‑diphosphate reductase 2 subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/9);  minor base A (1/0);  total (9/9)
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

AATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAA  >  NC_000913/2804276‑2804534
                                                                                                                                 |                                                                                                                                 
aaTAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGcccct                                                                                                                                >  2:455819/1‑133 (MQ=255)
aaTAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGcccct                                                                                                                                >  2:402760/1‑133 (MQ=255)
      gCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCCCGAATGCATTTCCGGTTCAGGCACCTCTTAtg                                                                                                                          >  2:232892/1‑133 (MQ=255)
        tttGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTAtgtg                                                                                                                        >  1:32247/1‑133 (MQ=255)
         ttGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGa                                                                                                                       >  2:376593/1‑133 (MQ=255)
          tggtgaATCTGGGCTACAAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGAt                                                                                                                      >  2:448112/4‑133 (MQ=255)
                   tGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGc                                                                                                             <  1:455819/133‑1 (MQ=255)
                      gCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGt                                                                                                          >  2:400942/1‑133 (MQ=255)
                            aaCCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAAc                                                                                                    <  1:402760/133‑1 (MQ=255)
                             aCCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACa                                                                                                   <  1:376593/133‑1 (MQ=255)
                                   aTTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGAt                                                                                             <  2:32247/133‑1 (MQ=255)
                                       cccGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAg                                                                                         >  1:393759/1‑133 (MQ=255)
                                           cAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTg                                                                                     <  1:400942/133‑1 (MQ=255)
                                                        aGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAggg                                                                        >  1:387500/1‑133 (MQ=255)
                                                                            ccgcGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCa                                                    >  1:469870/1‑133 (MQ=255)
                                                                                      gCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCgcag                                          <  2:387500/133‑1 (MQ=255)
                                                                                       ccGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCgcagc                                         >  2:90130/1‑133 (MQ=255)
                                                                                         gAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGcagcag                                       <  1:448112/133‑1 (MQ=255)
                                                                                                                             aGGCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCa   <  2:477678/133‑1 (MQ=255)
                                                                                                                              ggCCCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCaa  <  1:90130/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                 
AATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGATGAAAATCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAA  >  NC_000913/2804276‑2804534

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