Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,329,204 C→A 16.6% P300P (CCC→CCA PP_2043 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,329,2040CA16.6% 73.4 / 2.9 24P300P (CCC→CCAPP_2043hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/9);  new base A (2/2);  total (13/11)
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.83e-01

GGCGCCGAACATCGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCCACATCGACCCGGCACCGACCCGCTACGAACTGCACCGGGTGTGGGTCGTGGAGGCGACG  >  NC_002947/2329069‑2329324
                                                                                                                                       |                                                                                                                        
ggCGCCGAACATCGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGccc                                                                                                                          <  1:382298/136‑1 (MQ=255)
  cgcCGAACATCGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTa                                                                                                                        <  2:123185/136‑1 (MQ=255)
        aCATCGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTg                                                                                                                  <  1:296615/136‑1 (MQ=255)
            cGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTAc                                                                                                              <  1:215952/136‑1 (MQ=255)
                     caccacTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCATACAACTGTTACCGCCTGGCa                                                                                                     >  1:400994/1‑136 (MQ=255)
                     caccacTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCATACAACTGTTACCGCCTGGCa                                                                                                     >  1:52442/1‑136 (MQ=255)
                                cTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCt                                                                                          <  1:481819/136‑1 (MQ=255)
                                cTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCATACAACTGTTACCGCCTGGCATCGCCCGAGCt                                                                                          <  2:400994/136‑1 (MQ=255)
                                cTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCATACAACTGTTACCGCCTGGCATCGCCCGAGCt                                                                                          <  2:52442/136‑1 (MQ=255)
                                     aCCGCCGGCCTGTGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCt                                                                                     >  1:53985/1‑136 (MQ=255)
                                     aCCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCt                                                                                     >  1:367796/1‑136 (MQ=255)
                                        gccgGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCCCGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTAtg                                                                                  >  2:263224/1‑136 (MQ=255)
                                          cgGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTAtgtg                                                                                >  2:129070/1‑136 (MQ=255)
                                          cgGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTAtgtg                                                                                >  2:8275/1‑136 (MQ=255)
                                          cgGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGACAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTAtgtg                                                                                <  1:150913/136‑1 (MQ=255)
                                                              ctCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCca                                                            <  1:28695/135‑1 (MQ=255)
                                                                            aaTGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTa                                                                                                                        >  1:226685/1‑62 (MQ=255)
                                                                            aaTGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTa                                                                                                                        <  2:226685/62‑1 (MQ=255)
                                                                                                         gctgctGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCt                                                                                                         >  2:408498/1‑48 (MQ=255)
                                                                                                         gctgctGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCt                                                                                                         <  1:408498/48‑1 (MQ=255)
                                                                                                          ctgctgGGCAAGCCGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCCACATCGACCCGGCACCGACCCGCTACGAACTCCACCGGGAGTggg                >  1:223208/1‑136 (MQ=255)
                                                                                                                        ggCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGccc                                                                                               <  2:98423/43‑1 (MQ=255)
                                                                                                                        ggCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGccc                                                                                               >  1:98423/1‑43 (MQ=255)
                                                                                                                        ggCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCCACATCGACCCGGCACCGACCCGCTACGAACTGCACCGGGTGTGGGTCGTGGAGGCGACg  >  1:395020/1‑136 (MQ=255)
                                                                                                                        ggCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCCACATCGACCCGGCACCGACCCGCTACGAACTGCACCGGGTGTGGGTCGTGGAGGCGACg  >  1:159072/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                        
GGCGCCGAACATCGCCTATGACACCACTGGCCCTGGCACCGCCGGCCTGCGGACCGCCGACAGTCGCGACATGTTCAATGGCGCGCCTGACCGGTATGACTGGAAGCTGCTGGGCAAGCAGGCGCTGCACGTGCCCTACAACTGTTACCGCCTGGCATCGCCCGAGCTGCGCTATGTGGACCTGATCAAGCCAGGCCACATCGACCCGGCACCGACCCGCTACGAACTGCACCGGGTGTGGGTCGTGGAGGCGACG  >  NC_002947/2329069‑2329324

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