Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,474,8131.T26.7% 20.3 / 15.4 15coding (593/861 nt)ydbCputative NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (7/4);  new base T (2/2);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAACAACAAGGCCTGGTAAAACATATTGGCCTGAGCAACGTCACACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGT  >  NC_000913/1474671‑1474941
                                                                                                                                               |                                                                                                                                
caacaacaaGGCCTGGTAAAACATATTGGCCTGAGCAACGTCTCACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGc                                                                                                                            <  8:12379/149‑1 (MQ=255)
                         ttGGCCTGAGCAACGTCACACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCg                                                                                                   >  7:105850/1‑149 (MQ=255)
                             ccTGAGCAACGTCACACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATCCTTTGGCCCACGATGGCATTGCCTACGTGCCGttct                                                                                               >  1:141898/1‑149 (MQ=255)
                                        tCACACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCggg                                                                                    <  8:105850/149‑1 (MQ=255)
                                            acCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCGGATGATGCAATGATTGATGCTTTGCCCCACGATGGCATTGCCTACTTGCCGTTCTTCCCGCTCGGGGGCt                                                                                >  3:86923/1‑149 (MQ=255)
                                                aCGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTa                                                                             >  6:4931/1‑149 (MQ=255)
                                                aCGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTa                                                                             >  6:7036/1‑149 (MQ=255)
                                                aCGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTa                                                                             >  6:86058/1‑148 (MQ=255)
                                                    aGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCACGGGGGCTTTACACCg                                                                        >  4:112439/1‑149 (MQ=255)
                                                                             cGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTa                                                                             <  6:135425/119‑1 (MQ=255)
                                                                                aaTTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGtt                                             <  5:4931/149‑1 (MQ=255)
                                                                                aaTTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGtt                                             <  5:7036/149‑1 (MQ=255)
                                                                                       tgtgtgCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAg                                     >  6:29342/1‑149 (MQ=255)
                                                                                                             tCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAg               >  2:54414/1‑149 (MQ=255)
                                                                                                                            atgatgCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTGCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGt  <  2:141898/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                
CAACAACAAGGCCTGGTAAAACATATTGGCCTGAGCAACGTCACACCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGT  >  NC_000913/1474671‑1474941

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: