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.A77.4% 81.8 / 20.2 31coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base A (11/13);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
*NC_000913_3_hsa_tpiA1,544,0220CA96.8% 95.9 / NA 31L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_hsa_tpiA/1543877‑1544169
                                                                                                                                                  |||                                                                                                                                                   
tATTGAGATTAACGGTAACTACGCTAAATAGCATACAGACAAAGAAGGCAAGAAGTAGACAACTGGCTAGTATCCAGAGATTTCTTCGAACAATATGTTTCCCTTTATTTTCCTAGAAGGCAGGATTTTCTGGTTTGCAGTCGAG‑a‑aaa                                                                                                                                                   >  2:286843/1‑149 (MQ=255)
 aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑a                                                                                                                                                 <  2:352933/149‑1 (MQ=255)
         tAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AAAGCGATAACg                                                                                                                                          >  1:1361/1‑149 (MQ=255)
                aaCCGCGGTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTAttt                                                                                                                                   <  1:186054/149‑1 (MQ=255)
                   tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTtc                                                                                                                                >  2:243052/1‑149 (MQ=255)
                      gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCa                                                                                                                             >  2:325391/1‑149 (MQ=255)
                             aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAGGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCa                                                                                                                      <  1:340123/149‑1 (MQ=255)
                             aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCa                                                                                                                      <  1:90326/149‑1 (MQ=255)
                               cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                    >  2:585602/1‑149 (MQ=255)
                                        acaGAAGGCAAGAAGTAGACAACTGACTGAAATCCCGAGATTGCTTCGAGCAATATGTATTCCTTTATTTTCCCAGAGGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGAATATGTTTCTCATTTTGCAGTGCCATATTg                                                                                                           >  2:431975/1‑149 (MQ=255)
                                         caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                          <  1:498970/149‑1 (MQ=255)
                                               gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGCTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                    <  2:151731/149‑1 (MQ=255)
                                                cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aa                                                                                                                                                >  2:229315/1‑102 (MQ=255)
                                                cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aa                                                                                                                                                <  1:229315/103‑2 (MQ=255)
                                                 aagaagTAGACAACTGACTGATATCCAGAGATTTATTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTTTCGAA‑AAA‑AGATAACGACTATTTTTTTCATTTTTCACTGCCATATTGTTACTcgac                                                                                                  >  2:202867/1‑146 (MQ=255)
                                                      gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                             <  1:298662/149‑1 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                       <  1:525651/149‑1 (MQ=255)
                                                             actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                      <  1:325391/149‑1 (MQ=255)
                                                               tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                    >  1:120081/1‑149 (MQ=255)
                                                                          cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                         >  1:188464/1‑149 (MQ=255)
                                                                           agagATTTCTTTGAGGAAATTGTGTTCCTTTATTATCCAAAAATGCCGGATTTTCTGATTTCAATAAACAAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                        <  2:664336/149‑1 (MQ=255)
                                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGTAGTGCCATATTGTTCCTCCAATGCACACATTGGTAATGAAAAACAGa                                                                        >  2:675485/1‑149 (MQ=255)
                                                                              gaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                     <  1:322040/149‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                               <  1:243052/149‑1 (MQ=255)
                                                                                                 ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCCTTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                  >  1:113595/1‑149 (MQ=255)
                                                                                                 ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                  >  1:233282/1‑149 (MQ=255)
                                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                            >  2:614284/1‑149 (MQ=255)
                                                                                                                 cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATCGCCAGTTATTAGaa                                  >  2:182117/1‑149 (MQ=255)
                                                                                                                     aGGCGGGATTGTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                              <  2:723373/149‑1 (MQ=255)
                                                                                                                         cGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                          >  1:793432/1‑149 (MQ=255)
                                                                                                                          ggATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGATATTAGAATTAAGGATg                         >  1:356115/1‑149 (MQ=255)
                                                                                                                               ttCTCGTTTCCAGTCGCGAC‑AAA‑AGATACCGACTATTTTTCTCATTTTGCAGTGCCACATTGGTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                    <  1:302132/149‑1 (MQ=255)
                                                                                                                                              cgtgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt      <  1:286843/143‑1 (MQ=255)
                                                                                                                                               gcgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt     >  1:594347/6‑149 (MQ=255)
                                                                                                                                                  aaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAGAGATGTATTGGGTGAAGTGCTGATTGaa  >  2:541414/3‑149 (MQ=255)
                                                                                                                                                  |||                                                                                                                                                   
TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_hsa_tpiA/1543877‑1544169

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