Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,115,751 T→G 10.4% H28P (CAC→CCC)  PP_4500 ← TrpH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,115,7510TG10.4% 118.0 / 7.4 48H28P (CAC→CCC) PP_4500TrpH family protein
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base G (2/3);  total (24/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.16e-01

ACCCAGCACATGGATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGC  >  NC_002947/5115603‑5115896
                                                                                                                                                     |                                                                                                                                                 
aCCCAGCACATGGATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTg                                                                                                                                                  >  2:37536/1‑150 (MQ=255)
  ccAGCACATGGATGGTCGCGCCGCCCCACGTGCCGGACAGCTCCACCCCGCGTACCCAGAGAATACCCAGGTCCGAGGACGAAGTGATGC‑CGGCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCt                                                                                                                                                >  1:810013/1‑150 (MQ=255)
            gatggTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCg                                                                                                                                      <  1:335123/150‑1 (MQ=255)
             atggTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCgg                                                                                                                                     <  1:774354/150‑1 (MQ=255)
               ggTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACgggg                                                                                                                                   <  1:711520/150‑2 (MQ=255)
                gTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggggt                                                                                                                                  <  1:468883/150‑3 (MQ=255)
                  cgcgCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc                                                                                                                                <  2:523806/150‑1 (MQ=255)
                   gcgcCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggg                                                                                                                                    >  2:724463/1‑145 (MQ=255)
                   gcgcCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggg                                                                                                                                    <  1:724463/145‑1 (MQ=255)
                             aCGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAg                                                                                                                     >  1:584777/1‑150 (MQ=255)
                                              aCCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc                                                                                                                                <  2:889274/122‑1 (MQ=255)
                                              aCCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc                                                                                                                                >  1:889274/1‑122 (MQ=255)
                                                 ccGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAg                                                                                                 <  2:584777/150‑1 (MQ=255)
                                                         ccAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGcc                                                                                                        <  2:312070/135‑1 (MQ=255)
                                                         ccAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGcc                                                                                                        >  1:312070/1‑135 (MQ=255)
                                                              gcATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGAc                                                                                    >  1:196095/1‑150 (MQ=255)
                                                                 tACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGt                                                                                 <  1:715938/150‑1 (MQ=255)
                                                                  aCCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTg                                                                                <  1:920760/150‑1 (MQ=255)
                                                                    ccAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGa                                                                                                  >  2:145349/1‑130 (MQ=255)
                                                                    ccAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGa                                                                                                  <  1:145349/130‑1 (MQ=255)
                                                                    ccAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGcc                                                                              >  2:381089/1‑149 (MQ=255)
                                                                    ccAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCACGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCa                                                                              >  1:962833/1‑150 (MQ=255)
                                                                        gTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATc                                                                          >  1:425644/1‑150 (MQ=255)
                                                                          cAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCaa                                                                        <  1:37536/150‑1 (MQ=255)
                                                                           aaCGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCaa                                                                                                                                 >  1:727647/1‑92 (MQ=255)
                                                                           aaCGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCaa                                                                                                                                 <  2:727647/92‑1 (MQ=255)
                                                                               cAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCATATCAACATTc                                                                   <  2:196095/150‑1 (MQ=255)
                                                                                   cc‑TNGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATag                                                               >  2:701646/1‑150 (MQ=255)
                                                                                      tGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATagag                                                             >  2:947817/1‑150 (MQ=255)
                                                                                              ttCNGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCtt                                                     >  2:697589/1‑150 (MQ=255)
                                                                                                      ggCCTTCGAGGGTGTCATGATCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGa                                             <  2:1097934/150‑1 (MQ=255)
                                                                                                         cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATaga                                                              >  1:1089517/1‑130 (MQ=255)
                                                                                                         cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATaga                                                              <  2:1089517/130‑1 (MQ=255)
                                                                                                          ttCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCt                                                <  2:138132/143‑1 (MQ=255)
                                                                                                          ttCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCt                                                >  1:138132/1‑143 (MQ=255)
                                                                                                                      aTGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGttatt                             >  1:140110/1‑150 (MQ=255)
                                                                                                                      aTGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGttatt                             <  2:867811/150‑1 (MQ=255)
                                                                                                                              tcAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTc                     >  1:365701/1‑150 (MQ=255)
                                                                                                                                aGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACAtt                                                                    >  2:2402/1‑101 (MQ=255)
                                                                                                                                aGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACAtt                                                                    <  1:2402/101‑1 (MQ=255)
                                                                                                                                      aGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCg             >  2:376379/1‑150 (MQ=255)
                                                                                                                                      aGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCg             >  1:390750/1‑150 (MQ=255)
                                                                                                                                           ttGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCtgt                                                                                       <  1:837084/71‑1 (MQ=255)
                                                                                                                                           ttGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCtgt                                                                                       >  2:837084/1‑71 (MQ=255)
                                                                                                                                            tGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTg       <  1:947817/150‑1 (MQ=255)
                                                                                                                                             gCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGTGATAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTgg      <  2:181424/150‑1 (MQ=255)
                                                                                                                                             gCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTgg      <  2:810013/150‑1 (MQ=255)
                                                                                                                                                 cccGGGCTCATGAGCACGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGc  >  1:1063008/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCCAGCACATGGATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCC‑TGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGC  >  NC_002947/5115603‑5115896

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