Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,565 C→T D81D (GAC→GAT deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,5650CT95.7% 63.3 / NA 23D81D (GAC→GATdeoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/12);  minor base G (0/1);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCA  >  NC_000913/4617332‑4617800
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
cTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCgcgc                                                                                                                                                                                                                              >  1:366178/1‑250 (MQ=255)
                  cGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                                                                                                                                           >  2:112977/1‑251 (MQ=255)
                      cACTGAAATTGAGGGCCCTGCCCCCCCTGAAGGACGACGACACCGACGAGAATGTAATTGCTTTGTTTCTCCGGGCCAAAACCCCGGTCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                                                                       <  1:221741/251‑1 (MQ=255)
                                     ccatgcccccccTGAATGACGACGACCCCGACGAAAAAGTAATTGCTCTGTGTCATCGGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgat                                                                                                                                                                                        <  1:112977/242‑1 (MQ=255)
                                                               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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
CTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCA  >  NC_000913/4617332‑4617800

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