Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,544,021:1 +A 68.0% coding (40/1389 nt) narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0211.A68.0% 51.9 / 23.5 25coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base A (10/7);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAT  >  NC_000913_3_hsa_tpiA/1543879‑1544145
                                                                                                                                                 |                                                                                                                              
ttGAGATTAACGGTAACTGCGCCAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaa                                                                                                                        <  1:390528/149‑2 (MQ=255)
 tGAGATTAACGGTAACTGCGCTAAATAGCCTCCAGACACAGAAGGCAAGAAGTAGACACCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCTAGAAGGCCGGAGTTTCTGCTCTCCAGTCGCGCataata                                                                                                                       >  2:347072/1‑147 (MQ=255)
      tgaaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTAATTATCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAAc                                                                                                                   <  1:266966/147‑1 (MQ=255)
                   cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtctc                                                                                                      >  1:423399/1‑149 (MQ=255)
                     cTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGAGTATTTTtctctt                                                                                                    >  1:47808/1‑147 (MQ=255)
                                      acaGAAGGCAAGAAGTAGACGACTGACCGACATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCAGTTTACAGTGCCATATTg                                                                                   >  1:671085/1‑149 (MQ=255)
                                         gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                <  1:451754/149‑1 (MQ=255)
                                            ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAAGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCTTTTTGCAGTGCCATTTTGTTCCTc                                                                             >  2:41732/1‑149 (MQ=255)
                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTTTCTTTTTTCCTATCCATAATGTTCCTaaca                                                                          >  2:129459/1‑149 (MQ=255)
                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                          >  1:293070/1‑149 (MQ=255)
                                                                 tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACACATTGGTAATGa                                                        >  1:280479/1‑149 (MQ=255)
                                                                      tcgagagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGCACACGTTGGTAATGtaaaaa                                                   >  2:571286/4‑149 (MQ=255)
                                                                      tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGAACACATTGGTAATGaaaaaa                                                   >  1:619686/1‑149 (MQ=255)
                                                                      tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCAGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCGCACATACACACATTGCTAATGaaaaaa                                                   >  2:93732/1‑149 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                >  2:559490/1‑149 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                <  1:83771/149‑1 (MQ=255)
                                                                           agaTTTTTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                              <  1:180221/149‑1 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                       <  1:150946/149‑1 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGATTCCAGTGGCTGAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                       >  1:366443/1‑149 (MQ=255)
                                                                                    cGAGTAATATGTTTACCTTTATTTATCCAGAAGGCCGGGTTTTGTGGTTTCCAGTCGCAAAAATAGATAACGACTATTTTTTTCATTTGGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg                                     <  2:621384/149‑1 (MQ=255)
                                                                                                tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                         >  1:167065/1‑149 (MQ=255)
                                                                                                    cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                     <  1:639513/149‑1 (MQ=255)
                                                                                                              ccaagatgCCGGATTTTTTGGTTTCCAGTCGGTAAAAAAGATAACGACTATTTTTCGCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa           <  2:436669/142‑1 (MQ=255)
                                                                                                              ccAGAGGGCCGGTTTGTGTGGTTTCCAGGCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTTCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa           <  2:436666/149‑1 (MQ=255)
                                                                                                                       cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGCGGTAAGGCGCAATAGCCAGTTATTAAAATTAAGGAt  >  2:657310/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                              
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAT  >  NC_000913_3_hsa_tpiA/1543879‑1544145

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