Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,812 C→A 10.9% intergenic (‑204/+25) PP_2010 ← / ← PP_2011 cytochrome b561/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,8120CA10.9% 126.4 / 7.3 46intergenic (‑204/+25)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/15);  new base A (3/2);  total (29/17)
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.65e-02

AAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAG  >  NC_002947/2282706‑2282959
                                                                                                          |                                                                                                                                                    
aaaaTGTGTTCCCCGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCCCAGTCCAACTGCACAGTCTTGGATGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGAt                                                                                                            <  6:156052/148‑1 (MQ=255)
                                     gCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcc                                                                      >  8:255228/1‑149 (MQ=255)
                                             ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                                                 >  8:71959/1‑67 (MQ=255)
                                             ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                                                 <  7:71959/67‑1 (MQ=255)
                                             ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCATCAGAAACACCCGTTCTGCCTCCTGc                                                                >  6:225929/1‑147 (MQ=255)
                                             ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcct                                                                  <  4:181527/145‑1 (MQ=255)
                                             ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcct                                                                  >  3:181527/1‑145 (MQ=255)
                                             ccACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGCAGGCGACATCAGTTACCGGGCAGATAGGGc                                                                                                       <  5:7477/108‑1 (MQ=14)
                                             ccACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGCAGGCGACATCAGTTACCGGGCAGATAGGGc                                                                                                       >  6:7477/1‑108 (MQ=14)
                                             ccACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACGGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTTCCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGc                                                                >  5:257569/1‑147 (MQ=25)
                                                 aGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTg                                                                 >  7:46766/1‑142 (MQ=255)
                                                 aGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTg                                                                 <  8:46766/142‑1 (MQ=255)
                                                    ccAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCa                                                                                                                          >  2:127469/1‑82 (MQ=14)
                                                    ccAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCa                                                                                                                          <  1:127469/82‑1 (MQ=14)
                                                      aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCg                                                                             <  2:5562/125‑1 (MQ=21)
                                                       aCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                    >  1:24695/1‑149 (MQ=255)
                                                        cTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGGGGGAGCGGGCGTTCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACAAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCgg                                                   >  3:61716/1‑149 (MQ=37)
                                                         tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCCTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCACGCACCATCgg                                                   >  7:237529/1‑148 (MQ=37)
                                                          gCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                                                 >  7:83760/1‑149 (MQ=37)
                                                           cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCATCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                                                <  5:225929/149‑1 (MQ=255)
                                                           cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                                                >  3:51410/1‑149 (MQ=255)
                                                           cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                                                >  1:137978/1‑149 (MQ=255)
                                                             caGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                                                >  4:141773/1‑147 (MQ=255)
                                                                  ttGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAc                                         >  8:127063/1‑149 (MQ=255)
                                                                    ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAc                                         >  3:156073/1‑147 (MQ=255)
                                                                    ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAc                                         <  2:24695/147‑1 (MQ=255)
                                                                        gCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCg                                   >  3:160316/1‑149 (MQ=255)
                                                                            gtgCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTAc                                >  7:18730/1‑148 (MQ=255)
                                                                              gCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAg                             >  7:301205/1‑149 (MQ=255)
                                                                               cAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAg                             >  8:182147/1‑148 (MQ=255)
                                                                                    ggACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTg                       >  7:95477/1‑149 (MQ=255)
                                                                                       cTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTgcgc                    <  8:95477/149‑1 (MQ=255)
                                                                                       cTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGACACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTgcg                     >  7:70241/1‑149 (MQ=38)
                                                                                         gtgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCg                  >  1:132772/1‑149 (MQ=255)
                                                                                          tgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCgg                 >  4:90193/1‑149 (MQ=255)
                                                                                          tgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCg                  <  2:137978/148‑1 (MQ=255)
                                                                                            ttcgCGGGCACGCCCGCGAAGGGGACGCAGCAGGCGACATCAGATACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCACCTGCCAGGAAACATCGCCATCAGGCACGAGTc                                    >  7:201247/4‑128 (MQ=255)
                                                                                             gAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGt                        >  4:305438/1‑139 (MQ=255)
                                                                                             gAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGt                        <  3:305438/139‑1 (MQ=255)
                                                                                                cGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCctg           <  7:255228/149‑1 (MQ=255)
                                                                                                cGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCctg           <  8:18730/149‑1 (MQ=255)
                                                                                                  ggCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                                                         >  2:295161/1‑85 (MQ=255)
                                                                                                  ggCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                                                         <  1:295161/85‑1 (MQ=255)
                                                                                                       gCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTcc    <  6:257569/149‑1 (MQ=255)
                                                                                                       gCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTc     >  7:59851/1‑148 (MQ=255)
                                                                                                       gCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTc     <  8:59851/148‑1 (MQ=255)
                                                                                                        cccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCa   <  3:90193/149‑1 (MQ=255)
                                                                                                         ccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAg  >  8:129037/1‑149 (MQ=255)
                                                                                                          |                                                                                                                                                    
AAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAG  >  NC_002947/2282706‑2282959

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