Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,198,910 A→T 15.1% intergenic (‑254/+23) 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,9100AT15.1% 77.7 / 6.7 33intergenic (‑254/+23)ethA/PP_2806FAD‑containing monooxygenase EthA/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (20/8);  new base T (3/2);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

ACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGAT  >  NC_002947/3198762‑3199049
                                                                                                                                                    |                                                                                                                                           
aCTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgca                                                                                                                                             >  3:4812/1‑149 (MQ=255)
      aTTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACgg                                                                                                                                       >  7:52361/1‑149 (MQ=255)
                  tttGCGCCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACCCAGCACGGCGAATCCTGTgg                                                                                                                           >  7:302698/1‑149 (MQ=255)
                  tttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGATGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTgg                                                                                                                           >  7:134216/1‑149 (MQ=255)
                   ttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg                                                                                                                          >  5:35443/1‑149 (MQ=255)
                   ttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg                                                                                                                          >  3:76480/1‑149 (MQ=255)
                               cTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGCGGCGATGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCACCACGGCGAATCCTCTGGGGCCTCAGTAACt                                                                                                              >  7:317298/1‑149 (MQ=255)
                                              accAGACTGACCCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACtt                                                                                               >  7:222571/3‑149 (MQ=255)
                                                  gACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACt                                                                                           >  5:100465/1‑149 (MQ=255)
                                                  gACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACt                                                                                           >  7:39427/1‑149 (MQ=255)
                                                      gggTCAGTCTTTTATATTTCTTCATCGTAACTTTCGCTGGGGGGCGCTGCGGATTTCAGAGCCGGATGGGGTGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGcc                                                                                        <  8:134216/148‑1 (MQ=255)
                                                             tCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGc                                                                                >  4:121768/1‑149 (MQ=255)
                                                                 gTAGATATCGTCCTCGGAACGTTCGCCGGGGGGCGCTGCGGCTTTCCGAGCCGGAGGGGGCGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtc                                                                            <  6:168515/149‑1 (MQ=255)
                                                                 gTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtc                                                                            >  1:25633/1‑149 (MQ=255)
                                                                   aGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctc                                                                          >  2:60941/1‑149 (MQ=255)
                                                                             aTCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACt                                                                                                              >  3:263406/1‑103 (MQ=255)
                                                                                  aaCGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGCCTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAg                                                           >  8:309839/1‑149 (MQ=255)
                                                                                                          cTTTCAGAGCCGGAGGGGGAGGGGATCGCTGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACAAGCAACTTGACTGCCTGCCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGt                                          >  7:275441/1‑142 (MQ=255)
                                                                                                          cTTTAAGAGGCGGAGGGGGAGGGGATCGCTGTGAAGCGACGCAGCACGGAGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGt                                          <  8:275441/142‑1 (MQ=255)
                                                                                                           tttCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGAc                                  >  8:88624/1‑149 (MQ=255)
                                                                                                              cAGCGCCGGCGGGGGCGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCa                                <  8:61071/148‑1 (MQ=255)
                                                                                                                 agCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCGTACGTAAATGGCCACCAGCAAATTGAATCCATGGCAGCCGGTCACGGTCGCGGCCTCCAGCCGTCCACGGACGCCGTGACCGCCCCAGGCg                            >  8:67486/1‑149 (MQ=255)
                                                                                                                    cGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAg                               >  7:239322/1‑143 (MQ=255)
                                                                                                                    cGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAg                               <  8:239322/143‑1 (MQ=255)
                                                                                                                         ggggaggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAg                    <  8:52361/149‑1 (MQ=255)
                                                                                                                                 attCGCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGt                                                       <  1:58176/104‑1 (MQ=39)
                                                                                                                                 attCGCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGt                                                       >  2:58176/3‑106 (MQ=39)
                                                                                                                                    cGCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa           >  7:188418/1‑147 (MQ=255)
                                                                                                                                    cGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAAc          <  8:222571/148‑1 (MQ=255)
                                                                                                                                     gCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGc                                      >  8:136100/1‑119 (MQ=255)
                                                                                                                                     gCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGc                                      <  7:136100/119‑1 (MQ=255)
                                                                                                                                            aGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGAt  >  5:214740/1‑148 (MQ=255)
                                                                                                                                              cGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGa            >  8:242876/1‑136 (MQ=255)
                                                                                                                                              cGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGa            <  7:242876/136‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
ACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGAT  >  NC_002947/3198762‑3199049

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