Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,473,614 C→T 100% W156* (TGG→TGA rpsG ← 30S ribosomal subunit protein S7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,473,6140CT90.9% 21.6 / ‑3.8 12W156* (TGG→TGArpsG30S ribosomal subunit protein S7
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (6/4);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCA  >  NC_000913/3473478‑3473701
                                                                                                                                        |                                                                                       
tGTTACGGTAGCGTGCGATGGGTGTTGTCCGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGGTCAATTTACGTAGCCCAAAGCGGGCTGCTTCCTGGAGGCGCCCGCCTGGTGCCTAAACCTCCGACGGGATAATCACCgc                                                                                   >  1:186641/1‑142 (MQ=255)
                                  aaTTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATATAAGTAGCCCAAAGCGGGATGCTTACTGGAAGCACCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGAGCGCGAACGCCTTGATGGATTCGGCCAtca                                                 >  2:63255/2‑141 (MQ=255)
                                        gttCCTCGTTTATCTTTTCGGCGTTCAATTTAAGGAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGaa                                           <  1:176196/142‑1 (MQ=255)
                                         ttCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAAGCAACGGTAGTGTGCGAACGACTTGATGTCTGCGACCATACGGTGAAc                                          >  1:136636/1‑143 (MQ=255)
                                         ttCCTCGTTTATCTTTTAGGCATTCAATTTATGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAACACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTCGTTGGCTTCGGCCATACGGTGAAc                                          >  2:123574/1‑143 (MQ=255)
                                                   aTCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGt                                >  1:133900/1‑143 (MQ=255)
                                                          aGGCGTTTAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTaa                         <  2:116595/143‑1 (MQ=255)
                                                               tttaaTTTAAGTAGCCTAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCTGAAGTGTTAATCAACGGTAGTGTGCGAACGCCTGGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCa                    <  1:164719/140‑1 (MQ=255)
                                                               gtCAATGTAAGTAGCCCAAAGCGGGCTGCGTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATCACCAACGGTAGTGTGCGAACGCCTTGTGGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCa                    <  1:38137/142‑1 (MQ=255)
                                                                   aTTTAAGTAGTCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAggtc                >  1:164729/1‑140 (MQ=255)
                                                                        aGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGCACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTg           >  1:101969/1‑143 (MQ=255)
                                                                                 aaGCGGGCTGCTTACTGGAAGCGCCTTCCGGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTgca  <  2:87552/143‑1 (MQ=255)
                                                                                                                                        |                                                                                       
TGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCA  >  NC_000913/3473478‑3473701

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: