Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,052,8900AC20.0% 68.3 / 16.5 35L210R (CTG→CGG) ubiH2‑octaprenyl‑6‑methoxyphenol hydroxylase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base C (4/3);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGC  >  NC_000913/3052749‑3053038
                                                                                                                                             |                                                                                                                                                    
aTTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGc                                                                                                                                                <  6:38444/148‑1 (MQ=255)
      ccAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGtt                                                                                                                                          <  5:58776/148‑1 (MQ=255)
       ccaggcGCCAGCCCAAGGCCGTCTGGTGTTCCCGGCAAAACGTCCCGTCACGCCACGCCCACCCCCCTTCGCGCCGTTCCAGGGGCTGCCCCCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGc                                                                                                                                        <  5:102039/144‑1 (MQ=255)
          gTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAATGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGccga                                                                                                                                     <  3:13368/149‑4 (MQ=255)
                     aaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGc                                                                                                                                        >  7:34661/1‑135 (MQ=255)
                     aaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGc                                                                                                                                        >  5:14154/1‑135 (MQ=255)
                     aaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGc                                                                                                                                        <  6:14154/135‑1 (MQ=255)
                     aaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGc                                                                                                                                        <  8:34661/135‑1 (MQ=255)
                      aGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCaaaa                                                                                                                         <  5:116167/149‑1 (MQ=255)
                           gACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgcgcg                                                                                                                    <  7:9167/149‑1 (MQ=255)
                                gAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCCGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCtt                                                                                                               >  8:130266/1‑149 (MQ=255)
                                   ttCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCAt                                                                                                            <  7:130266/149‑1 (MQ=255)
                                                            cACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTgg                                                                                    >  7:14537/1‑148 (MQ=255)
                                                                     acCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACgg                                                                                                   <  6:9798/124‑1 (MQ=255)
                                                                     acCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGc                                                                          >  3:86937/1‑149 (MQ=255)
                                                                      ccTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGcc                                                                         <  8:14537/149‑1 (MQ=255)
                                                                       ctctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGtt                                                                                                                              >  7:57332/1‑95 (MQ=255)
                                                                       ctctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGtt                                                                                                                              <  8:57332/95‑1 (MQ=255)
                                                                       ctctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgcgc                                                                                                                     <  2:111662/104‑1 (MQ=255)
                                                                                   ttCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAg                                                            >  3:116898/1‑149 (MQ=255)
                                                                                      cAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCt                                                         >  3:47834/1‑149 (MQ=255)
                                                                                           gATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTc                                                                                                                             >  2:145798/1‑76 (MQ=255)
                                                                                           gATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTc                                                                                                                             <  1:145798/76‑1 (MQ=255)
                                                                                            aTGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTc                                                                                                              >  2:179800/1‑90 (MQ=255)
                                                                                            aTGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTc                                                                                                              <  1:179800/90‑1 (MQ=255)
                                                                                                cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCa                                               <  4:47834/149‑1 (MQ=255)
                                                                                                 accagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                              >  8:73283/1‑149 (MQ=255)
                                                                                                     gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCtgctgc                                                 <  2:121532/142‑1 (MQ=255)
                                                                                                     gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCtgctgc                                                 >  1:121532/1‑142 (MQ=255)
                                                                                                                   gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAg                                                                                                >  7:35263/1‑81 (MQ=255)
                                                                                                                   gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAg                                                                                                <  8:35263/81‑1 (MQ=255)
                                                                                                                    tccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTgg                           >  5:116219/1‑149 (MQ=255)
                                                                                                                           gctttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGc                    >  4:153231/8‑149 (MQ=255)
                                                                                                                            aCATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                   <  7:73283/149‑1 (MQ=255)
                                                                                                                                   cAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATggg            >  1:104601/1‑149 (MQ=255)
                                                                                                                                             aGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTGAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgc  >  4:101386/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
ATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGC  >  NC_000913/3052749‑3053038

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