Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,960,142 T→A 13.4% intergenic (+50/‑78) bolA → / → PP_1758 DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,960,1420TA13.4% 141.8 / 12.2 52intergenic (+50/‑78)bolA/PP_1758DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base A (4/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGC  >  NC_002947/1959996‑1960288
                                                                                                                                                  |                                                                                                                                                  
gaTCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                                 <  2:128019/150‑1 (MQ=255)
 aTCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAACATTTTTTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCg                                                                                                                                                <  1:521140/150‑1 (MQ=255)
    cATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAA‑GGGGCCGCTAAGCGGccc                                                                                                                                            >  1:1116242/1‑150 (MQ=255)
          cTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGc                                                                                                                                       >  1:537914/1‑150 (MQ=255)
           tGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAA‑GGGGCCGCTAAGCGGccccctttgg                                                                                                                                     <  2:1116242/150‑7 (MQ=255)
               cATCCACACCTACACCGCCGAAGAGTGCGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCttttt                                                                                                                                  >  1:635221/1‑150 (MQ=255)
                    acacCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAGCGGccccctttg                                                                                                                                      <  2:107043/141‑6 (MQ=255)
                    acacCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAGCGGccccctttg                                                                                                                                      >  1:107043/1‑136 (MQ=255)
                        cTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACg                                                                                                                         >  1:1084171/1‑150 (MQ=255)
                          acacCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTc                                                                                                                       <  1:325546/150‑1 (MQ=255)
                             ccgccgAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCaaa                                                                                                                    >  1:529593/1‑150 (MQ=255)
                                 cgAAGAGTGGGCCAAGGTCGGCGTTCCACCGGCCTCGCCCGTGTGTGTGGGCGGCGGGCACTGAAACCGTTCTGTTGTTCTGGTAGAATCCGCATCACCCAAAGGAGGCCGCTTAGCGGCCCCTTGCTCTTTATTCAACCGCCAAACCCg                                                                                                                >  1:431330/1‑150 (MQ=255)
                                            ccAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGTTCCGGTTACCGTCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCtt                                                                                                     <  2:347617/150‑1 (MQ=255)
                                                 gTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGc                                                                                                                              <  2:779346/120‑1 (MQ=255)
                                                 gTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGc                                                                                                                              >  1:779346/1‑120 (MQ=255)
                                                  tCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAgg                                                                                               <  2:537914/150‑1 (MQ=255)
                                                     gCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTAc                                                                                                   <  1:1060367/143‑1 (MQ=255)
                                                     gCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTAc                                                                                                   >  2:1060367/1‑143 (MQ=255)
                                                            aCCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCt                                                                                     >  2:588245/1‑150 (MQ=255)
                                                                       ccGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                          <  1:588245/150‑1 (MQ=255)
                                                                         gtgtgtGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTCGCACCAGCACCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTa                                                                        <  2:635221/150‑1 (MQ=255)
                                                                              tgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                          <  2:1089585/143‑1 (MQ=255)
                                                                              tgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                          >  1:1089585/1‑143 (MQ=255)
                                                                              tgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGt                                                                   >  1:351293/1‑150 (MQ=255)
                                                                                     ggcggGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAg                                                                                                >  2:1164238/1‑114 (MQ=255)
                                                                                     ggcggGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAg                                                                                                <  1:1164238/114‑1 (MQ=255)
                                                                                        gggCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATNACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCAt                                                         >  1:766610/1‑150 (MQ=255)
                                                                                              tGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  >  2:956588/1‑119 (MQ=255)
                                                                                              tGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  <  1:956588/119‑1 (MQ=255)
                                                                                                   ccGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  <  1:104412/114‑1 (MQ=255)
                                                                                                   ccGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  >  2:104412/1‑114 (MQ=255)
                                                                                                              ttCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAg                                   >  2:949970/1‑150 (MQ=255)
                                                                                                               tCTGGTAGACTGCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGACAACGCCCAGCTCTGTCCGCACCTTACGAGGGCAACCACCTGGAGATTCCGTTACATTTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGt                                  <  1:390662/150‑1 (MQ=255)
                                                                                                               tCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATCTATTCGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGt                                  <  2:1069008/150‑1 (MQ=255)
                                                                                                                 tGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTc                                >  2:111952/1‑150 (MQ=255)
                                                                                                                   gTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                          <  2:34028/106‑1 (MQ=255)
                                                                                                                   gTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                          >  1:34028/1‑106 (MQ=255)
                                                                                                                   gTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGt                              >  1:903000/1‑150 (MQ=255)
                                                                                                                           cGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTgg                      <  1:111952/150‑1 (MQ=255)
                                                                                                                             cATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                    <  2:903000/150‑1 (MQ=255)
                                                                                                                             cATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                    <  2:916523/150‑1 (MQ=255)
                                                                                                                              aTCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  <  2:1162501/87‑1 (MQ=255)
                                                                                                                              aTCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGga                                                                                  >  1:1162501/1‑87 (MQ=255)
                                                                                                                                 aCCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAgccc                >  1:117981/1‑147 (MQ=255)
                                                                                                                                      aaGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTc           <  1:949970/150‑1 (MQ=255)
                                                                                                                                       aaggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCg          >  2:281667/3‑150 (MQ=255)
                                                                                                                                          gggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGtc                                                                  >  2:143336/1‑91 (MQ=255)
                                                                                                                                          gggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGtc                                                                  <  1:143336/91‑1 (MQ=255)
                                                                                                                                           ggCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCt      <  2:1084171/150‑1 (MQ=255)
                                                                                                                                            gCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTg     >  1:636964/1‑150 (MQ=255)
                                                                                                                                            gCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTg     >  1:172600/1‑150 (MQ=255)
                                                                                                                                               gCTTAGCGGCCCCTTGGTTTTTATGCAAACGGCAAACCCAGTCCACCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCAGCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTgcgc  <  2:83764/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGC  >  NC_002947/1959996‑1960288

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