Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,756,4970GT63.6% 30.2 / 34.4 33C267F (TGT→TTT) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (5/7);  new base T (12/9);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCC  >  NC_000913_3_bme_pgi/1756354‑1756642
                                                                                                                                               |                                                                                                                                                 
cAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCg                                                                                                                                             >  2:377232/1‑150 (MQ=255)
cAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCg                                                                                                                                             >  2:176620/1‑150 (MQ=255)
           aCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAg                                                                                                                                  >  4:83374/1‑150 (MQ=255)
            ccAGGAAGGCCTCAACAACTTCGACGAAATCGTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGt                                                                                                                                 >  2:219252/1‑150 (MQ=255)
                       tCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                                                      <  3:259790/150‑1 (MQ=255)
                        caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                                                                                     <  1:176620/150‑1 (MQ=255)
                              cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                               <  1:39941/150‑1 (MQ=255)
                                                aGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGAt                                                                                             <  1:377232/150‑1 (MQ=255)
                                                  ccTCTGACGGCATCATGGGTGCGCGTGGCGCCCGGGGAGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAACATGATGATCGAAAAATGTATCCGTGCACGAAAAGACGTTATCACTGCGACCCCGATGCTGGATTCCATGaacg                                                                                           >  4:130019/1‑147 (MQ=255)
                                                              tcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGAt                                                                                                            >  3:238391/1‑121 (MQ=255)
                                                              tcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGAt                                                                                                            <  4:238391/121‑1 (MQ=255)
                                                                                  cTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCa                                                                                                 <  4:154282/112‑1 (MQ=255)
                                                                                  cTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCa                                                                                                 >  3:154282/1‑112 (MQ=255)
                                                                                        gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGtt                                                     <  1:219252/150‑1 (MQ=255)
                                                                                         tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGAGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                    >  4:276398/1‑150 (MQ=255)
                                                                                                       gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATc                                         >  3:137979/1‑147 (MQ=255)
                                                                                                       gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATc                                         <  4:137979/147‑1 (MQ=255)
                                                                                                         agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                                       >  4:238863/1‑147 (MQ=255)
                                                                                                         agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                                       <  3:238863/147‑1 (MQ=255)
                                                                                                             gTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt                                >  2:360891/1‑150 (MQ=255)
                                                                                                                         cAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              >  4:241770/1‑44 (MQ=255)
                                                                                                                         cAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              <  3:241770/44‑1 (MQ=255)
                                                                                                                           gaagaTGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAt                                                                                                >  1:235603/1‑72 (MQ=255)
                                                                                                                           gaagaTGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAt                                                                                                <  2:235603/72‑1 (MQ=255)
                                                                                                                                tgatgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGt             <  3:83374/150‑1 (MQ=255)
                                                                                                                                 gatgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCt                                                                                                         <  1:97595/57‑1 (MQ=255)
                                                                                                                                 gatgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCt                                                                                                         >  2:97595/1‑57 (MQ=255)
                                                                                                                                   tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg          <  4:74341/150‑1 (MQ=255)
                                                                                                                                     atCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt        >  2:105028/1‑150 (MQ=255)
                                                                                                                                         aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGTTGTCTGGTGAAt    >  2:287036/1‑150 (MQ=255)
                                                                                                                                         aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCATATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             >  2:286069/1‑93 (MQ=255)
                                                                                                                                         aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGACGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             <  1:286069/93‑1 (MQ=255)
                                                                                                                                           aaaTGTATCCGTCCACGTAAATTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTTACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc  <  1:102955/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCC  >  NC_000913_3_bme_pgi/1756354‑1756642

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: ‑