Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,443,461 C→T G213S (GGC→AGC)  secY ← Sec translocon subunit SecY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,4610CT80.0% 12.6 / ‑2.5 10G213S (GGC→AGC) secYSec translocon subunit SecY
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (7/1);  total (8/2)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATT  >  NC_000913/3443238‑3443677
                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
gagAACAAAATAATACTGAAAGCGAATATGGCCGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCTCTTTGACGCGCTTGCTCGATAGTATGGGc                                                                                                                                                                                               <  1:269537/249‑1 (MQ=255)
          taataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAACTGTGTGCTCTGTGCAGCATAGACACGACGACATTGCTGACGTTTCGCGTAGTTTACCACAATGAGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCCtt                                                                                                                                                                                                                        >  1:298807/1‑216 (MQ=255)
                                                                             ctGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGACGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTa                                                                                                                   >  1:10298/1‑250 (MQ=255)
                                                                                                                          tACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTccc                                                                                                                                                                            >  1:9354/1‑148 (MQ=255)
                                                                                                                             cacaATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACtt                                                                   >  1:23122/1‑250 (MQ=255)
                                                                                                                                                  cTCAACAAATACAACAAAGAACGTCACTGCAAATACTACTACTGCAACCAACAGAAACACGAGGAAGTGCAGGTAGCTTTGACGCGCTTGCTCGATCGTATggg                                                                                                                                                                                                >  1:153749/1‑104 (MQ=255)
                                                                                                                                                                                     ctaatactGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATCGTATGGGCAATAGCTGGCGGGAGACCCGCGACACTCCCGGCGAAGATAATGATTGCAATACCGTTGCCGATACATCGTTCAGTAATTTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCCGGCCTTGCATAcc           >  1:319534/1‑250 (MQ=255)
                                                                                                                                                                                             ncaaccaacAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATAtt  >  1:81807/2‑251 (MQ=255)
                                                                                                                                                                                                 tcaacAGCAACACGAGGAAGTGCAGGTCGCTTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCNNCGAAGa                                                                                                                                                         >  2:3968/2‑96 (MQ=255)
                                                                                                                                                                                                                     tGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGGTCGCCCAACCGCATCAGGAACATGGTTCCTGGGACCAGACTTACAACAGCg                                                          <  1:277072/171‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
GGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTGCAACCAACAGCAACACGAGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATT  >  NC_000913/3443238‑3443677

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