Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,446,939 C→T 15.4% Y70Y (TAC→TAT PP_3066 → lytic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,446,9390CT15.4% 50.3 / 4.1 26Y70Y (TAC→TATPP_3066lytic enzyme
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base T (2/2);  total (12/14)
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.84e-01

TTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGTGACGGTCGCCT  >  NC_002947/3446791‑3447085
                                                                                                                                                    |                                                                                                                                                  
tttATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAATTACCTGTCGAAGTACg                                                                                                                                                   >  1:199558/1‑150 (MQ=255)
    tGTTCCCGGACTCAACGGTACTATGGGGAAGTTCGCCATCATCCCGCGCCTGCGTCTGGCCGCGTTTCTCGCCCAGATCGGTCATGAGGCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGacac                                                                                                                                               <  1:245649/150‑1 (MQ=255)
     gTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACAcc                                                                                                                                              >  2:117602/1‑150 (MQ=255)
         ccGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACAccggcc                                                                                                                                          >  2:209194/1‑150 (MQ=255)
                             gggAAGTTCGCCATCATCACGCGCCTGCGGATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGc                                                                                                                      <  1:130028/150‑1 (MQ=255)
                                                    cccgCGTATGGCGGCGTTTGTCGCCCAGTTCGGTCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGAGAGCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGa                                                                                               <  2:174733/147‑1 (MQ=255)
                                                     cTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGAt                                                                                              <  1:117602/150‑1 (MQ=255)
                                                         gTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTc                                                                                          <  2:199558/150‑1 (MQ=255)
                                                                            ccAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTAtcgtggtcgtgg                                                                       <  2:139531/150‑1 (MQ=255)
                                                                            ccAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTAtcgtggtcgtgg                                                                       >  1:565450/1‑150 (MQ=255)
                                                                              aGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGc                                                                     >  2:240795/1‑150 (MQ=255)
                                                                                    ggCATGAGTCGGGTCAGTTGCAGTACGTGCGTGCGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACGGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCTAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCGTATcc                                                               >  1:424984/1‑150 (MQ=255)
                                                                                          aGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGa                                                         <  1:209194/150‑1 (MQ=255)
                                                                                                gTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGcc                                                   >  2:562500/1‑150 (MQ=255)
                                                                                                                        ttGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGTTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCcgcg                           >  2:229371/1‑150 (MQ=255)
                                                                                                                              aTGACAAGTACCTGTCGAAGTATGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAgg                                                            <  1:183923/111‑1 (MQ=255)
                                                                                                                              aTGACAAGTACCTGTCGAAGTATGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAgg                                                            >  2:183923/1‑111 (MQ=255)
                                                                                                                              aTGACAAGTACCTGTCGAAGTATGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAgg                                                            >  2:183925/1‑111 (MQ=255)
                                                                                                                              aTGACAAGTACCTGTCGAAGTATGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAgg                                                            <  1:183925/111‑1 (MQ=255)
                                                                                                                                      tACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGt             <  2:205260/150‑1 (MQ=255)
                                                                                                                                        ggtGTCGTACTTCGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAgg                                    <  1:54111/123‑1 (MQ=255)
                                                                                                                                             cGAAGTACAACACCGGCTGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGTGACGGTc      <  2:73429/150‑1 (MQ=255)
                                                                                                                                              gAAGTACGACACCTGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGTCTTATCCAGGTGACGAGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGTGACGGTCg     >  1:326450/1‑150 (MQ=255)
                                                                                                                                                 gTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGa                                      <  1:405457/114‑1 (MQ=255)
                                                                                                                                                 gTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGa                                      >  2:405457/1‑114 (MQ=255)
                                                                                                                                                 gTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGTGACGGTCGCCt  <  1:19406/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGCCGCGTTTCTCGCCCAGATCGGGCATGAGTCGGGTCAGTTGCAGTACGTGCGTGAGCTTGGCAATGACAAGTACCTGTCGAAGTACGACACCGGCCGGCTGGCTGAGCGCCTGGGCAACACGCCCGAAGCAGATGGCGATGGTCAGCTCTATCGTGGTCGTGGGCTTATCCAGGTGACGGGCCGCTTCAATTACGAGGCCTGCAGCGAGGCTCTGTTCGGTGACGGTCGCCT  >  NC_002947/3446791‑3447085

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: