Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,161,581 G→A 8.7% V117I (GTC→ATC)  PP_4544 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,161,5810GA 8.7% 131.1 / 2.4 46V117I (GTC→ATC) PP_4544hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (2/2);  total (22/24)
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.61e-01

AACGTATGAGCTCCTTTGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTGGGTCGACCTGGCCTACCTC  >  NC_002947/5161435‑5161724
                                                                                                                                                  |                                                                                                                                               
aaCGTATGAGCTCCTTTGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTATc                                                                                                                                               >  1:184884/1‑149 (MQ=255)
aaCGTATGAGCTCCTTTGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTATc                                                                                                                                               >  1:189644/1‑149 (MQ=255)
             cTTTGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAg                                                                                                                                    <  5:296735/147‑1 (MQ=255)
              tttGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGc                                                                                                                                   >  3:133033/1‑147 (MQ=255)
                      aaaaGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTCCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAAcc                                                                                                                         >  5:228889/1‑149 (MQ=255)
                         aGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGc                                                                                                                      <  4:133033/149‑1 (MQ=255)
                         aGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGc                                                                                                                      <  3:132460/149‑1 (MQ=255)
                                tGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGa                                                                                                               <  8:165288/149‑1 (MQ=255)
                                   gccgccAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGc                                                                                                            >  7:272680/1‑149 (MQ=255)
                                       ccAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTggg                                                                                                        >  8:160112/1‑149 (MQ=255)
                                         aaGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATGCGAACAGCAGCAGCGGCAACCAGGTCCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCCGGTCGACGCTCGGGt                                                                                                      >  2:201196/1‑149 (MQ=255)
                                          aGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTATCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTc                                                                                                     <  2:184884/149‑1 (MQ=255)
                                          aGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTATCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTc                                                                                                     <  2:189644/149‑1 (MQ=255)
                                              ccAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTa                                                                                                 >  4:271873/1‑149 (MQ=255)
                                              ccAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTa                                                                                                 >  2:248688/1‑149 (MQ=255)
                                               cAGAACAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTAt                                                                                                >  3:274383/1‑149 (MQ=255)
                                                        ggTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTgg                                                                                       <  1:248688/149‑1 (MQ=255)
                                                            aCGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGgcagc                                                                                   <  7:160112/149‑1 (MQ=255)
                                                            aCGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGgcagc                                                                                   <  7:95440/149‑1 (MQ=255)
                                                               tGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCagcagc                                                                                <  4:274383/149‑1 (MQ=255)
                                                               tGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCagaagc                                                                                >  6:179170/1‑149 (MQ=255)
                                                                     tGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGAc                                                                          >  6:310495/1‑149 (MQ=255)
                                                                       cGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGAc                                                                          <  8:272680/147‑1 (MQ=255)
                                                                       cGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGAc                                                                          <  3:271873/147‑1 (MQ=255)
                                                                       cGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACATGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGAc                                                                          <  1:126838/147‑1 (MQ=255)
                                                                               aGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGa                                                                <  6:34293/149‑1 (MQ=255)
                                                                                 aCGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAgcg                                                                               <  1:135340/132‑1 (MQ=255)
                                                                                 aCGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAgcg                                                                               >  2:135340/1‑132 (MQ=255)
                                                                                   gTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGaccaa                                                            >  2:259584/1‑149 (MQ=255)
                                                                                               gaTCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGa                                                <  5:310495/149‑1 (MQ=255)
                                                                                               gaTCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCAACAGCGAGATCAAGCCGCAGGTCGACGCTGGGGTCCGTATCGCGACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGa                                                <  1:201196/149‑1 (MQ=255)
                                                                                                     aaCAGAAGCAGCGGCAA‑CAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGCCAactac                                         >  5:229355/1‑149 (MQ=255)
                                                                                                            gcagcGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGa                                    >  3:304940/1‑148 (MQ=255)
                                                                                                                cGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACaaga                               <  4:304940/149‑1 (MQ=255)
                                                                                                                cGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACaag                                >  6:15514/1‑148 (MQ=255)
                                                                                                                  gCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGAc                                                                          <  8:251163/104‑1 (MQ=255)
                                                                                                                           aCCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTgg                    <  6:228889/149‑1 (MQ=255)
                                                                                                                            ccGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTggg                   <  5:15514/149‑1 (MQ=255)
                                                                                                                             cGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGt                                                                    >  6:177708/1‑99 (MQ=255)
                                                                                                                             cGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGt                                                                    <  5:177708/99‑1 (MQ=255)
                                                                                                                                    gAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTc                                                                  <  3:171339/94‑1 (MQ=255)
                                                                                                                                    gAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTc                                                                  >  4:171339/1‑94 (MQ=255)
                                                                                                                                    gAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTGGGTCGACCTg           >  5:200043/1‑149 (MQ=255)
                                                                                                                                            gCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGc                         <  2:205128/127‑1 (MQ=255)
                                                                                                                                            gCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGc                         >  1:205128/1‑127 (MQ=255)
                                                                                                                                              ccGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTGGGTCGACCTGGCCTACCTc  >  6:38310/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AACGTATGAGCTCCTTTGAACAAAAAGTGGCATGGGCCGCCAAGACCCAGATCAAGGGTGACGTGCGCCTGCGTTACGAAGACGTCAACGTCGACGATCCGAACAGCAGCAGCGGCAACCAGGACCGCCAGCGAGTACGTGCCCGTGTCGGCTTCTACAGCGAGATCAACCCGCAGGTCGACGCTGGGGTCCGTATCGCCACTGGCAGCAGCGCTGACCGACGTTCGACCAACCAGAGCTTCGACAACTACTTCGACAAGAAGTCGCTGTGGGTCGACCTGGCCTACCTC  >  NC_002947/5161435‑5161724

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