Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,177,447 G→A 10.7% E115K (GAG→AAG)  PP_5485 → GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,177,4470GA10.7% 92.3 / 2.6 38E115K (GAG→AAG) PP_5485GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base A (2/2);  total (21/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 = 6.57e-01

GTTGTCGTACCTGCCCAGCGAAAGTCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAACTCGCCTCCA  >  NC_002947/2177306‑2177586
                                                                                                                                             |                                                                                                                                           
gTTGTCGTACCTGCCCAGCGAAAGTCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCaga                                                                                                                                       <  1:210209/148‑1 (MQ=255)
gTTGTCGTACCTGCCCAGCGAAAGTCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATAAGCa                                                                                                                                         >  8:26753/1‑146 (MQ=255)
gTTGTCGTACCTGCCCAGCGAAAGTCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATAAGCa                                                                                                                                         >  8:21524/1‑146 (MQ=255)
                        tCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAgg                                                                                                               >  2:325227/1‑148 (MQ=255)
                          aGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCa                                                                                                             <  2:328300/148‑1 (MQ=255)
                                         gctgctCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCa                                                                                                                          >  6:181061/1‑120 (MQ=255)
                                               cAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTc                                                                                        >  1:89569/1‑148 (MQ=255)
                                                agagCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCaga                                                                                                                                       >  1:312265/1‑100 (MQ=255)
                                                agagCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCaga                                                                                                                                       <  2:312265/100‑1 (MQ=255)
                                                agagCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATAAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTc                                                                                        <  7:26753/147‑1 (MQ=255)
                                                agagCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATAAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTc                                                                                        <  7:21524/147‑1 (MQ=255)
                                                agagCACATGAGCGAATAGCGGAAGCTGACGGAGCCGAAAGCCTTCATCTGCTGACAACAACAGCACCCGACTTCTTCGCGGCGCGGGGATATGAGCAGAGAGACCTATCCCaa                                                                                                                         >  8:222418/1‑114 (MQ=255)
                                                             gAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGccc                                                                          >  1:13676/1‑148 (MQ=255)
                                                                  gCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTc                                                                    >  8:311581/1‑149 (MQ=255)
                                                                     gAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGc                                                                 <  7:311581/149‑1 (MQ=255)
                                                                      aaGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCg                                                                <  2:13676/149‑1 (MQ=255)
                                                                                    cAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTc                                                                    >  8:102474/1‑131 (MQ=255)
                                                                                    cAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTc                                                                    <  7:102474/131‑1 (MQ=255)
                                                                                         ccTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGc                                             >  1:101485/1‑149 (MQ=255)
                                                                                           ttCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCtt                                           >  4:254397/1‑149 (MQ=255)
                                                                                                           aacaGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGaa                                                  <  8:267863/126‑1 (MQ=255)
                                                                                                           aacaGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGaa                                                  >  7:267863/1‑126 (MQ=255)
                                                                                                            acaGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGgt                          <  3:254397/149‑1 (MQ=255)
                                                                                                                cACCCGACTTCTTCGCGGCTCGGGGATATTAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGa                                      <  5:19221/133‑1 (MQ=255)
                                                                                                                cACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGa                                      >  6:19221/1‑133 (MQ=255)
                                                                                                                   ccGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTcacca                   >  2:158874/1‑149 (MQ=255)
                                                                                                                   ccGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTcac                     <  7:172970/147‑1 (MQ=255)
                                                                                                                   ccGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTcac                     >  8:172970/1‑147 (MQ=255)
                                                                                                                     gACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTcacca                   >  2:302410/1‑147 (MQ=255)
                                                                                                                      aCTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAg                >  1:301869/1‑149 (MQ=255)
                                                                                                                        ttcttcGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGAc                                                                       <  2:202560/92‑1 (MQ=255)
                                                                                                                        ttcttcGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGAc                                                                       >  1:202560/1‑92 (MQ=255)
                                                                                                                          cttcGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGt              <  1:325227/147‑1 (MQ=255)
                                                                                                                          cttcGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTCAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTaa            >  2:339335/1‑149 (MQ=255)
                                                                                                                           ttcGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTAGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAAc           >  2:32003/1‑149 (MQ=255)
                                                                                                                              gcgGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAACTc         <  2:301869/148‑1 (MQ=255)
                                                                                                                               cgGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAACTCGc       <  2:89569/149‑1 (MQ=255)
                                                                                                                                       cgATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCg                            <  6:194795/119‑1 (MQ=255)
                                                                                                                                             gAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAACTCGCCTCCa  <  4:87290/140‑1 (MQ=255)
                                                                                                                                             gAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGAGTCACCACAGGTAACTCGCCTCCa  >  3:87290/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                           
GTTGTCGTACCTGCCCAGCGAAAGTCAGGGATCGGTAACCTGCTGCTCAGAGCACTTGAGCGAATAGCGGAAGCTGACGGAGCCCAAAGCCTTCATCTGCTGACCACAACAGCACCCGACTTCTTCGCGGCTCGGGGATATGAGCAGAGAGACCGATCCCAAGCCCCAAAGGCAATTTCCAGTTCGGAGGAGTTCAAAAGCCTCTGCCCGACGTCTGCGGTTTACATGGTGAAGGAGCTTGGCGACATAAACTCGGTGTCACCACAGGTAACTCGCCTCCA  >  NC_002947/2177306‑2177586

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: ‑