Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,561,787 G→A 7.8% P116S (CCC→TCC)  pctA ← methyl‑accepting chemotaxis protein PctA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,561,7870GA 7.8% 143.7 / 2.6 51P116S (CCC→TCC) pctAmethyl‑accepting chemotaxis protein PctA
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base A (2/2);  total (26/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 = 8.87e-01

ACGCCAACGCTGCTGGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTCCACCAGTT  >  NC_002947/1561642‑1561936
                                                                                                                                                 |                                                                                                                                                     
aCGCCAACGCTGCTGGAGACCTTGCTGAGGATGCCGAATCAGATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATc                                                                                                                                                   <  1:12189/150‑1 (MQ=255)
 cGCCAACGCTGCTGGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCt                                                                                                                                                  >  2:1106165/1‑150 (MQ=255)
            ctgGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCggg                                                                                                                                       >  2:1059103/1‑150 (MQ=255)
             tgGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTcgc                                                                                                                                              <  1:1212258/142‑1 (MQ=255)
             tgGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTcgc                                                                                                                                              >  2:1212258/1‑142 (MQ=255)
               gagaCCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTc                                                                                                                                                >  1:847126/1‑138 (MQ=255)
               gagaCCTCGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTc                                                                                                                                                <  2:847126/138‑1 (MQ=255)
                    cTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCtt                                                                                                                                                 >  2:82568/1‑132 (MQ=255)
                    cTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCtt                                                                                                                                                 <  1:82568/132‑1 (MQ=255)
                        cTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGaaa                                                                                                                           <  2:700186/150‑1 (MQ=255)
                             gATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAAccgcc                                                                                                                      >  2:1043788/1‑150 (MQ=255)
                               ggCCGATGACCATCTTGTTGGTGGACATGTCGATGTAGGGTTCCGCCAGCGTCGCGCCGCTGGCGTTCATGCCCTCCTTGTACCAAGGGCGGATGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGt                                                                                                                    <  2:153840/149‑1 (MQ=255)
                                gCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTc                                                                                                                   >  1:897750/1‑150 (MQ=255)
                                    ttGACCATCTTGGTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACt                                                                                                               <  2:1016929/149‑1 (MQ=255)
                                        ccATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGc                                                                                                           <  2:258255/150‑1 (MQ=255)
                                                          tGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGa                                                                                         <  2:897750/150‑1 (MQ=255)
                                                               aTGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTg                                                                                    <  2:68874/150‑1 (MQ=255)
                                                                tGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTgg                                                                                   >  1:295166/1‑150 (MQ=255)
                                                                tGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCAAGGTAGACCGTGAGGAAGCTgt                                                                                   <  2:1244946/150‑2 (MQ=255)
                                                                    gggTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGatat                                                                               <  1:1106165/150‑1 (MQ=255)
                                                                       ttCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAg                                                                            >  1:386226/1‑150 (MQ=255)
                                                                        tCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAgc                                                                           <  2:665721/150‑1 (MQ=255)
                                                                                          cTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAgg                                                         <  1:846170/150‑1 (MQ=255)
                                                                                               gTTCATGCCTTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTc                                                    >  2:443412/1‑150 (MQ=255)
                                                                                               gTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTAGTAGCCGTCGGGCATCTTCGCGTCCGGGCACACGCTGAAACCGCCCCCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGCCCCATCAGGGTTTc                                                    >  1:1054501/1‑150 (MQ=255)
                                                                                                ttCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCg                                                   >  1:330869/1‑150 (MQ=255)
                                                                                                   aTGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTAGAGGCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCg                                                >  2:113982/1‑150 (MQ=255)
                                                                                                      ccGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTg                                             >  2:69364/1‑150 (MQ=255)
                                                                                                       cGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTgg                                            >  1:94907/1‑150 (MQ=255)
                                                                                                       cGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGACATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAgc                                                                           >  1:98394/1‑119 (MQ=255)
                                                                                                       cGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGACATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAgc                                                                           <  2:98394/119‑1 (MQ=255)
                                                                                                       cGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGACATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAgc                                                                           <  2:98320/119‑1 (MQ=255)
                                                                                                       cGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGACATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAgc                                                                           >  1:98320/1‑119 (MQ=255)
                                                                                                        gTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCt                                                                                     <  1:468750/108‑1 (MQ=255)
                                                                                                        gTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCt                                                                                     >  2:468750/1‑108 (MQ=255)
                                                                                                        gTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCAt                                                             <  2:1278382/132‑1 (MQ=255)
                                                                                                        gTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCAt                                                             >  1:1278382/1‑132 (MQ=255)
                                                                                                               tACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGc                                                 <  2:334897/137‑1 (MQ=255)
                                                                                                               tACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGc                                                 >  1:334897/1‑137 (MQ=255)
                                                                                                                aCCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCt                                   <  2:94907/150‑1 (MQ=255)
                                                                                                                 ccAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGa                                                                  <  2:27/118‑1 (MQ=255)
                                                                                                                 ccAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGa                                                                  >  1:27/1‑118 (MQ=255)
                                                                                                                   aaGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTggg                                <  2:295166/150‑1 (MQ=255)
                                                                                                                            tGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCt                       >  2:811961/1‑150 (MQ=255)
                                                                                                                               gcgGGTCGTAGCCGTCTGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTgcg                    >  2:137070/1‑150 (MQ=255)
                                                                                                                                      gTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGcc                  <  1:1040951/145‑1 (MQ=255)
                                                                                                                                      gTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGcc                  >  2:1040951/1‑145 (MQ=255)
                                                                                                                                         gCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTc          <  1:69364/150‑1 (MQ=255)
                                                                                                                                            gTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTccac       >  1:894748/1‑150 (MQ=255)
                                                                                                                                            gTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTccac       >  2:542597/1‑150 (MQ=255)
                                                                                                                                              cGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTccacca     <  1:443412/150‑1 (MQ=255)
                                                                                                                                                 gCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTCCACCAGtt  >  2:1223816/1‑150 (MQ=255)
                                                                                                                                                 gCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTCCACCAGtc  <  1:1059103/150‑2 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ACGCCAACGCTGCTGGAGACCTTGCTGAGGATGCCGATGACCATCTTGTTGGTGGTCATGTCGATGTAGGGTTCGGTCAGCGTCGCGCCGCTGGCGTTCATGCCGTCCTTGTACCAAGGGCGGGTGCGCGGGTCGTAGCCGTCGGGCATCTTCGCGTCCGGGCGCACGCTGAAACCGCCGTCGACTTTGCCCAGGTAGACCGTGAGGAAGCTGGATATCAGCGCATTCTGACCCATCAGGGTTTCGGCGTTGGCCGGGTCCTGGGCAATGTTCTGCGCCAGGTTTTCCACCAGTT  >  NC_002947/1561642‑1561936

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