Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,604,263 C→A 7.9% G145V (GGC→GTC)  PP_3179 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,604,2630CA 7.9% 146.3 / 3.5 51G145V (GGC→GTC) PP_3179LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (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 = 5.69e-01

ACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGCCGTGATG  >  NC_002947/3604116‑3604412
                                                                                                                                                   |                                                                                                                                                     
acGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGAACGGCAGCGGTAGCCTGCGCATGAAGCCg                                                                                                                                                     >  1:681904/1‑150 (MQ=255)
  gAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGAt                                                                                                                                                   >  2:6525/1‑150 (MQ=255)
  gAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATAAAGCCGAt                                                                                                                                                   <  2:1037569/150‑1 (MQ=255)
      tCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACTATCTGCTCGTCAAGCAGCAAGCGAGCCTCCATGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTcc                                                                                                                                               <  1:725419/150‑1 (MQ=255)
      tCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGGAGCCTGCGCATGAAGCCGATGtac                                                                                                                                               >  2:106126/1‑148 (MQ=255)
            tCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGGGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCg                                                                                                                                         >  1:11739/1‑150 (MQ=255)
                 aaaCTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGc                                                                                                                                    <  1:6525/150‑1 (MQ=255)
                   aCTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGAt                                                                                                                                                   >  2:1035166/1‑133 (MQ=255)
                   aCTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGAt                                                                                                                                                   <  1:1035166/133‑1 (MQ=255)
                              tCGACGGCCTCGAAACGCGCCATGGAATGGCTGGCGGGCAGCGCCATGACCCTCTGTTAGTTAAGCAGCAAGCGCGCCTCCAGGGTGCCGAGCGGCAGCGGTAGCCTGCGGATGAAGCGGATCTCCGCCTTGCCGGCAAGAAAGGCCTTg                                                                                                                       >  2:540484/1‑150 (MQ=255)
                                gatctccTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGCCGATCTGCTCGTCAAGCAGCAAGCGAGCATTCAGGGTGCCGACCGGCAGCAGCTGCCTGCGCAGGATGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGc                                                                                                                     <  1:1041007/145‑1 (MQ=255)
                                      cTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCaa                                                                                                                              <  1:371395/135‑1 (MQ=255)
                                      cTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCaa                                                                                                                              >  2:371395/1‑135 (MQ=255)
                                      cTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGccagatc                                                                                                               >  2:1041007/1‑145 (MQ=255)
                                               cgcCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATgg                                                                                                      <  2:709796/150‑1 (MQ=255)
                                                gcCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     >  2:165654/1‑150 (MQ=255)
                                                          ggctggcGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGa                                                                                           <  1:999956/150‑1 (MQ=255)
                                                          ggctggcGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGa                                                                                           <  1:999805/150‑1 (MQ=255)
                                                            ctggcGGGCAGCGCCATGCCGATCTGCTCGTCAATCAGCAAGCCAGCCTCCAGGCGGCCGACCGGCAGCGGGAGCCTGGGCATGTAGCCGATGTCCGCGTTCCCGGCCAGCAAGACATTGGTCTGAGAGGCAGATGGCACCTCACGCAtc                                                                                         <  2:205725/150‑2 (MQ=255)
                                                            ctggcGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATg                                                                                         <  1:659490/150‑1 (MQ=255)
                                                                           aTGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGttt                                                                          <  1:614621/150‑1 (MQ=255)
                                                                                gATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATa                                                                     >  2:277561/1‑150 (MQ=255)
                                                                                gATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATa                                                                     >  1:615872/1‑150 (MQ=255)
                                                                                          tCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGa                                                           >  1:332249/1‑150 (MQ=255)
                                                                                            aaGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGaaa                                                         >  1:533175/1‑150 (MQ=255)
                                                                                               cagcaAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAgcg                                                      <  1:165654/150‑1 (MQ=255)
                                                                                                agcaAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAgcgg                                                     <  1:277561/150‑1 (MQ=255)
                                                                                                       gAGNCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGc                                              >  2:699913/1‑150 (MQ=255)
                                                                                                       gAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGc                                              <  2:533175/150‑1 (MQ=255)
                                                                                                          ccTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCggg                                           >  1:189775/1‑149 (MQ=255)
                                                                                                          ccTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCCCAAGCGGt                                           >  1:1047462/1‑150 (MQ=255)
                                                                                                           cTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGtt                                          >  1:152100/1‑150 (MQ=255)
                                                                                                                gggTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGagca                                     <  2:332249/150‑1 (MQ=255)
                                                                                                                gggTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGagca                                     >  1:532418/1‑150 (MQ=255)
                                                                                                                       gACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATc                              >  1:954003/1‑150 (MQ=255)
                                                                                                                       gACCGCCAGCGGTAGCCTGCGCATGAAGCCCATGTCCGCCTCGCCGCCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATc                              <  1:106126/150‑1 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     <  1:802138/77‑1 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     >  2:802138/1‑77 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGACGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     <  1:802784/77‑1 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGACGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     >  2:802784/1‑77 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGACGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     >  2:242919/1‑77 (MQ=255)
                                                                                                                         ccGGCAGCGGTAGCCTGCGCATGAAGACGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGc                                                                                                     <  1:242919/77‑1 (MQ=255)
                                                                                                                          cGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGta                           <  1:646268/150‑1 (MQ=255)
                                                                                                                           ggCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGtaa                          <  2:954003/150‑1 (MQ=255)
                                                                                                                                      ccTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTTAGAGGTACGCTTGGTTAATCCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGa               >  1:407554/1‑150 (MQ=255)
                                                                                                                                           gcATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATgg          >  2:628196/1‑150 (MQ=255)
                                                                                                                                            cATGAGGCGGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGAATCTCGCGGGTGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGc         <  2:504543/150‑1 (MQ=255)
                                                                                                                                            cATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGAt                                                                      <  2:290651/89‑1 (MQ=255)
                                                                                                                                            cATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGAt                                                                      >  1:290651/1‑89 (MQ=255)
                                                                                                                                                 aGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGCCGTGa    <  1:540484/150‑1 (MQ=255)
                                                                                                                                                  gCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGCCGTGAt   >  1:358416/1‑150 (MQ=255)
                                                                                                                                                   ccGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGCCGTGATg  <  2:681904/150‑1 (MQ=255)
                                                                                                                                                   ccGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGATGTACGTTTGGCTCCGTCTCGCGCAGGAGGCGCGGGCGGTTGTGCAACGAACCCTAATAAACAATTCAAATGGCCGTGATg  <  2:11739/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAAACGCGCCATGGGATGGCTGGCGGGCAGCGCCATGACGATCTGCTCGTCAAGCAGCAAGCGAGCCTCCAGGGTGCCGACCGGCAGCGGTAGCCTGCGCATGAAGCCGATGTCCGCCTCGCCGGCCAGCAAGGCCTTGGCCTGAGAGGCAGATGGCATCTCGCGGATGGTCAGCTGTACGTTTGGATAAGTCTCCCGAAAGCGGCGCAAGCGGTTGAGCAGCGAATCGTAATAAGCGATGGAAATGGCCGTGATG  >  NC_002947/3604116‑3604412

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