Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,544,022 1 bp→AA coding (39/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.A85.3% 96.9 / 11.7 34coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base A (15/14);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
*NC_000913_3_hsa_tpiA1,544,0220CA100.0% 112.0 / NA 34L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/17);  total (17/17)

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543879‑1544170
                                                                                                                                                 |||                                                                                                                                                    
ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑aaaaa                                                                                                                                                  <  1:320425/149‑2 (MQ=255)
        aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGa                                                                                                                                          <  1:108169/149‑1 (MQ=255)
          cGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACt                                                                                                                                        <  1:108474/149‑1 (MQ=255)
          cGGGAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCAAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCCGTCG‑CGAA‑AAAA‑GATACCGGCt                                                                                                                                        >  2:131222/1‑149 (MQ=255)
                         aTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTg                                                                                                                         <  1:566658/149‑1 (MQ=255)
                              aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    >  1:267928/1‑149 (MQ=255)
                                  aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCCGTGCCat                                                                                                                >  2:177599/1‑149 (MQ=255)
                                     cacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTGTTTTCCCAGAAGGCCGGATTTTCTGGTTTCGAGTCG‑CGAA‑AAAA‑GATAACGACTTTTTTTCTCATTTTGCAGTGCCATAtt                                                                                                             <  1:606345/149‑1 (MQ=255)
                                         gAAGGCAAGTAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  1:14264/149‑1 (MQ=255)
                                         gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         >  1:385531/1‑149 (MQ=255)
                                             gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     >  1:89306/1‑149 (MQ=255)
                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                   >  1:217194/1‑149 (MQ=255)
                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        >  1:20617/1‑149 (MQ=255)
                                                            ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                                         <  2:115895/146‑1 (MQ=255)
                                                            ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                                         >  1:115895/1‑146 (MQ=255)
                                                            ctcactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                      <  1:112427/146‑1 (MQ=255)
                                                             tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTTTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  1:101636/1‑149 (MQ=255)
                                                               actgaTATCCATGGATTTCTTCGAGCAATATGTTTGCCTTTAGTATCCCAGAAGGCCGGAGTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGTCTATTTTTCTCATTTTGCAGCGCCATATTGTTCCGCGCATGCACACATTGGTAAt                                                                                   <  2:585522/149‑1 (MQ=255)
                                                                            gaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      >  2:325274/1‑149 (MQ=255)
                                                                              tttCTTCGAGCAATATGTTTTCCTTCATGTTCACAGAAGGCCGGATTTTCTGGTTTCCAGTCGCCA‑A‑AAAA‑GATAAAGACTATTGTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    <  2:460187/149‑1 (MQ=255)
                                                                                              gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAAAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTATTGAAAAAAAGACAAAACACGAGGTAAGGCGCa                                                    >  1:435249/1‑149 (MQ=255)
                                                                                                  tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAAAGACTCTTTTTCTCATTTTGCAGTGCCATATTGTTCGCCTCATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAg                                                >  2:528576/1‑149 (MQ=255)
                                                                                                    cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                              >  1:43686/1‑149 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  1:533912/149‑1 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCACAATAGCCa                                             >  1:114696/1‑149 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATCCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  2:498473/1‑149 (MQ=255)
                                                                                                               cAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑AATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCTCAATAGCCAGTTATTAGaa                                   <  2:626209/149‑1 (MQ=255)
                                                                                                                   aGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:43686/149‑1 (MQ=255)
                                                                                                                    ggCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  2:126265/149‑1 (MQ=255)
                                                                                                                    ggCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCACAATAGCCAGTTATTAGAATTAAg                              <  2:114696/149‑1 (MQ=255)
                                                                                                                    ggCCGGATTTCTAGGTTTCCAGTCG‑CA‑A‑AAAATCATAACCACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  1:354538/149‑1 (MQ=255)
                                                                                                                       cGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTCTTTTTCTCATTTTGCAGTGCCATGTTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGCGGAAAGGCCAAATAGCCACGTATTAGAATTTAGGAt                           >  1:385134/1‑149 (MQ=255)
                                                                                                                       cGGATTTTCTGGTTTCCAGTCG‑CG‑AAAAAA‑GAGAACGACTGTTTTTCTCATTTTGCAGTGCAAAATTTTTCCCAACATGCAGACATTGGTAATGAAAAAAAGACAAAACAGTAGAAAAGGCGCACTAGCCAGCGATTAGAATTaacgcc                           >  2:142698/1‑145 (MQ=255)
                                                                                                                             ttCTGGTTTCCAGTCG‑CGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                     <  1:325274/149‑1 (MQ=255)
                                                                                                                                            tcgcgaaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        >  1:477725/8‑149 (MQ=255)
                                                                                                                                            tcgcgaaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        <  2:217194/142‑1 (MQ=255)
                                                                                                                                                gaaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa    >  1:508266/4‑149 (MQ=255)
                                                                                                                                                  aaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:95558/148‑1 (MQ=255)
                                                                                                                                                 |||                                                                                                                                                    
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543879‑1544170

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