Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,148,446 G→A 13.1% E184K (GAG→AAG)  gpr → L‑glyceraldehyde 3‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,148,4460GA13.1% 55.4 / 3.3 31E184K (GAG→AAG) gprL‑glyceraldehyde 3‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base A (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

TCTCATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTACTCCTCTGGCTCAGGGATTGCTGAC  >  NC_000913/3148305‑3148588
                                                                                                                                             |                                                                                                                                              
tctcATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTCCTCGCCAGGGCGGGCGCACAACCTGGTCGGGTTGCTGCGCGa                                                                                                                                               >  1:225515/1‑143 (MQ=255)
   cATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAACAAATGGTCGAGTTGCTCCGCGAGTg                                                                                                                                            >  2:264112/1‑143 (MQ=255)
             aTGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGCCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGCCGAAAATACCg                                                                                                                                  >  2:552849/1‑143 (MQ=255)
                aaaaTACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGCGCGGACGCAAACAATGGTCGAGTTGCTCCGCGAGCGGAAAATTCCGCTg                                                                                                                               >  2:318054/1‑143 (MQ=255)
                            tGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAGAATCCGTCGTTAATtcaccca                                                                                                                   >  2:491375/1‑139 (MQ=255)
                              gaagaaACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCAAGTGGAAAATTCCGCTGTTAATTCATCAAcc                                                                                                                 >  1:172181/1‑143 (MQ=255)
                              gaagaaACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCAAGTGGAAAATTCCGCTGTTAATTCATCAAcc                                                                                                                 >  1:314470/1‑143 (MQ=255)
                                 gaaACCGCCTCTGCGCCGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTc                                                                                                              <  1:286200/143‑1 (MQ=255)
                                                       aTGCGGTACCAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGGCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTgg                                                                                        <  1:117047/143‑1 (MQ=255)
                                                       aTGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGATAGAGTTGCTGCGCGAGTGGAAACTTCCGCTGTTAATTCATCAACCTTCGTACAATGTACTGAACCGCTgg                                                                                        >  1:253553/1‑143 (MQ=255)
                                                          cGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTg                                                                                     <  1:320722/143‑1 (MQ=255)
                                                                aaaGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGCACAATTTACTGAACCGCGGGGTGGATaaa                                                                               >  2:261295/1‑143 (MQ=255)
                                                                     ggTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCCGTTAATTCATCAAACTTCGTCCAATTTACTTAACCGCTGGGTGGATAAAAACgg                                                                          >  2:438069/1‑143 (MQ=255)
                                                                       tAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGcc                                                                        >  2:123496/1‑143 (MQ=255)
                                                                          ggCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCctgc                                                                     <  1:287284/143‑1 (MQ=255)
                                                                              cTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAGATGGTCGAGTTGCTGACCGAGTGGAAAATTCCGCTGTTAATTCATCACCCTTCGTACAAGTTACTGCACCGGTGGGTGGCTAAAAGCGGCCTGCTGGa                                                                 >  2:381737/1‑143 (MQ=255)
                                                                                 taggTCGGGATCTACTCGTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATAc                                                              <  2:225515/140‑1 (MQ=255)
                                                                                    gTCGGGATCTCCGCTAACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCTTCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCt                                                           <  1:80985/143‑1 (MQ=255)
                                                                                               ctctTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCAAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACg                                                <  2:314470/143‑1 (MQ=255)
                                                                                               ctctTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCAAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACg                                                <  2:172181/143‑1 (MQ=255)
                                                                                                   cgcccGCCAGAGCGGACGCAAACAATGGTCGAGTTGCTGCGCGAGTGGAACATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCGGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGt                                            <  1:138149/139‑1 (MQ=255)
                                                                                                        gccacctcGGACGCAAATAATGGCCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCt                                       <  1:220161/136‑1 (MQ=255)
                                                                                                         ccAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACTCTGCAAAATAACGGCGTGGGCTg                                      >  1:84449/1‑143 (MQ=255)
                                                                                                           agagCCGACGCAAAAAATCGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGCACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTa                                    <  1:464766/143‑1 (MQ=255)
                                                                                                             agCGGACGCCAAACCGGGTCGAGTTGCTGCGCGAGTGGAAAACTCCGCGGTTAATTCATCAACCTTCGTCCAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTAtt                                  <  1:295904/143‑1 (MQ=255)
                                                                                                              gCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTg                                 <  1:139446/143‑1 (MQ=255)
                                                                                                                    gcaacgaaTGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTa                           <  2:426976/137‑1 (MQ=255)
                                                                                                                     cAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGAGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATCGCGGCGTGGGCGGGCTTGCCTTTAc                          >  1:289838/1‑143 (MQ=255)
                                                                                                                          aaTGGTCGAGTGGCTGCGCGAGTGGAACATTCCGCTGTTAACTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTActcctc                     <  1:200771/143‑1 (MQ=255)
                                                                                                                                cGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTGTACTCCTCTGGCTc               >  1:469863/1‑143 (MQ=255)
                                                                                                                                        tgcGCGAGTGGAAAATTCCGCTGTCACTTCATCCACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTACTCCTCTGGCTCAGGGATTg       <  1:260985/143‑1 (MQ=255)
                                                                                                                                            cgAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTACTCCTCTGGCTCAGGGATTGCTGa   <  2:73317/143‑1 (MQ=255)
                                                                                                                                             gAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATATGTGGCAAAATAACGGCGTGGGCTGTATTGCCTTTACTCCTCTGGCCCAGGGCTTGCTgcc  >  2:500806/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TCTCATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTAAGGCGCTGTATGTCGGGATCTCCTCTTACTCGCCAGAGCGGACGCAAAAAATGGTCGAGTTGCTGCGCGAGTGGAAAATTCCGCTGTTAATTCATCAACCTTCGTACAATTTACTGAACCGCTGGGTGGATAAAAGCGGCCTGCTGGATACCCTGCAAAATAACGGCGTGGGCTGTATTGCCTTTACTCCTCTGGCTCAGGGATTGCTGAC  >  NC_000913/3148305‑3148588

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: