Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,328,502 G→A 9.1% G18D (GGC→GAC)  PP_2927 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,328,5020GA 9.1% 120.0 / 3.2 44G18D (GGC→GAC) PP_2927hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

AACGAAAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAAT  >  NC_002947/3328366‑3328644
                                                                                                                                        |                                                                                                                                              
aacgaaAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGAc                                                                                                                                    >  8:90069/1‑149 (MQ=255)
     aaaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTgg                                                                                                                               >  3:114650/1‑149 (MQ=255)
     aaaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTgg                                                                                                                               >  1:313747/1‑149 (MQ=255)
      aaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGc                                                                                                                              <  5:231297/149‑1 (MQ=255)
      aaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGc                                                                                                                              <  6:321560/149‑1 (MQ=255)
      aaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGc                                                                                                                              >  6:231297/1‑149 (MQ=255)
             gAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATg                                                                                                                         <  4:287364/147‑1 (MQ=255)
                       ccGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGc                                                                                                              >  5:282543/1‑148 (MQ=255)
                       ccGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGc                                                                                                              >  1:194630/1‑148 (MQ=255)
                        cGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGa                                                                                                            <  7:90069/149‑1 (MQ=255)
                                cGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTgg                                                                                                    >  5:90157/1‑149 (MQ=255)
                                  gAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTgggg                                                                                                  <  4:200931/149‑1 (MQ=255)
                                            cccTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCggg                                                                                        >  1:231405/1‑149 (MQ=255)
                                             ccTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGc                                                                                       >  2:256375/1‑149 (MQ=255)
                                                       gATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGa                                                                              >  8:51687/1‑148 (MQ=255)
                                                        aTGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATTGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACa                                                                            <  7:263486/149‑1 (MQ=255)
                                                         tGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAg                                                                           <  4:208273/149‑1 (MQ=255)
                                                            ccTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGc                                                                          >  2:186176/1‑147 (MQ=255)
                                                             cttcttGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAg                                                                       >  3:279260/1‑149 (MQ=255)
                                                                       aggaggGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCAccc                                                             <  4:114650/149‑1 (MQ=255)
                                                                        ggaggGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCAccc                                                             <  6:90157/148‑1 (MQ=255)
                                                                          aggGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTg                                                          <  1:13020/149‑1 (MQ=255)
                                                                           gggCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCgg                                                                                         >  5:118132/1‑117 (MQ=255)
                                                                           gggCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCgg                                                                                         <  6:118132/117‑1 (MQ=255)
                                                                           gggCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGa                                                         <  2:231405/149‑1 (MQ=255)
                                                                                ggTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAAc                                                    <  4:279260/149‑1 (MQ=255)
                                                                                 gTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAgcggc                                                                                                                                       <  2:173957/65‑1 (MQ=255)
                                                                                 gTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAgcggc                                                                                                                                       >  1:173957/1‑65 (MQ=255)
                                                                                    aTGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGc                                                >  6:115422/1‑149 (MQ=255)
                                                                                    aTGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCCACAGCCAGGCGGTCACCCCTGAGCAACTGGc                                                <  6:282543/149‑1 (MQ=255)
                                                                                       gACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCgg                                              >  1:140853/1‑148 (MQ=255)
                                                                                        aCGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGaa                                            <  1:186176/149‑1 (MQ=255)
                                                                                            ggCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGACGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGa                                        >  7:34217/1‑149 (MQ=255)
                                                                                            ggCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGACGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGa                                        >  7:39149/1‑149 (MQ=255)
                                                                                             gCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGAc                                       <  5:115422/149‑1 (MQ=255)
                                                                                                    ttCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGAGATGAACTGGGGCTGTACCGGGCCATGGCCGCCAGCCAGGCGGTCACCCATGAGCAACTGGCGGAAAGGACCGGCTGc                                >  6:55566/1‑149 (MQ=255)
                                                                                                      cACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCAccc                                                             >  4:174147/1‑118 (MQ=255)
                                                                                                      cACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCAccc                                                             <  3:174147/118‑1 (MQ=255)
                                                                                                       aCGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCa                              <  7:51687/148‑1 (MQ=255)
                                                                                                                    cAAGCTGGTGGGCGACATGGGCACAGCGGAGACGCTGGCCAATGTCATTCTCGGCGAAGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGt                  >  6:321520/1‑147 (MQ=255)
                                                                                                                            tggGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATgg        <  5:310028/149‑1 (MQ=255)
                                                                                                                                  aCATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAAt  <  8:194559/149‑1 (MQ=255)
                                                                                                                                  aCATGGACGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAAt  <  8:34217/149‑1 (MQ=255)
                                                                                                                                  aCATGGACGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAAt  <  8:39149/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
AACGAAAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTGGCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAAT  >  NC_002947/3328366‑3328644

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