Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,233,504 T→A V470D (GTT→GAT)  prkC → non‑specific serine/threonine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,233,5040TA96.5% 91.7 / NA 29V470D (GTT→GAT) prkCnon‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (12/16);  minor base C (1/0);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATC  >  CP000730/1233370‑1233618
                                                                                                                                      |                                                                                                                  
tGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                                              <  1:852086/141‑1 (MQ=255)
 gAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAtt                                                                                                             <  1:649173/141‑1 (MQ=255)
      tgtTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                                                        >  1:116846/1‑141 (MQ=255)
         tGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGt                                                                                                     <  1:745057/141‑1 (MQ=255)
           aaCGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtat                                                                                                   <  1:109891/141‑1 (MQ=255)
                    aCAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATGAAAAGCGTTAGGTCTTAAAGATGCTACGATTGAAAAAGTATATAATAATc                                                                                          >  1:255248/1‑141 (MQ=255)
                        tgtTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAgc                                                                                      <  1:36251/141‑1 (MQ=255)
                            gATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCa                                                                                  >  1:856637/1‑141 (MQ=255)
                            gATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCa                                                                                  >  1:744365/1‑141 (MQ=255)
                                      tatCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACa                                                                        >  1:838174/1‑141 (MQ=255)
                                           aaGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATACTCAAGCGCCAAAAGGATACATTGCa                                                                   <  1:777151/141‑1 (MQ=255)
                                            aGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCaa                                                                  >  1:892740/1‑141 (MQ=255)
                                                      aaaGGTTAAAAAGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAgtgt                                                        <  1:638227/141‑1 (MQ=255)
                                                               aaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaaa                                               >  1:414267/1‑141 (MQ=255)
                                                                  gCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATAc                                            <  1:298016/141‑1 (MQ=255)
                                                                    cAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATATGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                          <  1:611789/141‑1 (MQ=255)
                                                                    cAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                          >  1:488068/1‑141 (MQ=255)
                                                                    cAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGGAACCGCAAATACTg                                          <  1:455958/141‑1 (MQ=255)
                                                                      aaTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaa                                        >  1:384988/1‑141 (MQ=255)
                                                                       aTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaaa                                       <  1:216195/141‑1 (MQ=255)
                                                                       aTGTCATTACTTGACCGAAGGAGGAAGCCTTGAAAAAAATAAAATCGCTAGGTCTTAAAGATGATCCGAATGAAAAAGTATATAAGAATCAATCGCCAAAAGGATACATTGCAAATCAAAGTGAAACCGCAAATACTGaaa                                       <  1:321488/141‑1 (MQ=255)
                                                                              tGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                <  1:275798/141‑1 (MQ=255)
                                                                              tGGGTTACCTAAGCAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                <  1:276177/141‑1 (MQ=255)
                                                                                         aggaggAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTa                     >  1:523092/1‑141 (MQ=255)
                                                                                          ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTaa                    <  1:815117/141‑1 (MQ=255)
                                                                                           gaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAGAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCCTGATTCTAat                   >  1:886811/1‑141 (MQ=255)
                                                                                                 gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaaa             >  1:575512/1‑141 (MQ=255)
                                                                                                 gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaaa             >  1:575230/1‑141 (MQ=255)
                                                                                                            aTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATc  <  1:839935/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                  
TGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATC  >  CP000730/1233370‑1233618

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