Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0211.A73.0% 94.3 / 30.5 37coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base A (12/15);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGT  >  NC_000913_3_hsa_tpiA/1543878‑1544159
                                                                                                                                                  |                                                                                                                                            
aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCCGGCACATAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaa                                                                                                                                       >  1:203423/1‑149 (MQ=255)
    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaagat                                                                                                                                   >  1:59842/1‑145 (MQ=255)
                  tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTTTTTTCCCAGTAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtc                                                                                                                      >  1:160851/1‑149 (MQ=255)
                  tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtc                                                                                                                      >  1:123272/1‑149 (MQ=255)
                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCa                                                                                                                   <  1:607614/149‑1 (MQ=255)
                            aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATTCAGAGATTTATTCGATCAATATGTCTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                            <  2:637914/149‑1 (MQ=255)
                               aTCCAGACACAGAATTCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                         >  2:525980/1‑149 (MQ=255)
                                tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTTTGGGCCCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                        >  1:821365/1‑149 (MQ=255)
                                 ccAGACACAGAAGGCAAGAAGTAGACCACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTGCTCATTTTGCAGTGcc                                                                                                       <  1:90187/149‑1 (MQ=255)
                                        caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAGTTTCCCAGAAGGCCGGATTTTCTTGTTTCCATTCGCGAAAAAAGATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGt                                                                                                <  2:96447/149‑1 (MQ=255)
                                                  gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                      <  2:123272/149‑1 (MQ=255)
                                                   aagTAGACAACTGACTGATATCCAGAGCTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAATAAGATAACGACTATTTTTCTCCTTTTGCCGTCCCATCTTGTTCCTCAAATgt                                                                                     >  2:450141/1‑148 (MQ=255)
                                                    agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                    <  2:59842/149‑1 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                   <  1:132138/149‑1 (MQ=255)
                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                          >  1:599804/1‑149 (MQ=255)
                                                               tactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATACCCTCTATCGTGCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                         <  2:770037/148‑1 (MQ=255)
                                                                       tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                 >  1:368457/1‑149 (MQ=255)
                                                                                   ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                     >  1:214628/1‑149 (MQ=255)
                                                                                    tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                    <  2:368457/149‑1 (MQ=255)
                                                                                        gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                            >  1:440174/1‑121 (MQ=255)
                                                                                        gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                            <  2:440174/121‑1 (MQ=255)
                                                                                             atGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                            <  1:371756/116‑1 (MQ=255)
                                                                                             atGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                            >  2:371756/1‑116 (MQ=255)
                                                                                              tGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGCCGCCTGGCGAGAGAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcgc                                          <  1:450141/149‑1 (MQ=255)
                                                                                                ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAAAACGGGGTAAGGGGCaa                                        >  2:503994/1‑149 (MQ=255)
                                                                                                 tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                       >  2:746875/1‑149 (MQ=255)
                                                                                                  ttCCTTTATTTTCCCATAAGGCCGGATTTTCGGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                      <  1:767522/149‑1 (MQ=255)
                                                                                                            ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                              <  1:628216/67‑1 (MQ=255)
                                                                                                            ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                              >  2:628216/1‑67 (MQ=255)
                                                                                                            ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                            <  1:682814/149‑1 (MQ=255)
                                                                                                               ctagagGGCCGGATTTTCTGGTTTCCAGTGGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTACTCACATGCACATATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                         <  2:203423/147‑1 (MQ=255)
                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCAAAAAAAGATAACGACTATTTTTCTCATTTTGCATTGCCATATTGTTCCTCACATGCACCCATTGGTAATGAAAAAAAGACAAAACACAAGGTAAGGCGCAATAGCCAGTTATTAGa                         >  1:426707/1‑149 (MQ=255)
                                                                                                                     ggCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                   <  2:214628/149‑1 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             <  2:260952/77‑1 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             >  1:260952/1‑77 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa     >  1:684406/1‑149 (MQ=255)
                                                                                                                                      tCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGt  <  1:191253/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGT  >  NC_000913_3_hsa_tpiA/1543878‑1544159

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: