Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,042,376 C→T 19.0% R140R (CGC→CGT PP_3563 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,042,3760CT19.0% 34.3 / 2.4 21R140R (CGC→CGTPP_3563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base T (2/2);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

ACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCGCCAACCACA  >  NC_002947/4042244‑4042524
                                                                                                                                    |                                                                                                                                                    
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                                                                                                            <  1:206740/63‑1 (MQ=25)
                                                                                                             aggcGGCCGCTGAGCAGCGGCCGTCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                                            >  2:206740/4‑66 (MQ=25)
                                                                                                                   ccGCGGCTCACGGGCCGCCTGGGCGACTGGCAGGCCGGGGGACAGAAGGTGCTGAAGGAGGCCACCGCCGCGCGCAAGGTGCCACGGTAGTGGGGCGGCCCGCCTTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCggg                  <  2:455860/150‑1 (MQ=255)
                                                                                                                     gctgctCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGt                >  2:137783/1‑150 (MQ=255)
                                                                                                                               ggctcccTGGGCGACCTGCAGCACGAGGAAGAAAGGGCGCTGCAGGAGGTCAACGCCGCGCGAAAGGTCCCAAGGCAATGCGGCGGCCGGACGTTGGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCGGGCCGGGGTCGCCGCCAAc      <  1:220499/145‑1 (MQ=255)
                                                                                                                               ggCCGTCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCGCCAAc      <  2:9439/150‑1 (MQ=255)
                                                                                                                               ggCCGCCTTGGCGACTGGCAGGCCGAGGGACAGAAGGTGAGGCAGGAGGTCAACGCGGAGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCGCCAAc      >  2:178530/1‑150 (MQ=255)
                                                                                                                                 ccgccTGGGCGACGGGCAGGCCGAGGGACAGAAGGTGCGGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGTGGTTTGCCCCCGCCCCGCCACTAAGCTGGAGCGCAGTTCTGGCCGGGGTCGCCGCCAACca    <  1:420604/150‑1 (MQ=255)
                                                                                                                                   gccTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCGCCAACcaca  >  2:539922/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                                    
ACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCGCCAACCACA  >  NC_002947/4042244‑4042524

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: