Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,042,365 C→A 20.1% L137I (CTC→ATC) ‡ PP_3563 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,042,3650CA20.1% 36.4 / 5.1 20L137I (CTC→ATC) ‡PP_3563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/11);  new base A (2/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGT  >  NC_002947/4042219‑4042510
                                                                                                                                                  |                                                                                                                                                 
tCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCa                                                                                                                                                <  2:548197/150‑1 (MQ=255)
    ccTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACTGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCGCGGCCGCTGAGCAGCGgc                                                                                                                                           <  1:526939/150‑1 (MQ=255)
                         acgcGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGc                                                                                                                       <  1:298382/150‑1 (MQ=255)
                             gccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAg                                                                                                                   <  2:450998/150‑1 (MQ=255)
                              ccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                  <  1:523199/150‑1 (MQ=255)
                                     gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGaa                                                                                                           <  1:452768/150‑1 (MQ=255)
                                      catgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAg                                                                                                          <  2:536335/150‑1 (MQ=255)
                                                aGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGGCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCagg                                                                                                >  2:112077/1‑150 (MQ=255)
                                                        ttCTGCCAGGTGGTGCTCGACCCGCAATTCGTAGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAc                                                                                        >  1:63359/1‑150 (MQ=255)
                                                                                      gtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGc                                                                                                                       >  2:984/1‑89 (MQ=255)
                                                                                      gtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGc                                                                                                                       <  1:984/89‑1 (MQ=255)
                                                                                          atatCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGACACGGCAATGCGGCGGcc                                                      >  2:34291/1‑150 (MQ=255)
                                                                                                     aGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgcc                                           <  2:161705/150‑1 (MQ=255)
                                                                                                               gCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGTGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGccgccccgc                                 <  2:217784/150‑1 (MQ=255)
                                                                                                                         gcgcATCGTGCTGGCACGGCCCCTGCTCAGCGGCCCCCTGGGCAACTGGCAGGCCGAGGAACAGAAGGTCCTGCAGGAGTTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTg                       <  2:63359/150‑1 (MQ=255)
                                                                                                                             aTCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGc                   >  2:220499/1‑150 (MQ=255)
                                                                                                                                     caggcGGCCGCTGAGCAGCGGCCGTCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCt           >  1:9439/5‑150 (MQ=255)
                                                                                                                                      aggcGGCCGCTGAGCAGCGGCCGTCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                              >  2:206740/4‑66 (MQ=25)
                                                                                                                                      aggcGGCCGCTGAGCAGCGGCCGTCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                              <  1:206740/63‑1 (MQ=25)
                                                                                                                                            ccGCGGCTCACGGGCCGCCTGGGCGACTGGCAGGCCGGGGGACAGAAGGTGCTGAAGGAGGCCACCGCCGCGCGCAAGGTGCCACGGTAGTGGGGCGGCCCGCCTTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCggg    <  2:455860/150‑1 (MQ=255)
                                                                                                                                              gctgctCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGt  >  2:137783/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGT  >  NC_002947/4042219‑4042510

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