Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,841 T→G 33.2% intergenic (‑16/+108) PP_3081 ← / ← PP_5560 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,8410TG33.2% 12.1 / 10.1 21intergenic (‑16/+108)PP_3081/PP_5560hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/10);  new base G (4/3);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTC  >  NC_002947/3466702‑3466989
                                                                                                                                           |                                                                                                                                                    
cGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGGCGGccgcgggc                                                                                                                                             >  8:184724/1‑147 (MQ=255)
 gggCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCAACATCGCGCCTTGGCCGCCGCAGGCGGCCa                                                                                                                                            <  2:96848/149‑1 (MQ=255)
          gTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCtgtg                                                                                                                                    >  8:208470/1‑148 (MQ=255)
              cAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAggg                                                                                                                               <  7:95697/149‑1 (MQ=255)
                          ggggTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCCGCCCACGCTGTGCAGGGCAGCAAaccagc                                                                                                                   >  4:67572/1‑149 (MQ=255)
                             gTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGaa                                                                                                                <  3:67572/149‑1 (MQ=255)
                               ccGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTTAAAGACGAGCAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACa                                                                                                              <  7:184724/149‑1 (MQ=255)
                                              aCCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGt                                                                                               <  3:207899/149‑1 (MQ=255)
                                               ccGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGGGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGGGAACAGCCAGTGTGTGAGGTg                                                                                              >  8:216063/1‑149 (MQ=255)
                                                   aggcggggTCCTGGCCGACCTGTTGCTGGGCAATGACGCCGCCATCAGTTAATGACGAGGAACTGCCCGAGATCGCGCCTTGGCGGCCTGGGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATg                                                                                          <  1:125449/144‑1 (MQ=255)
                                                         ggTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAAt                                                                                    <  1:15740/149‑1 (MQ=255)
                                                                     ccTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGc                                                                        <  7:208470/149‑1 (MQ=255)
                                                                     ccTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCAGCGGCGGCCCACGCGGTGCAGGGCAGCAAACCAGCGCACAGCCAGTGTCGGGGGTGCATGGTCAATCCTTGGCAGGGc                                                                        >  8:12862/1‑149 (MQ=255)
                                                                                          cgcaATCAGTTAAACACTAGTAACTGACCTACATCGCGCCTTGGCCGCGTGCGGGGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCa                                                   <  2:302837/149‑1 (MQ=255)
                                                                                               tccgTGAAACACTAGGAACTGACCGACATCGCGCCTGGGCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCt                                              <  2:134354/146‑1 (MQ=255)
                                                                                                           gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGCCTAGCCCAGGTCTGCTTCAGCGGCt                                  <  6:174984/149‑1 (MQ=255)
                                                                                                                 cTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTg                            <  3:147436/149‑1 (MQ=255)
                                                                                                                              agcgttGGCCGCCGCAGGCGGCCAAGGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCGGGTCCTCGGCAGTg               >  3:44836/5‑149 (MQ=255)
                                                                                                                              agcgttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGGCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTg               >  4:183970/5‑149 (MQ=255)
                                                                                                                                cgttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGcc             >  1:122157/3‑149 (MQ=255)
                                                                                                                                    ggCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc           <  3:36047/147‑1 (MQ=255)
                                                                                                                                           tGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTc  >  8:81921/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTC  >  NC_002947/3466702‑3466989

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