Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,801,0560AT21.7% 85.0 / 23.6 46N31I (AAT→ATT) nrdINrdEF cluster assembly flavodoxin
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base T (5/5);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGC  >  NC_000913_3_hsa_tpiA/2800910‑2801203
                                                                                                                                                  |                                                                                                                                                   
gtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                                   >  2:960520/1‑149 (MQ=255)
   tccgtccgGACATTATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGc                                                                                                                                                <  2:280343/149‑1 (MQ=255)
     cgtccgCACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCgg                                                                                                                                              <  2:500940/149‑1 (MQ=255)
        ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTGTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCACATCCCGCTCATTGAGCGGGat                                                                                                                                           <  2:351467/149‑2 (MQ=255)
          ggACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACg                                                                                                                                         >  2:176707/1‑149 (MQ=255)
                 aTTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                                  <  1:960520/149‑1 (MQ=255)
                   tAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGt                                                                                                                                <  2:824842/149‑1 (MQ=255)
                     aCCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAg                                                                                                                              <  2:738710/149‑1 (MQ=255)
                      ccGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           >  1:127489/1‑134 (MQ=255)
                      ccGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           <  2:127489/135‑2 (MQ=255)
                       cGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGAc                                                                                                                            >  1:341077/1‑149 (MQ=255)
                           gggggccAGGGCAACACGCGGCCAGTGTATGAGGCAGCTTGTCTACTTCTCCAGGAGCTCCGAAAACACGCGGCGTTTTAGCGAACGTTGAGGTCTGCGCGCGGTGCGCATCCGGCTCAATGAGCGTGAACGGATTCAGGTAGACGAGc                                                                                                                        <  2:157564/144‑1 (MQ=255)
                            gCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAAACGAGcc                                                                                                                       >  2:926118/1‑149 (MQ=255)
                                     gcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCt                                                                                                              >  2:50426/1‑149 (MQ=255)
                                       cacaCGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGa                                                                                                            <  1:176707/149‑1 (MQ=255)
                                         cacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGCATGACGGTATGGCGCATCCCGCTCATAAAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                          <  2:924228/149‑1 (MQ=255)
                                          acGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCg                                                                                                         >  1:795900/1‑149 (MQ=255)
                                                 cAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCctc                                                                                                  <  2:337697/149‑1 (MQ=255)
                                                     gCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCtttc                                                                                              >  2:345347/1‑147 (MQ=255)
                                                     gCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAc                                                                                              <  2:707057/149‑1 (MQ=255)
                                                       aggAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAACCACGCCGCGTTTTAGCGAACGTTTAGGGCTCCCCGCGGTGCACATCCCGCTCAATGAGCGGAAACGAGATCAGGTGAACGAGCGTTTCATCCGGATCGTGCGCTCttcagg                                                                                            >  2:137720/3‑145 (MQ=255)
                                                               gCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                          >  1:432802/1‑127 (MQ=255)
                                                               gCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                          <  2:432802/127‑1 (MQ=255)
                                                                    tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGCCAGCCGTAcc                                                                               >  2:949068/1‑148 (MQ=255)
                                                                      tACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGGACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGGCGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGc                                                                             >  2:533387/1‑149 (MQ=255)
                                                                         ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTgg                                                                          >  2:318960/1‑149 (MQ=255)
                                                                            tcCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           >  1:103322/1‑80 (MQ=255)
                                                                            tcCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           <  2:103322/81‑2 (MQ=255)
                                                                             ccAGCATGTCCGAAACCCCGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATTCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGTcgcg                                                                      <  2:734541/149‑1 (MQ=255)
                                                                                    cTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCAc                                                               <  1:602010/149‑1 (MQ=255)
                                                                                           aaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGt                                                        >  1:213764/1‑149 (MQ=255)
                                                                                             cacGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTaa                                                      >  1:549137/1‑149 (MQ=255)
                                                                                                gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                                   >  1:594484/1‑149 (MQ=255)
                                                                                                     gtgttATCGAACGTTTAGGTCGGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGGGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCttt                                              <  2:644730/146‑1 (MQ=255)
                                                                                                         tATCGAACGCTTAGGTCTGCCCGCGGTGAGCATCCCGCTCAATGAGCGGGAACGAATTCAGGTAGCCGAGCCTTACATCCTGACCGTACCCTCTTACGGCGGCGGCCGTACGACTTACGCCGAGCCACGATAGCGAGTGCGCTTTTCaa                                          >  2:267058/1‑149 (MQ=255)
                                                                                                          aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                         >  1:301764/1‑149 (MQ=255)
                                                                                                          aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCATACCTCCTGATCGTGCCCTCTTACGGCTGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                         <  1:911260/149‑1 (MQ=255)
                                                                                                                  tttAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGGTCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAGTTCGCTTTTTAACcgccgacc                                 >  1:595701/1‑147 (MQ=255)
                                                                                                                  tttAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTACACGACGAGc                                 >  2:673881/1‑149 (MQ=255)
                                                                                                                  tttAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGc                                 >  2:852342/1‑149 (MQ=255)
                                                                                                                          tGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCgg                         <  1:945503/149‑1 (MQ=255)
                                                                                                                               gcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGAGCTCCGCCAGGCCCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCCGGCGtt                    >  1:225336/1‑149 (MQ=255)
                                                                                                                                   tGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt                >  2:764834/1‑149 (MQ=255)
                                                                                                                                   atccgttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt                >  2:656876/7‑149 (MQ=255)
                                                                                                                                      cgttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgc             >  2:259593/4‑149 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg            >  1:282276/3‑149 (MQ=255)
                                                                                                                                                 aaTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGc  >  1:782447/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGC  >  NC_000913_3_hsa_tpiA/2800910‑2801203

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