Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,456,504 G→A 7.4% H175H (CAC→CAT PP_3950 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,456,5040GA 7.4% 154.2 / 2.9 54H175H (CAC→CATPP_3950methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base A (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

GCTGTATGCCAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGA  >  NC_002947/4456356‑4456649
                                                                                                                                                      |                                                                                                                                                 
gCTGTATGCCAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGt                                                                                                                                                  <  2:1134760/150‑1 (MQ=255)
gCTGTATGCCAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGt                                                                                                                                                  <  2:211354/150‑1 (MQ=255)
       gCCAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGTCTCAGATGATGATGTCCCGGCCTCTCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAg                                                                                                                                           <  1:142139/150‑1 (MQ=255)
        ccAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAgg                                                                                                                                          >  2:1059159/1‑150 (MQ=255)
          aGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTAGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAgggg                                                                                                                                        >  2:42544/1‑150 (MQ=255)
                          tGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAAc                                                                                                                        <  1:1059159/150‑1 (MQ=255)
                                    aCAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAAc                                                                                                                        >  2:575919/1‑140 (MQ=255)
                                    aCAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAAc                                                                                                                        <  1:575919/140‑1 (MQ=255)
                                     cAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGAATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCCTCATGCACGAgggcaaagt                                                                                                                                   <  2:832000/129‑7 (MQ=255)
                                     cAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGAATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCCTCATGCACGAgggcaaagt                                                                                                                                   >  1:832000/1‑123 (MQ=255)
                                        gTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTg                                                                                                          >  1:909731/1‑150 (MQ=255)
                                           gTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGaa                                                                                                       <  2:325886/150‑1 (MQ=255)
                                           gTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGaa                                                                                                       <  1:981409/150‑1 (MQ=255)
                                                    cGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGa                                                                                              <  1:606392/150‑1 (MQ=255)
                                                        tCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAgg                                                                                                                                          >  2:719695/1‑102 (MQ=255)
                                                        tCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAgg                                                                                                                                          <  1:719695/102‑1 (MQ=255)
                                                        tCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGAc                                                                                          <  1:42544/150‑1 (MQ=255)
                                                            gggTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACggg                                                                                       >  2:1071055/1‑150 (MQ=255)
                                                              gTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTg                                                                                                                                     >  2:1075323/1‑102 (MQ=255)
                                                              gTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTg                                                                                                                                     <  1:1075323/102‑1 (MQ=255)
                                                              gTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGcc                                                                                     >  1:446189/1‑150 (MQ=255)
                                                                 aaTGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGAt                                                                                  <  1:375557/150‑1 (MQ=255)
                                                                       cccTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGt                                                                                                                      >  1:517802/1‑108 (MQ=255)
                                                                       cccTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGt                                                                                                                      <  2:517802/108‑1 (MQ=255)
                                                                            ccTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCa                                                                       >  2:26837/1‑150 (MQ=255)
                                                                              tCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGc                                                                     >  2:176853/1‑150 (MQ=255)
                                                                                  gcTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGt                                                                 <  1:1071055/150‑1 (MQ=255)
                                                                                    tGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCa                                                                       >  1:724736/1‑142 (MQ=255)
                                                                                    tGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCa                                                                       <  2:724736/142‑1 (MQ=255)
                                                                                       aTCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTg                                                            >  2:90164/1‑150 (MQ=255)
                                                                                        tCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCt                                                                                <  2:191706/129‑1 (MQ=255)
                                                                                        tCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCt                                                                                >  1:191706/1‑129 (MQ=255)
                                                                                         cAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTgca                                                          <  1:90164/150‑3 (MQ=255)
                                                                                          aGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAg                                                         <  2:446189/150‑1 (MQ=255)
                                                                                          aGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAg                                                         >  2:342901/1‑150 (MQ=255)
                                                                                                 ccAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCt                                                                                                                                      >  2:461415/1‑66 (MQ=255)
                                                                                                 ccAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCt                                                                                                                                      <  1:461415/66‑1 (MQ=255)
                                                                                                        aTCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGaaa                                           >  2:244059/1‑150 (MQ=255)
                                                                                                           gTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACgg                                        >  1:1198650/1‑150 (MQ=255)
                                                                                                           gTTGGGCACCCGCCCACACAGACAACCCAGCACTCTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACATCATCGTGGAAGCGCGAGAACGACGGTCCGATCTGCAAGGTCACCGAGCCGCTGTAAGGCCCAGACAGAAACgg                                        >  2:327639/1‑150 (MQ=255)
                                                                                                                gCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCtt                                   >  2:561394/1‑150 (MQ=255)
                                                                                                                        ccACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGcccc                           >  2:770874/1‑150 (MQ=255)
                                                                                                                              aGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGcc                     >  2:607189/1‑150 (MQ=255)
                                                                                                                               gACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCa                    <  1:176853/150‑1 (MQ=255)
                                                                                                                                 cAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCaga                  <  1:244059/150‑1 (MQ=255)
                                                                                                                                      cAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGc             >  2:789820/1‑150 (MQ=255)
                                                                                                                                      cAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGc             >  2:201034/1‑150 (MQ=255)
                                                                                                                                              agccccTCATGCACGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTc                                 <  2:1107870/120‑1 (MQ=255)
                                                                                                                                              agccccTCATGCACGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTc                                 >  1:1107870/4‑123 (MQ=255)
                                                                                                                                               gCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTc    <  1:607189/150‑1 (MQ=255)
                                                                                                                                                  ccTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGa  <  1:721408/150‑1 (MQ=255)
                                                                                                                                                  ccTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGa  >  2:569017/1‑150 (MQ=255)
                                                                                                                                                  ccTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGa  >  1:1003298/1‑150 (MQ=255)
                                                                                                                                                  ccTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGa  >  2:279194/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GCTGTATGCCAGGGTCGAAGCGCGATTGCGCCATGAACAGGTAGTTGTCGCCCGACTCGGGGTAAATGTGCCCTGCCTCGCGCTGGATCAGGTCACCCAGCACATCGTTGGGCACCCGCCCACACAGACAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGA  >  NC_002947/4456356‑4456649

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