Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,946,3310CT47.8% 55.5 / 64.1 46noncoding (2628/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (8/14);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGG  >  NC_000913_3_hsa_tpiA/3946187‑3946377
                                                                                                                                                |                                              
ggggaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGgaa                                           <  4:58127/144‑2 (MQ=2)
ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGgag                                           >  1:587804/1‑150 (MQ=32)
  ggaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         <  1:124631/146‑1 (MQ=2)
  ggaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         <  1:311970/146‑1 (MQ=2)
  ggaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         >  1:559795/5‑150 (MQ=11)
  cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGaa                                         <  1:41523/150‑1 (MQ=32)
   ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAAc                                        >  3:373228/1‑150 (MQ=32)
   gaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaat                                        <  4:54918/147‑2 (MQ=2)
   gaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         <  4:230515/146‑1 (MQ=2)
   gaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         >  3:230515/4‑149 (MQ=2)
    gggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGgaaaatt                                       <  2:1027826/148‑6 (MQ=11)
    aggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGgaaaatt                                       <  4:726950/148‑6 (MQ=2)
    aggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGgaaaatt                                       <  3:96243/148‑6 (MQ=2)
    aggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaatt                                       <  3:7622/148‑3 (MQ=2)
     tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTg                                      >  2:582217/1‑150 (MQ=32)
     tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTg                                      >  2:51204/1‑150 (MQ=32)
       tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAg                                    >  1:59590/1‑150 (MQ=32)
       tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAg                                    <  2:595924/150‑1 (MQ=32)
           gCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGaaaaga                                >  3:317007/1‑145 (MQ=11)
           gCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGaaaaga                                >  1:382948/1‑145 (MQ=2)
           gCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGaaaaga                                >  3:360919/1‑145 (MQ=2)
             aCCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGaaaagagc                              >  1:947289/1‑143 (MQ=2)
                  ggtGTCGGCGCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGGGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTc                         <  4:589269/148‑1 (MQ=12)
                   aTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTcc                        <  2:757217/150‑1 (MQ=32)
                      tCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAg                     <  1:841488/150‑1 (MQ=32)
                        ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTa                   >  4:739443/1‑150 (MQ=32)
                         gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTAc                  <  3:739443/150‑1 (MQ=32)
                           tcatcgCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACga                >  1:1078398/1‑150 (MQ=2)
                           tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACga                <  1:409321/150‑1 (MQ=32)
                            catcgCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACgag               <  1:479126/150‑1 (MQ=2)
                            catcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACgag               >  4:94482/1‑150 (MQ=32)
                              tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGagag             >  2:25160/1‑150 (MQ=32)
                              tcaCATCCTGGAGCTGAAGTAGGGCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGagag             <  2:1226489/150‑1 (MQ=17)
                                 caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGga                                            >  2:794893/1‑116 (MQ=32)
                                 caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGga                                            <  1:794893/116‑1 (MQ=33)
                                   tCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCg        <  4:3051/150‑1 (MQ=2)
                                   tCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCg        <  3:401268/150‑1 (MQ=2)
                                     cTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGa      <  2:1165501/150‑1 (MQ=2)
                                     cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGGACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      >  4:113226/1‑150 (MQ=17)
                                     cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      >  4:407946/1‑150 (MQ=32)
                                     cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      >  4:407963/1‑150 (MQ=32)
                                     cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      <  4:430854/150‑1 (MQ=32)
                                      tGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaa     >  1:1065499/1‑149 (MQ=2)
                                       ggAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGt    <  1:916583/150‑1 (MQ=32)
                                        gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGACCTGTGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTg   <  3:359422/150‑1 (MQ=14)
                                         aGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGAATgg  <  2:1065499/150‑1 (MQ=2)
                                                                                                                                                |                                              
GGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGG  >  NC_000913_3_hsa_tpiA/3946187‑3946377

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