Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0221.A33.3% 31.2 / 38.8 27coding (39/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (8/10);  new base A (5/4);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_hsa_tpiA/1543894‑1544157
                                                                                                                                  |                                                                                                                                       
aCTGCGCTAAATAGCATCTAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGCTTTCTTCGAGCGATATGTTTTCCTTTATTTGCCCATAAGGCCGGATTTTCTAGTTTCCAGTCGCG‑CAAAAAGATAACGATTAtttt                                                                                                                      >  2:430410/1‑149 (MQ=255)
                       acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCTAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                               >  1:247105/1‑149 (MQ=255)
                          gagggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                            <  2:527339/146‑1 (MQ=255)
                          gAAGGAAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                            <  1:494314/149‑1 (MQ=255)
                             ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                         <  2:445812/149‑1 (MQ=255)
                                 agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTGTTTTCCCAGAAGGCCGGATTTGCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                     >  1:540493/1‑149 (MQ=255)
                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                              >  1:534789/1‑149 (MQ=255)
                                          cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                            >  2:112613/1‑149 (MQ=255)
                                                         cAGAGATGTCTGCTAGCAATATGTTGTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTGTCATTTTGCAGTGGCCAAGTTTGCCTCACATGCACCCATTGGTAATGAAAAAAAg                                                             >  2:559957/1‑149 (MQ=255)
                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAACAGGCCGGATTTTCTGGACTCCTGCCGCG‑CAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATTAAACAAAGa                                                            >  1:522508/1‑149 (MQ=255)
                                                           gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                           >  2:226356/1‑149 (MQ=255)
                                                               tttCTTAGAGCAAAATGTTTTCCTTTATTTTACCAGAAGAAAGAAAAATCTGGTTAACAGTCGCG‑AAAAAAGATAACGACTATTTTCCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACAAAAAAAAAAAAAAAAAGACaaaa                                                       <  2:491052/149‑1 (MQ=255)
                                                                ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCA‑TAAAAAGATAACGACTACGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                      <  1:527819/149‑1 (MQ=255)
                                                                       aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGCAAAAAAGACAAAACACGAGGt                                               >  1:58605/1‑149 (MQ=255)
                                                                        gCAATTTGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTATGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                              <  2:352398/149‑1 (MQ=255)
                                                                        gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                              <  1:218305/149‑1 (MQ=255)
                                                                               gTTTTCCTTTATTTCCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCa                                       <  1:76507/149‑1 (MQ=255)
                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAAACCAAACACGAGGTAAGGCGCAATAGCCa                                >  2:24556/1‑149 (MQ=255)
                                                                                         aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                             >  1:283818/1‑149 (MQ=255)
                                                                                                 aGACGGCCGGATTCTCTGGTTTCCAGTCGCGAC‑ACAAGATCCCGTCTTTTTTTCTCATTTTGCAGTGCCACATTGTCCCTCACATGCACAAATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                     <  1:559957/149‑1 (MQ=255)
                                                                                                     ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                 <  1:332710/149‑1 (MQ=255)
                                                                                                     ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                 <  1:306418/149‑1 (MQ=255)
                                                                                                      gccgttttttCTTTTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGTTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg                >  2:36115/7‑149 (MQ=255)
                                                                                                                    tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGCACACATTGGAAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGac  <  2:540493/149‑2 (MQ=255)
                                                                                                                    tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                  >  2:410741/1‑117 (MQ=255)
                                                                                                                    tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                  <  1:410741/117‑1 (MQ=255)
                                                                                                                    tttCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa  <  1:112613/149‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
ACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_hsa_tpiA/1543894‑1544157

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