Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,197,222 G→A 10.3% A493V (GCG→GTG)  ethA ← FAD‑containing monooxygenase EthA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,197,2220GA10.3% 95.4 / 3.2 39A493V (GCG→GTG) ethAFAD‑containing monooxygenase EthA
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base A (2/2);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

GGCCGACGATGAGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAGGGTTCGCGCGGGCCCTGGCGCGGCATGCTGTC  >  NC_002947/3197077‑3197349
                                                                                                                                                 |                                                                                                                               
ggCCGACGATGAGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGAAAGTATAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCAtt                                                                                                                             >  1:352170/1‑150 (MQ=255)
  ccGACGATGAGCAGGCAAGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCCCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTgg                                                                                                                           <  1:31820/150‑1 (MQ=255)
       gaTGAGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACt                                                                                                                      >  1:389823/1‑150 (MQ=255)
        aTGAGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTg                                                                                                                     <  1:284752/150‑1 (MQ=255)
           aGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCgg                                                                                                                  <  2:6859/150‑1 (MQ=255)
            gCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGt                                                                                                                 <  1:286777/150‑1 (MQ=255)
            gCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGt                                                                                                                 <  1:485463/150‑1 (MQ=255)
                  cTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACtt                                                                                                           <  2:224545/150‑1 (MQ=255)
                     ccGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAg                                                                                                        <  2:389823/150‑1 (MQ=255)
                        gggCAGCTTGAGGGCTATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCa                                                                                                     <  2:352170/150‑1 (MQ=255)
                                  atggccatggcTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCACATTGGGTACTGCGGTGAACTTCAGGCACTtgtctgtc                                                                                           >  2:304243/1‑150 (MQ=255)
                                  atggccatggcTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCACATTGGGTACTGCGGTGAACTTCAGGCACTtgtctgtc                                                                                           >  2:304192/1‑150 (MQ=255)
                                     gccatggcTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAAATGGACATCGGGGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCAcc                                                                                        >  2:11628/1‑150 (MQ=255)
                                            cTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGTCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCAc                                                                                 <  2:452750/150‑1 (MQ=255)
                                                    gCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTACGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCgcagc                                                                         >  1:179652/1‑150 (MQ=255)
                                                      agagCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCTGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGcagcag                                                                       >  1:40592/1‑150 (MQ=255)
                                                            tGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGAc                                                                 >  1:254270/1‑150 (MQ=255)
                                                                 tCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTcgc                                                            >  1:241782/1‑150 (MQ=255)
                                                                         cAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAaataa                                                    >  1:100898/1‑150 (MQ=255)
                                                                          aGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGttt                                                                    >  1:589216/1‑133 (MQ=255)
                                                                          aGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGttt                                                                    <  2:589216/133‑1 (MQ=255)
                                                                                 aGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGa                                            <  2:368806/150‑1 (MQ=255)
                                                                                       gATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCACATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCAc                                      <  1:304243/150‑1 (MQ=255)
                                                                                       gATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCACATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCAc                                      <  1:304192/150‑1 (MQ=255)
                                                                                           ggTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAg                                  <  2:241782/150‑1 (MQ=255)
                                                                                            gTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCCGAAATAATCCCTCGACATCGCCCAgg                                 >  2:35519/1‑150 (MQ=255)
                                                                                                  aTCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACtt                                                                                                  <  1:496836/79‑1 (MQ=255)
                                                                                                  aTCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACtt                                                                                                  >  2:496836/1‑79 (MQ=255)
                                                                                                  aTCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTtgt                                                                                                >  1:395672/1‑81 (MQ=255)
                                                                                                  aTCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTtgt                                                                                                <  2:395672/81‑1 (MQ=255)
                                                                                                                    gctgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAg                                                                                                        >  1:415501/1‑55 (MQ=255)
                                                                                                                    gctgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAg                                                                                                        <  2:415501/55‑1 (MQ=255)
                                                                                                                     ctgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCgcag                                                                          >  2:401590/1‑84 (MQ=255)
                                                                                                                     ctgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCgcag                                                                          <  1:401590/84‑1 (MQ=255)
                                                                                                                      tgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAGGGTTCGCGCGGGCCCTGGCGCGGCATg       <  2:274333/150‑1 (MQ=255)
                                                                                                                      tgcAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGCCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAGGGTTCGCGCGGGCCCTGGCGCGGCATg       <  2:68712/150‑1 (MQ=255)
                                                                                                                           cTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATc                                              <  2:473134/106‑1 (MQ=255)
                                                                                                                           cTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATc                                              >  1:473134/1‑106 (MQ=255)
                                                                                                                           cTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAGGGTTCGCGCGGGCCCTGGCGCGGCATGCTGTc  <  2:40592/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                               
GGCCGACGATGAGCAGGACTGCCGGGGCGGCTTCATGGCCATGGCTGCGATAGCAGAGCCTGCGCTCGCCTGCCAGCTCGCAGAAGCGATCGGTGCTCATCGGCGGCTACCTTGCTGCTGCAACTGGACATCGGCGTGCAAGGGCGCATTGGGTACTGCGGTGAACTTCAGGCACTTGTCTGTCACCGCGCCACGGCGCAGCAGTTTGACGTCGCGGAAATAATCCATCGACATCACCCAGGGTTCGCGCGGGCCCTGGCGCGGCATGCTGTC  >  NC_002947/3197077‑3197349

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