Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,470,830 G→A 9.1% V118M (GTG→ATG)  PP_2165 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,470,8300GA 9.1% 115.5 / 3.2 44V118M (GTG→ATG) PP_2165two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base A (2/2);  total (20/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 = 8.93e-01

CGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGACCGCCGGTCGCCAGGTGCAGATGAACATG  >  NC_002947/2470683‑2470976
                                                                                                                                                     |                                                                                                                                                  
cGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTg                                                                                                                                                  >  1:137251/1‑150 (MQ=255)
     ggtcagCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTgcgccg                                                                                                                                             >  2:775485/1‑150 (MQ=255)
          gcgagcgTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGccc                                                                                                                                        >  2:330458/1‑150 (MQ=255)
            gagcgTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTc                                                                                                                                      >  1:225492/1‑150 (MQ=255)
              gcgTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGa                                                                                                                                    <  2:1043053/150‑1 (MQ=255)
                    gccgcAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGATGCGCCGCGCCCTCGACCGCTc                                                                                                                              >  2:689898/1‑150 (MQ=255)
                    gccgcAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGATGCGCCGCGCCCTCGACCGCTc                                                                                                                              >  2:689813/1‑150 (MQ=255)
                           ttGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCtt                                                                                                                       >  1:644422/1‑150 (MQ=255)
                                  tGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGGGGGAGCGGTGCGGCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGGGCTgg                                                                                                                >  1:892281/1‑150 (MQ=255)
                                                ggTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTAc                                                                                                  >  2:677148/1‑150 (MQ=255)
                                                   gCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACcgc                                                                                               <  2:532385/150‑1 (MQ=255)
                                                      ggCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGAc                                                                                            <  2:855451/150‑1 (MQ=255)
                                                         gTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCCCGCCTTGTGCTGGAGAGCCTCCGCTAGCGCGACAAg                                                                                         >  2:238617/1‑150 (MQ=255)
                                                          tCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGc                                                                                        <  2:769526/150‑1 (MQ=255)
                                                            aTGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTc                                                                                      >  1:1126305/1‑150 (MQ=255)
                                                                 tGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAgc                                                                                 <  2:372452/150‑1 (MQ=255)
                                                                   aTGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCAGTGCAGACGTGTGCTCTTCCGATACATTCGGTGCGCCTCGTCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAgccg                                                                               <  1:524002/150‑1 (MQ=255)
                                                                           gAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACcgc                                                                       <  1:1108944/150‑1 (MQ=255)
                                                                                 ttGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTg                                                                 >  1:986778/1‑150 (MQ=255)
                                                                                     gcCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAgg                                                             <  1:330458/150‑1 (MQ=255)
                                                                                           gcgcCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGcc                                                       <  2:184085/150‑1 (MQ=255)
                                                                                            cgccgcCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCt                                                      <  1:238617/150‑1 (MQ=255)
                                                                                               cgccgACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCa                                                   <  2:687792/150‑1 (MQ=255)
                                                                                                gccgACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGATGCGCCGCGCCCTCGCCCGCCTGCACATTGTGCTGCAACGCCTTCTCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCAc                                                  <  1:689898/150‑1 (MQ=255)
                                                                                                gccgACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGATGCGCCGCGCCCTCGACAAATACCCGATAGTGCTGGAACGCCGGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCAc                                                  <  1:689813/150‑1 (MQ=255)
                                                                                                       aCCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGc                                           <  2:117725/150‑1 (MQ=255)
                                                                                                         cTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCag                                         >  1:649362/1‑150 (MQ=255)
                                                                                                          tcatcaAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAg                                                              >  1:960817/1‑128 (MQ=255)
                                                                                                          tcatcaAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAg                                                              <  2:960817/128‑1 (MQ=255)
                                                                                                            atcaAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCaggag                                      <  2:137251/150‑1 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAAcc                                                                         <  2:715966/110‑1 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAAcc                                                                         >  1:715966/1‑110 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCaggag                                      >  2:586645/1‑145 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCaggag                                      <  1:586645/145‑1 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCa                                 <  2:225492/150‑1 (MQ=255)
                                                                                                                 gCCGCTGGAAGACCTGGCCGTACCCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCa                                 <  1:50544/150‑1 (MQ=255)
                                                                                                                    gCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGAc                              <  2:986778/150‑1 (MQ=255)
                                                                                                                       ggAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCaggag                                      >  2:1113696/1‑139 (MQ=255)
                                                                                                                       ggAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCaggag                                      <  1:1113696/139‑1 (MQ=255)
                                                                                                                            aCCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGACCGCCGGTc                      >  1:343532/1‑150 (MQ=255)
                                                                                                                              cTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGACCGCCGGTCGc                    >  2:348526/1‑150 (MQ=255)
                                                                                                                                                 ttCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGa                               >  2:3150/1‑122 (MQ=255)
                                                                                                                                                 ttCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGa                               <  1:3150/122‑1 (MQ=255)
                                                                                                                                                  tCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGACCGCCGGTCGCCAGGTGCAGATGAACATg  <  1:775485/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCAAGCCGCTGGAAGACCTGGCCGTACTCGAGCATTCGGTGCGCCGCGCCCTCGACCGCTCGCGCCTTGTGCTGGAAAACCAGCGCTACCGCGACAAGCTCGAAGCCGCCAACCGCGAACTGGAGGCCAGCCTGCACCTGTTGCAGGAGGACCAGACCGCCGGTCGCCAGGTGCAGATGAACATG  >  NC_002947/2470683‑2470976

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: