Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,033,610 G→A 62.5% G156S (GGT→AGT)  trkH → K(+) transporter TrkH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,033,6100GA62.5% 9.5 / 11.3 22G156S (GGT→AGT) trkHK(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (11/2);  minor base T (3/5);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 5.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGGAGCCACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCTGTGCGCTGGCGTTGTGG  >  NC_000913/4033480‑4033747
                                                                                                                                  |                                                                                                                                         
gggAGCCACTACGCTGGTGTGGCTGGACTCGCTCCCCCACGCCATCCTCTGTGATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGTGTGGGATGCAGc                                                                                                                               <  1:151836/143‑1 (MQ=255)
      cACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTGTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGCTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGGGGGATGCAGCTCTATc                                                                                                                         >  2:195254/1‑143 (MQ=255)
                      cTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGACTTGGATGCAGCTCTATCGCGCAGAAATgcccgg                                                                                                         >  1:536515/1‑143 (MQ=255)
                          aTTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATgcccggcccg                                                                                                     >  1:334849/1‑143 (MQ=255)
                                tCCTTCCCGCCATCCTATTGTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATGATCGTGTTGGCGGTTGCGATACTGCCTGTACTCGGCGTGAGTGGGATGCAGTTCTATCGCGCAGAAATGCCCGGCCCGCTGaaa                                                                                               <  1:446436/143‑1 (MQ=255)
                                    tCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTAGCCGTAAGAGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATa                                                                                           >  1:443268/1‑143 (MQ=255)
                                           aTCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGCTACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCCCAGAAATGCCCGGCCCGCTGAAAGATAACAAAAt                                                                                    >  1:403814/1‑143 (MQ=255)
                                            tCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTCTCCTCGGCGTGAGGGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGGTAACAAAATg                                                                                   >  1:524144/1‑143 (MQ=255)
                                               tctTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGTGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAGAATgcgc                                                                                >  1:150163/1‑143 (MQ=255)
                                                         ccAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGTGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGGTGAACGATTACAAAATGCGCCCGCGAAtgt                                                                      >  2:195322/1‑141 (MQ=255)
                                                          cAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGTGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGc                                                                     <  1:257591/143‑1 (MQ=255)
                                                                          ttGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAACTGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGCAACGGCGAAAAcc                                                     >  1:143611/1‑143 (MQ=255)
                                                                                   tGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGGGGTTg                                            >  1:362717/1‑143 (MQ=255)
                                                                                       gATCATCGTGTTGGCGGTTGCGATACTGCCTAACCTCACCGTGAGTGGGATGCAGCGCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGAgtt                                        >  1:158290/1‑143 (MQ=255)
                                                                                        atcatcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCt                                                                                                   >  2:140372/1‑83 (MQ=255)
                                                                                        atcatcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCt                                                                                                   <  1:140372/83‑1 (MQ=255)
                                                                                                            gATACTGCCGATCCTCGGCGTGTGTGGGGTGCAGCTCTATCGCGCAGAAACGCCCGGCCCGCTGAAAGACAACTAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCt                   <  2:46227/143‑1 (MQ=255)
                                                                                                              aaCTGTCGATCCTCGGCGTGTGGGGGATGCAGCGCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCtgt                 <  1:381807/142‑1 (MQ=255)
                                                                                                                  gCCTATCCTCGGCGTGTGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGACTACCCTGTGGTTGATTTATGTCTTGCTGACCGTCTCCTGTccgc             >  1:243005/1‑143 (MQ=255)
                                                                                                                          tCGGCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCTGTGCGCTGGCGTtg     >  1:499072/1‑143 (MQ=255)
                                                                                                                            ggCGTGTGTGGGGCGCAGCTCTAGCGCGCAGCCATGCCCGGCCCGCTGAACGATAACAAAATGCGCTCGTGAATCGCGGAAAGGGCGAAGACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCTGTGCGCTGGCGTtgtg   <  2:142561/143‑1 (MQ=255)
                                                                                                                             gCGTGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGGAACGGCGAAAACCCAGTGGTTGAAGAATGTCTTGCTGACCGTCGCCTGCGCGCTGGCGtggtgg  >  2:458912/1‑143 (MQ=255)
                                                                                                                                  |                                                                                                                                         
GGGAGCCACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTCGCCTGTGCGCTGGCGTTGTGG  >  NC_000913/4033480‑4033747

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: