Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,390 G→A 7.7% G171R (GGG→AGG)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,3900GA 7.7% 137.0 / 3.0 52G171R (GGG→AGG) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base A (2/2);  total (28/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 = 3.14e-01

CAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTC  >  NC_002947/4301242‑4301538
                                                                                                                                                    |                                                                                                                                                    
cAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACgg                                                                                                                                                     <  1:1138889/150‑1 (MQ=255)
 aaCATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACggg                                                                                                                                                    <  2:433710/150‑1 (MQ=255)
    aTCACCCTGGGGTCTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGAAGTGCCACCGGCGGCGACGTTCAACATCAGCCTTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTa                                                                                                                                                 >  2:81028/1‑150 (MQ=255)
      cACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGCCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATc                                                                                                                                               <  2:874390/150‑1 (MQ=255)
         ccTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGTATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGTAAGGGGGTATCCGc                                                                                                                                            >  1:1151415/1‑150 (MQ=255)
         ccTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGACGTCTTCAAGCGGCATGAAGTGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGc                                                                                                                                            >  1:209847/1‑150 (MQ=255)
                  tCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCggg                                                                                                                                   <  2:209847/150‑1 (MQ=255)
                     gcgcACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACCAACCTATGCACACCACGTCTGAAAGAGCACAAGTCTTAACTGCCGTATGCAACGGGGTATCCGCTTCGTCGGGCGa                                                                                                                                >  1:870868/1‑150 (MQ=255)
                     gcgcACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACAGGGTATCCGCTTCGTCGGGCGa                                                                                                                                >  1:406644/1‑150 (MQ=255)
                         aCCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTg                                                                                                                            >  2:814337/1‑150 (MQ=255)
                            ttCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAg                                                                                                                         >  1:201138/1‑150 (MQ=255)
                              cAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCt                                                                                                                       <  2:201138/150‑1 (MQ=255)
                                 ccAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGa                                                                                                                    >  2:492821/1‑150 (MQ=255)
                                         accCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCt                                                                                                            >  1:75549/1‑150 (MQ=255)
                                                    tCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGc                                                                                                 >  1:623074/1‑150 (MQ=255)
                                                      cacaGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCt                                                                                               >  2:475836/1‑150 (MQ=255)
                                                       acaGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTAACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTg                                                                                              >  1:920318/1‑150 (MQ=255)
                                                          gTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGtta                                                                                           <  1:475836/150‑2 (MQ=255)
                                                          gTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTc                                                                                           <  1:492821/150‑1 (MQ=255)
                                                                aaGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTc                                                                                     >  2:207339/1‑150 (MQ=255)
                                                                aaGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTc                                                                                     <  2:804357/150‑1 (MQ=255)
                                                                         ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtcgtcg                                                                                                                                     <  2:2048/93‑1 (MQ=255)
                                                                         ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTtcgtcg                                                                                                                                     >  1:2048/1‑93 (MQ=255)
                                                                         ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCggg                                                                            <  2:709167/150‑1 (MQ=255)
                                                                         ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCggg                                                                            <  1:207339/150‑1 (MQ=255)
                                                                         ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACAGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCggg                                                                            >  1:191168/1‑150 (MQ=255)
                                                                          cACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCgggg                                                                           >  1:5694/1‑150 (MQ=255)
                                                                            cTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                                                               <  1:548159/96‑1 (MQ=255)
                                                                            cTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGAt                                                                                                                               >  2:548159/1‑96 (MQ=255)
                                                                             tGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCt                                                                        <  2:623074/150‑1 (MQ=255)
                                                                              gggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTa                                                                       <  2:75549/150‑1 (MQ=255)
                                                                               ggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTAc                                                                      <  1:814337/150‑1 (MQ=255)
                                                                                  gCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACAGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAAc                                                                   <  2:191168/150‑1 (MQ=255)
                                                                                     aCGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAg                                                                <  2:989571/150‑1 (MQ=255)
                                                                                        ttCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACAGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTc                                                             <  2:406644/150‑1 (MQ=255)
                                                                                               tCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCAcggcgg                                                                             <  1:1012661/127‑1 (MQ=255)
                                                                                               tCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCAcggcgg                                                                             >  2:1012661/1‑127 (MQ=255)
                                                                                               tCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTaca                                                      >  1:51849/1‑150 (MQ=255)
                                                                                                     ttCACACAACGTCTGAAGGCGTGGTAACCGTGACTGCTGAATGCAACGGGGTATCCGCTGCGTCGGGCGCTTGGCAGCTCGACAGTACCGGGGCTTCACGCCTGTTCGAGTTCCACAGCGGGGGGTACAACCAGGTCGTCTGCACCGAcc                                                >  2:1039268/1‑149 (MQ=255)
                                                                                                        acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCaa                                             <  1:201653/150‑1 (MQ=255)
                                                                                                                    aaGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCccgcc                                     <  2:399682/146‑1 (MQ=255)
                                                                                                                    aaGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCccgcc                                     >  1:399682/1‑146 (MQ=255)
                                                                                                                     aGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGa                                <  2:5694/150‑1 (MQ=255)
                                                                                                                      ggCGTGGTAACCGTAACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGAc                               <  2:920318/150‑1 (MQ=255)
                                                                                                                         gTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGa                                                                                                                    >  1:211410/1‑62 (MQ=255)
                                                                                                                         gTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGa                                                                                                                    <  2:211410/62‑1 (MQ=255)
                                                                                                                           ggTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTAGCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTATAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTc                          >  2:868833/1‑150 (MQ=255)
                                                                                                                                  gTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGc                   >  2:543447/1‑150 (MQ=255)
                                                                                                                                    gACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGcat                 >  2:21774/1‑150 (MQ=255)
                                                                                                                                            aaTGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCg         >  1:425109/1‑150 (MQ=255)
                                                                                                                                             aTGCAGTGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTCAAAAGTCCGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGa        <  1:81028/150‑1 (MQ=255)
                                                                                                                                                 aaCGGGGTATCCGGTTCGTCGGGCGCGTTGCAGGTTGACAGTACCTGCGCGTCACGCCTGTTCGAGTTCCCCGGCGTGGATTTCAACCAGCTCGACTACACCGACACAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGctc    >  2:799542/1‑150 (MQ=255)
                                                                                                                                                 aaCGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGctc    >  2:472633/1‑150 (MQ=255)
                                                                                                                                                   cGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCtctc  <  1:472633/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACATCACCCTGGGGTTTCGGCGCACCTTCAACCAGACCGACCCGGTGAATTCCACAGTCCTGAAGGGGGTGCCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTC  >  NC_002947/4301242‑4301538

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