Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,9050TC20.0% 81.0 / 15.0 40N568D (AAC→GAC) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (3/5);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTCGACGACATATTC  >  NC_000913/860759‑861052
                                                                                                                                                  |                                                                                                                                                   
cGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTg                                                                                                                                                    <  5:55526/148‑1 (MQ=255)
        cGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCg                                                                                                                                           >  5:17193/1‑149 (MQ=255)
              tCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     <  7:103485/149‑3 (MQ=255)
              tCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     <  6:8852/149‑3 (MQ=255)
                    gCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCa                                                                                                                               <  4:97188/149‑1 (MQ=255)
                        aCGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGaa                                                                                                                            <  6:17193/148‑1 (MQ=255)
                          gCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAAcc                                                                                                                         >  8:107209/1‑149 (MQ=255)
                               tgctCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGc                                                                                                                                 >  1:164323/1‑136 (MQ=255)
                               tgctCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGc                                                                                                                                 <  2:164323/136‑1 (MQ=255)
                                        gTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     <  4:5801/123‑3 (MQ=255)
                                        gTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     >  3:5801/1‑121 (MQ=255)
                                         tCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCt                                                                                                                                       <  1:39358/120‑1 (MQ=255)
                                         tCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCt                                                                                                                                       >  2:39358/1‑120 (MQ=255)
                                                       tCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGa                                                                                             >  2:147773/1‑148 (MQ=255)
                                                         aTCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGa                                                                                                                                          <  4:74458/101‑1 (MQ=255)
                                                         aTCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGa                                                                                                                                          >  3:74458/1‑101 (MQ=255)
                                                         aTCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATAc                                                                                          >  1:136714/1‑149 (MQ=255)
                                                             gCCAGTGCGGCAGCAAGCTGTTGCTGCCCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACtttt                                                                                      >  8:128140/1‑149 (MQ=255)
                                                              ccAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     >  4:13651/1‑99 (MQ=255)
                                                              ccAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTCGGCATTGCCAACCtcg                                                                                                                                     <  3:13651/101‑3 (MQ=255)
                                                               cAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATa                                                                                           <  5:130925/142‑1 (MQ=255)
                                                               cAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATa                                                                                           >  6:130925/1‑142 (MQ=255)
                                                                     ggCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTcgctcgct                                                                              >  3:78461/1‑149 (MQ=255)
                                                                                gTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATAc                                                                        <  1:123969/144‑1 (MQ=255)
                                                                                gTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATAc                                                                        >  2:123969/1‑144 (MQ=255)
                                                                                gTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTc                                                                    >  4:89883/1‑148 (MQ=255)
                                                                                       aCCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATccac                                                              <  7:107209/147‑1 (MQ=255)
                                                                                                     gTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCaga                                                <  3:89883/147‑1 (MQ=255)
                                                                                                        cAAACACCAGTTTCTTCACTGCAGCCAGGATGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTAGATCGCTCAATAAAGTCATCCACCAGCGCCGAGCAGAGAata                                           >  8:87564/1‑149 (MQ=255)
                                                                                                         aaaCACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAatat                                          >  7:106842/1‑149 (MQ=255)
                                                                                                                aGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCAc                                   <  7:87564/149‑1 (MQ=255)
                                                                                                                     cttcACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTtc                               <  1:147773/148‑1 (MQ=255)
                                                                                                                           tGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTc                          >  2:33277/1‑147 (MQ=255)
                                                                                                                           tGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTc                          <  1:33277/147‑1 (MQ=255)
                                                                                                                                 cAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATc                                         <  2:23526/126‑1 (MQ=255)
                                                                                                                                 cAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATc                                         >  1:23526/1‑126 (MQ=255)
                                                                                                                                  aGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGgt                 <  8:106842/149‑1 (MQ=255)
                                                                                                                                       gTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTcga            <  7:128140/149‑1 (MQ=255)
                                                                                                                                          gCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTcgacga         >  3:103016/1‑149 (MQ=255)
                                                                                                                                               ggTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTCGACGACatat    <  4:78461/149‑1 (MQ=255)
                                                                                                                                                 ttGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTCGACGACATAttc  <  2:136714/149‑1 (MQ=255)
                                                                                                                                                 ttGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTCGACGACATAttc  <  4:103016/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCAGGCCGTCGAAGTCATCTGCCAGTGCGGCAGCAAGCTGTTGCTGACCAATCGCACCTTGTTCAAACACCAGTTTCTTCACTGCCGCCAGGCTGTTGCCGAGGTTGGCAATGCCGACCTGCAGGCCAGAAACCCAGTCATATTTCGCGCCGCCTTGCTTGATACTTTTCGCTCGCTCAATACAGTCATCCACCAGCGCCGAGCAGAGAATATCGTGCACGTTCTCTTCCAGCATGGTGTCGACGACATATTC  >  NC_000913/860759‑861052

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