Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,310,841 C→T 40.9% E110K (GAA→AAA)  pnp ← polynucleotide phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,310,8410CT40.9% 3.5 / 14.0 20E110K (GAA→AAA) pnppolynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (2/6);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATA  >  NC_000913/3310706‑3310982
                                                                                                                                       |                                                                                                                                               
cTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTAGCGGGTTAACAGAAACCACGGTGGCGATAACCTGCACTTCGTTGACg                                                                                                                                        >  2:153422/1‑143 (MQ=255)
    gcgGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGAGAAGATAGAAGTGTGACACAACTGCCCAACAGACACCACGGTGGCGAGACCCGGCACTTCGTGGCCGAAGc                                                                                                                                    >  2:67443/1‑143 (MQ=255)
       ggCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAACACTGCGGAAGCACCAATCATCGCGACGATATCCGGGATAACTAGCGGGATAACAGAAACAACGATGGCGATAACCCGAACCTTGTTGACGAAGCCtt                                                                                                                                 >  2:200730/1‑143 (MQ=255)
        gcagcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTc                                                                                                                                <  1:236231/143‑1 (MQ=255)
             aCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGa                                                                                                                           >  1:207379/1‑143 (MQ=255)
                     ggCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCggg                                                                                                                   <  2:204606/143‑1 (MQ=255)
                     ggCCAATCAACCGAATACCCGACAGCGACAGCGCTGCCGAAGCACCAATCCTCGCGACGATATCCGGGTTACCTTACCTGGTTATAGAATCCACGGTGGCGATAACCTGACCTTCGTTGACGAAGCCTTCCGGGAACAGCggg                                                                                                                   <  2:127991/143‑1 (MQ=255)
                          ttGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAAt                                                                                                              >  2:71916/1‑143 (MQ=255)
                                       cagacagaGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATc                                                                                                 <  1:57445/143‑1 (MQ=255)
                                               gaccgggcTGCGGAAGCACCAGATATCGCGACGTTATCCGGGGTAACTGGCGGGTTGACAGAAACGACGGTGGCGATAAGCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGAcgcgc                                                                                         <  2:223637/137‑1 (MQ=255)
                                                           gAAGCACCAATCAGCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTc                                                                             <  1:59276/143‑1 (MQ=255)
                                                                cccAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCtt                                                                        <  1:179427/142‑1 (MQ=255)
                                                                 ccAATCAGCGCGACGAGATCTGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTc                                                                       <  2:7485/143‑1 (MQ=255)
                                                                  cAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGCCAAAACCGACGTTATGACGTGCAATTAGTGTGGTTTTTGGCT‑‑CGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGGTTCGCCTTCg                                                                      >  1:235401/1‑143 (MQ=255)
                                                                   aaTCATCCCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTTGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGc                                                                     <  2:166683/143‑1 (MQ=255)
                                                                            gacgcTATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGggcggc                                                            <  1:156342/143‑1 (MQ=255)
                                                                               gaTATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGAGCAGGGTTTCGCCTTCGCTTGGGCGGCCtt                                                         >  2:24478/1‑143 (MQ=255)
                                                                                         ttAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGgaa                                               >  1:107715/1‑143 (MQ=255)
                                                                                                                                   aCTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTGCGCGTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTg     <  1:24478/143‑1 (MQ=255)
                                                                                                                                      tCGCTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCCATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATa  <  1:134362/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                               
CTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATA  >  NC_000913/3310706‑3310982

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