Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,256 A→G 21.9% intergenic (‑520/+40) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,2560AG21.9% 32.1 / 7.4 23intergenic (‑520/+40)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base G (2/3);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

TGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGC  >  NC_002947/3223113‑3223392
                                                                                                                                               |                                                                                                                                        
tGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGGCGCCtt                                                                                                                                    >  2:332425/1‑150 (MQ=255)
   aaGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGGCGCCTTccc                                                                                                                                 <  1:332425/150‑2 (MQ=255)
        gCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGaa                                                                                                                            <  2:506206/150‑1 (MQ=255)
          aaTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGGCGCCTTccctggaa                                                                                                                            <  1:120916/148‑7 (MQ=255)
          aaTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGGCGCCTTccctggaa                                                                                                                            >  2:120916/1‑142 (MQ=255)
               cgcCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGGCGCCTTCCCTGGAAGgcgtca                                                                                                                     <  1:109644/150‑4 (MQ=255)
                  cAGAGCCTGCAAAAAGCAGGGCAGGGGCTATTGGCGGCTGGTGCTCGGGTTGGATCGGCAATGCGCGGGGAGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGt                                                                                                                  >  1:94724/1‑150 (MQ=255)
                         tGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGGAAGGCGGGGGGCGGCGGTTGCGCGACGACGGAAAGCGGCGAAGACACGCGCGGTTTGGATGGGGTGAGGCCTTCTAGGTAAGGCGCGGGGTCTACCgg                                                                                                           >  2:32669/1‑150 (MQ=255)
                                  cagggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                  <  1:97380/150‑1 (MQ=255)
                                     ggcaggCGCTATCGGGGCCTGGTGTTTGGGTGGGATCGGGAAGGCGCGGGGGGGGAGTTGCGCGACGACCGAAGGCGGCGAAGAGACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTGCCGGGGGGCGGCTTac                                                                                               <  2:94724/150‑2 (MQ=255)
                                          gcgcTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCGGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACtt                                                                                          <  1:559083/150‑1 (MQ=255)
                                                        tGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTg                                                                            <  1:305858/150‑1 (MQ=255)
                                                               tGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAAt                                                                     <  1:334083/150‑1 (MQ=255)
                                                               tGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAAt                                                                     <  2:381401/150‑1 (MQ=255)
                                                                    tGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGa                                                                <  2:420570/150‑1 (MQ=255)
                                                                      gATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACa                                                              >  1:108795/1‑150 (MQ=255)
                                                                        tCGGCAAGGCGCGGGGAGGCAGTTGCGCGACGACGGAAGCCGGCGAAGGCACCAGCCGCAGTCTGGGGATGACGCCTTCCAGGGAAGGCGCCGGGTCTACGGGGAGGCGACTTAGGACTTGGAGGTTTTGGGTGGAGAAATCATGACAAt                                                            >  2:21457/1‑150 (MQ=255)
                                                                           gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGaa                                                         <  2:108795/150‑1 (MQ=255)
                                                                                                                            gAGCGGTTTGGTGGGTGTGACGCCTTGAAGGGGAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGTCGGTTTTGGGTGGAGAAATCATGACAATGAAAGCGTAAATCAAGGCTTTCCCACGCTCTGCGCTAAAGGGTTCGGACTCg        >  2:211279/1‑150 (MQ=255)
                                                                                                                            gAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGATAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCg        >  1:482005/1‑150 (MQ=255)
                                                                                                                             aGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                  >  2:129970/1‑91 (MQ=255)
                                                                                                                             aGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                  <  1:129970/91‑1 (MQ=255)
                                                                                                                                 gTTTGGTGGGGGTGACGCGTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAgg   >  2:358148/1‑150 (MQ=255)
                                                                                                                                               aCGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgc  >  1:332971/1‑137 (MQ=255)
                                                                                                                                               aCGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgc  <  2:332971/137‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGC  >  NC_002947/3223113‑3223392

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