Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,798,241 G→C 7.8% A156G (GCC→GGC)  qseB ← two‑component system DNA‑binding transcriptional activator QseB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,798,2410GC 7.8% 156.2 / 4.6 51A156G (GCC→GGC) qseBtwo‑component system DNA‑binding transcriptional activator QseB
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base C (2/2);  total (24/27)
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.20e-01

AATGCCGTTACCCAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGCAGGCGTGCG  >  NC_002947/4798093‑4798381
                                                                                                                                                    |                                                                                                                                            
aaTGGCGTTACCCAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGc                                                                                                                                             <  1:1175747/150‑1 (MQ=255)
   gCCGTTACCCAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTg                                                                                                                                          >  2:585439/1‑150 (MQ=255)
         actcAGCTTGCCACTCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACCCGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCGTGAGCAGGCCCTGCTccag                                                                                                                                    <  2:81062/147‑4 (MQ=255)
         aCCCAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcg                                                                                                                                    >  1:998112/1‑150 (MQ=255)
           ccAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggc                                                                                                                                  >  1:841879/1‑150 (MQ=255)
           ccAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggc                                                                                                                                  >  1:805716/1‑150 (MQ=255)
            cAGCTTGCGACGCAGGTGGTGGGTGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGCCAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggcg                                                                                                                                 <  2:619046/150‑1 (MQ=255)
             aGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggcgg                                                                                                                                <  2:998112/150‑1 (MQ=255)
             aGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggcgg                                                                                                                                <  1:150638/150‑1 (MQ=255)
             aGCTGGCGAGGCAGGTGGTGGATGTGCACAGTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcggcgg                                                                                                                                <  2:582551/150‑1 (MQ=255)
                  gCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAg                                                                                                                           <  1:1083085/150‑1 (MQ=255)
                         aGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACa                                                                                                                    <  1:506239/150‑1 (MQ=255)
                           gtggtgGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATg                                                                                                                  <  2:980267/150‑1 (MQ=255)
                            tggtggATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGAAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGc                                                                                                                 <  1:781181/150‑1 (MQ=255)
                                  aTGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTggccg                                                                                                           <  1:687059/150‑1 (MQ=255)
                                    gtgCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGc                                                                                                                 <  1:22083/142‑1 (MQ=255)
                                    gtgCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGc                                                                                                                 >  2:22083/1‑142 (MQ=255)
                                      gCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCa                                                                                                       <  1:905250/150‑1 (MQ=255)
                                        acaTTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAgg                                                                                                     <  1:773429/150‑1 (MQ=255)
                                           ttCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTc                                                                                                  <  2:805716/150‑1 (MQ=255)
                                              aGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGcc                                                                                               >  1:536843/1‑150 (MQ=255)
                                                  cGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAg                                                                                                      <  2:722700/139‑1 (MQ=255)
                                                  cGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAg                                                                                                      >  1:722700/1‑139 (MQ=255)
                                                            gACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCgg                                                                                 >  1:1179354/1‑150 (MQ=255)
                                                                    cgCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTTCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAgc                                                                         <  2:1260167/150‑1 (MQ=255)
                                                                     gCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAgcg                                                                        <  2:841879/150‑1 (MQ=255)
                                                                         aaaCCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCAGGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGa                                                                                                                      >  1:1361814/1‑100 (MQ=255)
                                                                         aaaCCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGa                                                                                                                      <  2:1361814/100‑1 (MQ=255)
                                                                         aaaCCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGCCCAGGGTCGCCTCGCAGCTCCTCGGGCCATCGCGCAGc                                                                    >  2:583263/1‑150 (MQ=255)
                                                                         aaaCCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTCGCCGGCCAGGATCGCCTCGCACCTAATCGGGCCATAGCGCAGc                                                                    >  2:583208/1‑150 (MQ=255)
                                                                          aaCCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCg                                                                   >  2:34059/1‑150 (MQ=255)
                                                                                  acGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcgg                                                                                                                                   >  1:1255267/1‑78 (MQ=255)
                                                                                  acGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCAcgg                                                                                                                                   <  2:1255267/78‑1 (MQ=255)
                                                                                   cGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCAGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGGCATGCTGGCCGCCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCACGGGGCCGTGCt                                                          >  1:209998/1‑150 (MQ=255)
                                                                                    gCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTc                                                         >  1:230892/1‑150 (MQ=255)
                                                                                      aGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGcc                                                                                               >  2:417998/1‑110 (MQ=255)
                                                                                      aGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGcc                                                                                               <  1:417998/110‑1 (MQ=255)
                                                                                        tCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCGTGAGCAGGCCCTGCTccag                                                                                                                                    >  1:893524/1‑68 (MQ=37)
                                                                                        tCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCGTGAGCAGGCCCTGCTccag                                                                                                                                    <  2:893524/71‑4 (MQ=37)
                                                                                                ctgctCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCAc                                             >  1:1168143/1‑150 (MQ=255)
                                                                                                           aCAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCgg                                  <  2:1383386/150‑1 (MQ=255)
                                                                                                               cACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAg                                                                                                                           >  1:834793/1‑57 (MQ=255)
                                                                                                               cACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAg                                                                                                                           <  2:834793/57‑1 (MQ=255)
                                                                                                                aCGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCAc                                                   <  2:147512/128‑1 (MQ=255)
                                                                                                                aCGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCAc                                                   >  1:147512/1‑128 (MQ=255)
                                                                                                                            ttCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAg                 >  2:1292588/1‑150 (MQ=255)
                                                                                                                                gcagcaGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTAt             <  1:34059/150‑1 (MQ=255)
                                                                                                                                  agcagcGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGc           >  2:76774/1‑150 (MQ=255)
                                                                                                                                      gcgcCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGCAGGc       <  1:1021920/150‑1 (MQ=255)
                                                                                                                                      gcgcCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGCAGGc       >  2:296050/1‑150 (MQ=255)
                                                                                                                                           gtgAGCAGGCCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGCAGGCGTGCg  >  2:730209/3‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AATGCCGTTACCCAGCTTGCGACGCAGGTGGTGGATGTGCACATTCAGGGCGTTGCTCTCGACTTCGTCGCTGAAACCATACACGCAGTCTTTCAGCTGCTCGCTGGACAGCACGCGCCCGGGGTTCTGCAGCAGCGCCTGGAGCAGGGCCTGCTCACGGCGGGACAGGTCGACATGCTGGCCGGCCAGGGTCGCCTCGCAGCTACTCGGGTCATAGCGCAGCGGGCCGTGCTCGATCACGTTCACCGCCCGCCCGGCCACCCGCCGCAGCAAGGTATGCAGGCGTGCG  >  NC_002947/4798093‑4798381

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