Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,198,906 A→T 16.2% intergenic (‑250/+27) ethA ← / ← PP_2806 FAD‑containing monooxygenase EthA/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,198,9060AT16.2% 57.4 / 6.8 31intergenic (‑250/+27)ethA/PP_2806FAD‑containing monooxygenase EthA/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base T (3/2);  total (19/12)
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.23e-01

ACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAA  >  NC_002947/3198762‑3199040
                                                                                                                                                |                                                                                                                                       
aCTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACg                                                                                                                                       >  3:265223/1‑147 (MQ=255)
 cttcTATTTCAGGATTATTTGCACCTGTTACTTCGCCGATCAACGCTCAGACTGGGTACGTCTGTTAGATATCGTCATCGGAACGCTCGACGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgcag                                                                                                                                    <  6:233/146‑1 (MQ=255)
                 atTTGCACGTGTTACTTCGCCGTCCTTCCATCTGTCTGGGGCTGTCTGGGAGATCTCGGCCTCCGACCGTGCGACGGGGGGGGCTGCGGCGTTCCGCGCCGGTGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtgt                                                                                                                     <  8:45442/148‑1 (MQ=255)
                   ttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg                                                                                                                  >  7:141050/1‑149 (MQ=255)
                         cTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCt                                                                                                              >  3:63312/1‑147 (MQ=255)
                          tGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAg                                                                                                           >  7:18706/1‑149 (MQ=255)
                                  cGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGt                                                                                                          >  3:9456/1‑142 (MQ=255)
                                                 aGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑c                                                                                   >  6:151667/1‑149 (MQ=255)
                                                         tCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTgg                                                                            >  2:195813/1‑148 (MQ=255)
                                                                            cATCGGCACGTTCGACGGGGGGCGCGGCGGCCTGCCGCGCCGGCGGGGGCGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGc                                                        <  8:80501/149‑1 (MQ=255)
                                                                                    cGTGCGACGGGGGGCGCGGCGGCTTTCTGAGCCGGAGGGGGCGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGc                                                <  8:55079/149‑1 (MQ=255)
                                                                                     gTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCg                                               >  1:49807/1‑149 (MQ=255)
                                                                                      ttcgtcggggggggGTGCGGCTTTCCGAGCGGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGt                                              <  4:63312/143‑1 (MQ=255)
                                                                                         gAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGTGGTCTCGGTCTCGGCCTCCAGGCGTGCa                                           >  2:87733/1‑149 (MQ=255)
                                                                                          acggggggCGCGGCGGCGTTCCGAGCCGGAGGGGCAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCAc                                          <  8:141050/147‑1 (MQ=255)
                                                                                           aGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACg                                         >  6:82982/1‑149 (MQ=255)
                                                                                               gggCGCGGCGGCGTCCCCCGCCGGCCGGGGCGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACCTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTc                                      <  1:195813/149‑1 (MQ=255)
                                                                                                    ctccGGCTCTGAAAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTGAGTAACTGGCCACCAGCAACTTGA‑CTGGCTGGCAGGCGGTCACGGTCTCGGCGACCAGGCGTGCACGGTCGCCGt                                 >  6:42436/4‑148 (MQ=255)
                                                                                                                    cGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCATGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAg                 >  8:218144/1‑148 (MQ=255)
                                                                                                                          gggaggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAgg          <  2:49807/149‑1 (MQ=255)
                                                                                                                              ggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGt      <  5:82982/149‑1 (MQ=255)
                                                                                                                              ggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGt      >  7:63579/1‑149 (MQ=255)
                                                                                                                                 gaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCg    >  3:91818/1‑148 (MQ=255)
                                                                                                                                 attCGCCGTGCTGCGTCGCTGCACGGCG‑ATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGa                  >  2:226611/3‑133 (MQ=39)
                                                                                                                                 attCGCCGTGCTGCGTCGCTGCACGGCG‑ATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGa                  <  1:226611/131‑1 (MQ=39)
                                                                                                                                  aTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa  <  5:151667/149‑1 (MQ=255)
                                                                                                                                      ccGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCa                                           >  5:3512/1‑104 (MQ=255)
                                                                                                                                      ccGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCa                                           <  6:3512/104‑1 (MQ=255)
                                                                                                                                      ccGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa  >  1:134116/1‑145 (MQ=255)
                                                                                                                                      ccGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAgg                                                                        >  4:68798/1‑75 (MQ=255)
                                                                                                                                      ccGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa  >  2:225726/1‑145 (MQ=255)
                                                                                                                                                |                                                                                                                                       
ACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGA‑CTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAA  >  NC_002947/3198762‑3199040

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