Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,135,915 T→C L349P (CTA→CCA)  SAUSA300_RS10870 → potassium transporter KtrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,135,9150TC96.2% 61.7 / ‑5.2 26L349P (CTA→CCA) SAUSA300_RS10870potassium transporter KtrB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/2);  minor base G (1/0);  total (24/2)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATATACGTAAAGAAAATAATGTTTCAGTATTCAATAAAGAAATATCTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCTATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAC  >  NC_007793/2135821‑2136008
                                                                                              |                                                                                             
acacaccGTAAAGAAAATAATGTTTCAGTATTCAATAAAGAAATATCTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTa                                                                                          >  1:2044963/7‑100 (MQ=255)
  tataCGTAAAGAAAATAATGTTTCAGTATTCAATAAAGAAATATCTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATc                                                                                        <  1:2583983/100‑1 (MQ=255)
                         cAGTATTCAATAAAGAAATATCTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATCGGCATCATTACTTTTATATTATc                                                                 >  1:6342542/1‑100 (MQ=255)
                                              cTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCGATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATc                                            >  1:3263595/1‑100 (MQ=255)
                                                gACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCAt                                          >  1:5300214/1‑100 (MQ=255)
                                                gACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCAt                                          >  1:4192015/1‑100 (MQ=255)
                                                gACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCAt                                          >  1:343145/1‑100 (MQ=255)
                                                  gccaCATATAAAACTATCTATTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTa                                        >  1:1670936/4‑100 (MQ=255)
                                                                     atTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCGAGTTATTATTCGAAg                     >  1:4225623/1‑100 (MQ=255)
                                                                     atTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAg                     >  1:5804021/1‑100 (MQ=255)
                                                                     atTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAg                     >  1:1948652/1‑100 (MQ=255)
                                                                     atTGTTACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAg                     >  1:2022827/1‑100 (MQ=255)
                                                                        gTTGCCATTAATATCTCATTTCCATTTATCAGTATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGCTATTATTCGAAGTgg                  >  1:6644006/1‑100 (MQ=255)
                                                                          tACCATTAATATCTCATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGtt                >  1:5719796/1‑100 (MQ=255)
                                                                                tAATATCTCATTTCCATTTATCAGCGTCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCa          >  1:1653836/1‑100 (MQ=255)
                                                                                  cgatCTCATTCCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCAtt        <  1:4601157/98‑1 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTTTTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:3033721/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:2911685/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:2455616/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:2297088/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:4411650/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:4430428/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:1860734/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:6283594/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:1549874/1‑100 (MQ=255)
                                                                                        cATTTCCATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAc  >  1:1531183/1‑100 (MQ=255)
                                                                                              |                                                                                             
TATATACGTAAAGAAAATAATGTTTCAGTATTCAATAAAGAAATATCTGACAAACATATAAAACTATCTATTGTTACCATTAATATCTCATTTCTATTTATCAGCATCATTACTTTTATATTATCGATAATTAATCCGAACATATCATTAATCAAGTTATTATTCGAAGTGGTTTCTGCATTCGGAAC  >  NC_007793/2135821‑2136008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑