Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0211.A70.6% 77.2 / 31.1 34coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (7/3);  new base A (13/11);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATG  >  NC_000913_3_hsa_tpiA/1543874‑1544146
                                                                                                                                                    |                                                                                                                               
ttttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcga                                                                                                                                 <  2:521720/149‑2 (MQ=255)
     ttGAGATGAAAGGTACCGGCTCATAATAGTATCCATATACAGAGAGCAAGAATAAGATAACTACCTGATATCCAGAGATGTCTTTGAAAAATTTGTTTTCCTTCATTTTCCCAGAAGGCCTGATTTTCTGGTTTCCAGTCGCG‑A‑a‑aaaa                                                                                                                         <  2:681484/149‑2 (MQ=255)
        agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑a‑aaaagat                                                                                                                      >  2:747338/1‑145 (MQ=255)
         gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑ACAAAGAtg                                                                                                                      >  1:226817/1‑148 (MQ=255)
          actAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATaa                                                                                                                     <  1:663109/147‑1 (MQ=255)
             aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGa                                                                                                                  <  1:352621/149‑1 (MQ=255)
                    aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTAtttt                                                                                                           >  2:24355/1‑149 (MQ=255)
                    aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTAtttt                                                                                                           >  1:717950/1‑149 (MQ=255)
                     cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTAttttt                                                                                                          <  1:747338/149‑1 (MQ=255)
                      tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAGGGCCGGATTTTATGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTtc                                                                                                         >  1:517187/1‑149 (MQ=255)
                      tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTtc                                                                                                         >  2:550676/1‑149 (MQ=255)
                                 gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAg                                                                                              <  2:250327/149‑1 (MQ=255)
                                               aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTTGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                >  2:74807/1‑149 (MQ=255)
                                                aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                               <  1:274886/149‑1 (MQ=255)
                                                     agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                          <  1:431973/149‑1 (MQ=255)
                                                            gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                   >  1:275758/1‑149 (MQ=255)
                                                             aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATTTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                  <  1:427361/149‑1 (MQ=255)
                                                             aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                 <  1:510202/134‑1 (MQ=255)
                                                             aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                 >  2:510202/1‑134 (MQ=255)
                                                               aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                <  1:24355/149‑1 (MQ=255)
                                                               aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTAGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTTTTTTTTTTCTTTTCCATATTGTTCCTCACATGCACACATTg                                                                >  2:449146/1‑149 (MQ=255)
                                                                      tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                         >  1:6317/1‑149 (MQ=255)
                                                                           tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTTTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                    >  1:361675/1‑149 (MQ=255)
                                                                           tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATCgaaaaa                                                    >  2:65497/1‑149 (MQ=255)
                                                                              agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                 >  2:323390/1‑149 (MQ=255)
                                                                              agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑CAA‑AAAGATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACATGCACACTTTGGTAATGAAAAAACGa                                                 >  1:334780/1‑149 (MQ=255)
                                                                              agagATTGCTGCCTGCAATATGGTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                 <  2:226817/149‑1 (MQ=255)
                                                                                  aTTTTTTCGAGCAATATTTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                             <  2:72053/149‑1 (MQ=255)
                                                                                  aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                             >  1:195285/1‑149 (MQ=255)
                                                                                       ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                        >  2:213829/1‑149 (MQ=255)
                                                                                       ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTAGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCTTATTGTTCCTCACATGCACACATTGGTAATGATGATATGATAAAacagg                                        >  2:212977/1‑147 (MQ=255)
                                                                                                          tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                     >  2:635270/1‑149 (MQ=255)
                                                                                                                   ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa            <  2:560576/149‑1 (MQ=255)
                                                                                                                         ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg      <  1:156926/149‑1 (MQ=255)
                                                                                                                             ggATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg  >  1:417017/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                               
TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATG  >  NC_000913_3_hsa_tpiA/1543874‑1544146

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