Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,158,8541.G27.3% 40.0 / 25.5 22coding (986/1434 nt)ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base . (12/4);  new base G (6/0);  total (18/4)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAA  >  NC_000913_3_hsa_tpiA/1158707‑1158999
                                                                                                                                                    |                                                                                                                                                 
cgggCGGGACCATGAGCTCCGCGGCGCAGACCTCGTTCCTTACCGGTATCACCGAGCCGATCGAGTTTTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATtcttctt                                                                                                                                                   <  2:547049/148‑2 (MQ=255)
       tattatGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑gggggatg                                                                                                                                           >  1:313755/1‑145 (MQ=255)
                    gcggcgCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTAcgtc                                                                                                                              >  1:283614/1‑149 (MQ=255)
                       gcgtTGACTCCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGtt                                                                                                                           <  2:56168/149‑1 (MQ=255)
                                        gACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCTCGCGACTCTGGCAGGCCTGACATTCCCAATCTTTATTCTACT‑GGGGATGCGGGACGGTACGTCGTTCTCGGACGGGCTGATCg                                                                                                          >  2:488608/1‑149 (MQ=255)
                                          ccGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGAc                                                                                                        <  1:451953/149‑1 (MQ=255)
                                                   ccGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGtt                                                                                               >  1:317469/1‑149 (MQ=255)
                                                          gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                        >  1:545314/1‑149 (MQ=255)
                                                                      cTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACagcagc                                                                             >  2:318033/1‑149 (MQ=255)
                                                                          aTGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAAca                                                                        >  1:96062/1‑148 (MQ=255)
                                                                          aTGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACt                                                                        >  1:97817/1‑149 (MQ=255)
                                                                           tGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACt                                                                        >  1:94575/1‑149 (MQ=255)
                                                                            gTTCGTTGCGCCGATCCTGTACATCACCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACtgt                                                                      <  2:223563/149‑1 (MQ=255)
                                                                             atCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑CGGGATGCGTGACGCTACGTCGTTCTCGCACGGTCTGATCAACTTCATCGTTCTGTCTGGTTACAGCAGCAAACtgtg                                                                     >  2:526852/2‑149 (MQ=255)
                                                                              tCGTTGCGCCGATCCTGTACATCATCCCCGCGATTCTGGCAGGTCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTgg                                                                    <  1:224352/149‑1 (MQ=255)
                                                                                     gcCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTcc                                                             >  1:27794/1‑149 (MQ=255)
                                                                                                    catcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGCCTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATc                                               >  1:212172/1‑149 (MQ=255)
                                                                                                    catcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATc                                               >  1:495937/1‑149 (MQ=255)
                                                                                                     atcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCg                                              >  1:369727/1‑149 (MQ=255)
                                                                                                                      aGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGAGCGTCGGTATCGGTTATGCGATTGTTtact                            >  1:88123/1‑149 (MQ=255)
                                                                                                                               aTTCCCAATCTGTATTCTTCT‑GGGGCTGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATGGGTTGTGCGATTGTTTTCTTCCGCATCt                   >  1:227343/1‑149 (MQ=255)
                                                                                                                                 tCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTGCTACCGCATCTTc                 >  1:293419/1‑149 (MQ=255)
                                                                                                                                    cAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGATCCTGCCGATCGTCGGTATCGGTTATGCGATTGTTTACTCCACCATCTTCcgc              >  1:157977/1‑149 (MQ=255)
                                                                                                                                         tGTATTCTTCTGGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGACCGACTTCATCGTCCTGTCTGGTAACAGCACCAAAATGTGGCTGTTCCCGATCGTCGCTAGCGGTTATGCGATTGTTTACTCAACCATCTTCCGCGTGc          >  2:50667/1‑149 (MQ=255)
                                                                                                                                                ttct‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTCTCGGTTACGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTtaa  >  2:465044/1‑149 (MQ=255)
                                                                                                                                                ttct‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTaaa  <  1:285969/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCT‑GGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAA  >  NC_000913_3_hsa_tpiA/1158707‑1158999

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