Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,587,557 G→T 20.6% intergenic (‑12/‑246) ttgA ← / → ttgR efflux pump periplasmic linker TtgA/HTH‑type transcriptional regulator TtgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,587,5570GT20.6% 55.4 / 3.5 24intergenic (‑12/‑246)ttgA/ttgRefflux pump periplasmic linker TtgA/HTH‑type transcriptional regulator TtgR
Reads supporting (aligned to +/- strand):  ref base G (5/14);  new base T (2/3);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GTCACGACGCCGACCTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGGGTGGTACT  >  NC_002947/1587424‑1587692
                                                                                                                                     |                                                                                                                                       
gTCACGACGCCGACCTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGt                                                                                                                                       <  2:517060/136‑1 (MQ=255)
gTCACGACGCCGACCTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGt                                                                                                                                       <  2:406825/136‑1 (MQ=255)
  cACGACGCCGACCTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCTGTAACGGCTTGATTGAATTGCATGAGGATCCTCGGGTCg                                                                                                                                     <  2:528557/136‑1 (MQ=255)
             ccTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAggg                                                                                                                          <  1:509078/136‑1 (MQ=255)
             ccTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAggg                                                                                                                          <  1:521359/136‑1 (MQ=255)
                  ggagCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACg                                                                                                                     <  2:265660/136‑1 (MQ=255)
                    agCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCt                                                                                                                   >  1:232076/1‑136 (MQ=255)
                    agCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCt                                                                                                                   >  1:279887/1‑136 (MQ=255)
                       cTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCTTGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCaa                                                                                                                <  2:378045/136‑1 (MQ=255)
                         gcgcCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCCGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGGAACGGCTGGCTTGAATTGCATGAGGATCCTCGTGTCGCTGGAGCATGGAGACGCTCAaga                                                                                                              >  2:454397/1‑136 (MQ=255)
                         gcgcCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGTGTCGCTGGAGCAGGGAGACGCTCAaga                                                                                                              >  2:255576/1‑136 (MQ=255)
                         gcgcCGCTGGCGCTGCTTCTTCTTTCATACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAaga                                                                                                              >  2:530767/1‑136 (MQ=255)
                            ccgcTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTa                                                                                                           >  2:470921/1‑136 (MQ=255)
                                tGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATgg                                                                                                       <  1:470921/136‑1 (MQ=255)
                                tGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATgg                                                                                                       <  1:269368/136‑1 (MQ=255)
                                                        gCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAAt                                                                               >  2:345137/1‑136 (MQ=255)
                                                                                       gAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGa                                                <  1:136139/136‑1 (MQ=255)
                                                                                                       cGGCTGGCTTGAATTGCATGAGGATCCTCGTGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCAt                                <  2:526802/136‑1 (MQ=255)
                                                                                                                           aGGATCCTCGTGTCGCTGGAGCATGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCTTGGTTGTTTGTAAATActgct            <  1:454397/136‑1 (MQ=255)
                                                                                                                           aGGATCCTCGTGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATActgct            <  1:255576/136‑1 (MQ=255)
                                                                                                                           aGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATActgct            <  1:530767/136‑1 (MQ=255)
                                                                                                                             gATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTgg          <  2:279887/136‑1 (MQ=255)
                                                                                                                                  tCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGggtggt     <  1:489199/136‑1 (MQ=255)
                                                                                                                                  tCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGggtggt     <  1:328607/136‑1 (MQ=255)
                                                                                                                                     ggTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGGGTGGTACt  >  2:130933/1‑136 (MQ=255)
                                                                                                                                     ggTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGGGTGGTACt  >  2:139468/1‑136 (MQ=255)
                                                                                                                                     ggTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGGGTGGTACt  >  2:254275/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GTCACGACGCCGACCTGAGGAGCCTGCGCCGCTGGCGCTGCTTCTTCTTTCTTACAGCCACTGAGCAGGGTTGCCAGGGCGACGGCGGAAACCAGAGCGGTAACGGCTGGCTTGAATTGCATGAGGATCCTCGGGTCGCTGGAGCAGGGAGACGCTCAAGAGTAATGGAAGAGATCTTGATCTGGAAAAAATGCTATCCGGTGGATAAATAGCTTGCTAAGGAATATACTTACATTCATGGTTGTTTGTAAATACTGCTGGGTGGTACT  >  NC_002947/1587424‑1587692

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