Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091367,8872.T38.1% 23.2 / 69.1 42coding (451/1503 nt)araAL‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base . (12/14);  new base T (6/10);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTG  >  NC_000913/67748‑68028
                                                                                                                                             |                                                                                                                                               
gCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCCTGGGCTTGCTTGGCGTGCCA‑‑GT‑‑g                                                                                                                                            >  2:2284/1‑143 (MQ=255)
        aTCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑gtgaccggt                                                                                                                                    <  2:302411/143‑9 (MQ=255)
        aTCGCCATCCGTGGCCGCCGCGTCACGCATGGTATCGCCAAATCGGCAGACTTTCAGATGGCGGGTATCCTGTTTAGAGGCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑gtgaccggt                                                                                                                                    <  2:545084/143‑9 (MQ=255)
           gCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAAc                                                                                                                                   <  1:573598/143‑1 (MQ=255)
              aTCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCAc                                                                                                                                <  1:659212/143‑1 (MQ=255)
               tCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCaagaa                                                                                                                             >  2:688029/1‑139 (MQ=255)
                   tGACCGCTACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGt                                                                                                                         <  2:250326/143‑1 (MQ=255)
                    gACCGCCACGTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGtt                                                                                                                        <  1:629122/143‑1 (MQ=255)
                            cgtcccGCCTGTTATCGCCAAATCGGCAGACTTTCAGATGACGTGTA‑CCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACg                                                                                                                 <  2:411847/138‑1 (MQ=255)
                               catggaTGTCATCGCTCAATCGGAAGACTTTGAGATGACGGGTACCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCa                                                                                                               <  1:121614/138‑1 (MQ=255)
                                       ttATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCCTGTTGCCGACGCATAcgcgcg                                                                                                       >  2:153660/1‑143 (MQ=255)
                                           cGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCACTGT‑‑GAAAGACAAACATGGCATGTTGCTGACGCATACGCGCGCCaa                                                                                                   <  1:593580/143‑1 (MQ=255)
                                             ccAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCCCCTGGGCTTGTTTATCCTGCCC‑‑GTGGGACCGGGAACCCCGGCCAGTTGCTGGCGTATACGTGCGCCGATg                                                                                                 >  1:512486/1‑143 (MQ=255)
                                               gaaggggtcggcTTTCGGATGCAGGGTATCCTGTTTCGAGAGTGCCTGACGCATGTGGGGGGCGTGGCTGGTATGGGGTTGTCGATGCTGCCAGTGT‑‑GACCGGTAACTACGGCAGTTTGCTGATGCATACGCGCGCCAATGaa                                                                                               <  2:355433/132‑1 (MQ=255)
                                                    ggCAGACTTTCAGATGACGGGTATCCTGTTGAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGa                                                                                          <  2:581269/143‑1 (MQ=255)
                                                     gaaGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGaa                                                                                         <  1:804759/141‑1 (MQ=255)
                                                     gCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGATGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACt                                                                                       >  1:325721/1‑143 (MQ=255)
                                                        gACTTTCAGAGGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTcgc                                                                                    <  1:411788/143‑1 (MQ=255)
                                                           tttCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCACGCCAATGAAACAGAACCcgcc                                                                                   >  2:560944/1‑142 (MQ=255)
                                                              cAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCAGCCAGGAGCGGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGAcc                                                                                <  1:103642/143‑1 (MQ=255)
                                                                   gACGGGCATCGTGTTTAGTGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGAACGCCAtgt                                                                         <  2:37650/143‑1 (MQ=255)
                                                                     cGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGc                                                                       <  1:165251/143‑1 (MQ=255)
                                                                        gTATCCTGTTTAGAGACCGCCTGACGCAGCCAGGAGCCGATCCGCTCATGCGCTAGTTTATCCTGCCA‑‑GT‑‑GACCGGTCACAACCGGCTGTTGCTGACGCATACGCGCGCCAATGAAGCCGACATCGCGCCCGCCATCGGCAGt                                                                    >  2:12589/1‑143 (MQ=255)
                                                                         tATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑AT‑‑GACCGTTAACCACGGCATGTTGCGGCCGCCTAGACACGCCTATGCAGTCGAACTCGAGACAGCCAGGGGCCGTc                                                                   >  2:189233/1‑143 (MQ=255)
                                                                               gTTTAGAGACCGCCTGTCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACCCGCGACCGCCATGTGCAGTCAggt                                                               >  1:797292/1‑143 (MQ=255)
                                                                               gTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttc                                                             >  1:156515/1‑143 (MQ=255)
                                                                                  tAGGGACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcagg                                                          <  2:48172/143‑1 (MQ=255)
                                                                                             tGACGCGTCCACGAGCCGATACGCTCTTGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCCGTCACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTc                                               <  2:658274/143‑1 (MQ=255)
                                                                                                cgcgTCCCGGGGCCGATACGCTCATGGGCTTGTTTTGCCTGCCA‑‑GT‑‑GGCGGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCat                                            <  2:479264/143‑1 (MQ=255)
                                                                                                      cAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACAGCCCTGTGCAGTCTGGTGGAGGGTCATAAAGTCAATATCg                                        >  1:158096/1‑143 (MQ=255)
                                                                                                       aGGAGCCGATACGCTCATGGGCTAGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGAAGCATACGCGCGCCAGTGAAGCCGGACTGGCGACAGCGATGTCCACTGTGGTTCCGCTCCATAACGTCCGAATCGa                                       >  2:402741/1‑143 (MQ=255)
                                                                                                         caGCCGATACGCTCTTGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATa                                     <  1:270342/142‑1 (MQ=255)
                                                                                                          ggctgATTCGCGCATGGGCTTGTTTCTCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGATGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTg                                  <  1:376402/139‑1 (MQ=255)
                                                                                                             cGATACGCTCCTGGGCTTGTTTATCCTTCCA‑‑GT‑‑GCCCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTcc                               <  1:542117/143‑1 (MQ=255)
                                                                                                               aTACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACACGTGCAGAGGAATCGGACCTAGGGTCAGCGAAGGTCAGTCTGGTTCAGGTTCATAAAGTCCATCTCGATACTGTcc                               >  2:169112/1‑143 (MQ=255)
                                                                                                                         gggCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATCAAGTCCATATCGATCCTGTCCCACGGCCGcgcc                   >  1:28735/1‑143 (MQ=255)
                                                                                                                            cTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCcgcg                >  1:340864/1‑143 (MQ=255)
                                                                                                                                     ccTGCCG‑‑GT‑‑GCCCGGTACCCACGGCATGTTCCTGCCGCATACGCGGGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTgg       <  2:334940/143‑1 (MQ=255)
                                                                                                                                     ccTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCGTAGTACCGGTACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTgg       <  1:584693/143‑1 (MQ=255)
                                                                                                                                      cTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCCCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTggg      >  1:504742/1‑143 (MQ=255)
                                                                                                                                       tGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCGCGGCAGCGCCGCGGTGAGCTGGgt     >  1:701017/1‑143 (MQ=255)
                                                                                                                                          cA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCACGGCAATGAACCCGCACTCGCGACGGCCATGAGCCTTCGGGTTCAGCTTCATAAAGTCGAGATCGGTAGTGTAGCACGGCAGCCCCGAGCTGTACttcgtggt  >  2:310116/1‑136 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTG  >  NC_000913/67748‑68028

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