Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,631 T→G V103V (GTT→GTG deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,6310TG96.3% 75.7 / ‑2.5 27coding (309/780 nt)deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (12/14);  minor base . (0/1);  total (12/15)
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.46e-01

ACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCAT  >  NC_000913/4617395‑4617876
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
cccGACGATAAAGTAATTGCGCTGTGACATCAGCCCAAACCCCCGGTCGCAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACGTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCTTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcg                                                                                                                                                                                                                                           <  1:152271/250‑1 (MQ=255)
  cgacgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcgcg                                                                                                                                                                                                                                         >  1:271164/1‑251 (MQ=255)
     cgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                                                                                                                                                                                                                                      >  1:344025/1‑251 (MQ=255)
                            caccagtcAAAACCCCGGGTCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt                                                                                                                                                                                                                <  1:262823/244‑1 (MQ=255)
                                                 ggCAATACCGCCGCTATCTGTATCTATCCTAGATTTCTCCCGTTTGCTCGCAAAACACTGCACGAGCAGGCCCCCCAGGAAATCCGTATTCCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGAGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGCCGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGTTTACGCGCAggttcgcttttacctctgcgactc                                                                                                                                                                                           >  2:182439/1‑231 (MQ=255)
                                                                atcggTATCTAACCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAACCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAAGCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGAGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTg                                                                                                                                                                           <  2:262054/247‑1 (MQ=255)
                                                                      atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTGAAGCGCAGGTTGGGTTTACACTGGTGAAAGCCttgtaagaggcgaggcgagg                                                                                                                                                                     >  2:64970/1‑232 (MQ=255)
                                                                               cGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCggaagt                                                                                                                                                            >  1:382374/1‑246 (MQ=255)
                                                                                   ttATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAATGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTAct                                                                                                                                                         <  1:338816/250‑1 (MQ=255)
                                                                                           aatgctcgataaccactGAAAGAGCAGGGCACCCCAGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGGAACGATGACATCGAAAGCGCGCTGGCAGAAACCCGTGCGGCAAGCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCATCGAATGTACTGCTGAAAga                                                                                                                                               <  2:120141/236‑2 (MQ=255)
                                                                                              gCTCGCAAAACACTGAAAGAGCCGGGCACCCCGGAAATCCGTATTGCGACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGa                                                                                                                                             <  2:245686/250‑1 (MQ=255)
                                                                                                     aaaCACTGAAAGAGCAGGGCACCCCGGAAATCCGTAGTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaa                                                                                                                                     <  1:182439/251‑1 (MQ=255)
                                                                                                      aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaaa                                                                                                                                    <  1:231837/247‑1 (MQ=255)
                                                                                                                             cGGAAATCCGTTTTTCTACGGTAACCAACCTCCCACACGGTTACGAGGACCTCGAAATCGCTCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAgcg                                                                                                             <  1:263416/251‑1 (MQ=255)
                                                                                                                                 aaTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGCTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGCATGTACTGATGGAAGGGAACATCGAAACCGGCGCACTggaacacgcacggctg                                                                                                          >  1:367953/1‑235 (MQ=255)
                                                                                                                                  aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGAt                                                                                                        <  1:211004/251‑1 (MQ=255)
                                                                                                                                    ccgtactgcgtcGGTAACCAACTTCCCACACGGTAATGATGACATCGAAATCGCGCTGGCAGAAACACTTGCGGCAATCTCCTGCTGTTCCGATGAAGTTGACGTGGTGTTGCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTGTTGGCCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGGTGAAAGTGATCATCGAAACCGGCGAACTGAAAGGCGAATCGCTGATcc                                                                                                      <  2:293604/240‑1 (MQ=255)
                                                                                                                                                         cTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAAACGGCGCACTGAAAGACGAAACGCTGATCCGTAAAGCGACTGAAATCTcc                                                                                  >  1:53209/1‑250 (MQ=255)
                                                                                                                                                             ccACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAACTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGGCGCGAATGTACTGCTGAAAGAGGTCATCGGAAACGGCGAAATGAAAAACGAAGCGCTGATCCGGAAAGCGGCTGAAATCTCCCTCa                                                                              >  1:239719/1‑250 (MQ=255)
                                                                                                                                                                        cgatgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGGTTTGAACTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGGACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTATGAAATCTCCATCAAAGCGGGGGCgg                                                                  >  1:72879/5‑251 (MQ=255)
                                                                                                                                                                               atcgaaatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGCAAAGGGTCATCGAAACCGGCGCACTGAAAAACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCACTCCACGCGGGTGAGGACTtcat                                                           >  1:257366/1‑251 (MQ=255)
                                                                                                                                                                                                   aaCCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGGCCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGAGAATGTACTGCTGAAAGTGATCATCGGAACCGGAGAACTGAAAGACGAAGCGCTGGTCCGTAAAACGTCTGAAATCTCCCTCAAAACGGGTGCGGACTTCATCAAAACCCCTACCCGTaaa                                        >  1:96097/1‑250 (MQ=255)
                                                                                                                                                                                                         ggCGGCAAGCGCCTGCGGGTCCGATGAAGTTGACGTGTTGTTCCCGGACCGCGCGCAGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGCATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCTTAAAGCGTCTGAAATCTCCATCAAAGCGTGTGCGGACTTCATCAAAACCTCTTCCGGTAAAGTGGCtg                                 <  2:175052/250‑1 (MQ=255)
                                                                                                                                                                                                               aaTCGCCTACGGTGCCGATGAAGTTGACGTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAgc                                                                                                              <  1:48977/167‑1 (MQ=255)
                                                                                                                                                                                                                      tACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGGAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCCTCCAAGCGGGGGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCgccgcc                    >  1:354218/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGcatcat  <  1:227430/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gACGTGGTGGTCCCGTACCGAGCGGTGATGGCGGGTAACGAGCAGTGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACTCCGGAAAGCGCGCGcatcat  <  2:239719/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
ACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCAT  >  NC_000913/4617395‑4617876

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