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.A83.3% 62.6 / 7.7 24coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base A (7/13);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 3.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_hsa_tpiA1,544,0220CA100.0% 72.4 / NA 24L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/13);  total (11/13)

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543879‑1544171
                                                                                                                                                 |||                                                                                                                                                     
ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGCTTCTAGTCG‑CG‑A‑aaa‑aa                                                                                                                                                  >  1:554229/1‑148 (MQ=255)
      ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGCTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGATTCCAGTGG‑CGAA‑AAA‑AGATAAc                                                                                                                                             >  1:481039/1‑149 (MQ=255)
        aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGa                                                                                                                                           <  1:97052/149‑1 (MQ=255)
                              aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                     >  2:94892/1‑149 (MQ=255)
                               tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCACTCG‑AGAA‑AAA‑ACATACCGAATATTTTTTTCATTTTGCAGTGc                                                                                                                    >  2:181164/1‑149 (MQ=255)
                                  aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATAGTTTCCTTTATTTTCCCAGAAGACTGGATTTTCTGGTTTCCAGTCGACA‑A‑AAA‑AGATAACGACTATTTTTCTAATATTCAAGTGCAat                                                                                                                 >  2:467557/1‑149 (MQ=255)
                                   aacacAGAATGCAAGAAGTAGACAACTCACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                <  1:440772/148‑1 (MQ=255)
                                               aagcagTAGACAACTCACTGATATCCAGAGATTTCTTCGAGCAATATGCTTTCCTTTATTTACCCAGAAGTCCGGATTTTCTGGTTGCCAGAAG‑AAAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    <  2:304263/149‑1 (MQ=255)
                                                        aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:197950/1‑149 (MQ=255)
                                                                 tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑AAAGAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                                  >  2:345380/1‑149 (MQ=255)
                                                                      tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATGTTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                             >  2:389509/1‑149 (MQ=255)
                                                                       ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  1:377273/149‑1 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          <  1:325876/149‑1 (MQ=255)
                                                                              tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                     <  1:225478/149‑1 (MQ=255)
                                                                              attcttcGAGCACAAGGTTTGCCTGTATTTTACCAGACGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGAGTTTCTTTCTCCTGTTTCAGGTCCATATTGTTCCTCACAGGCACACATGGGTAATGAAAAAAAGACaaaa                                                                     <  2:36187/148‑1 (MQ=255)
                                                                                 cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                  >  1:239468/1‑149 (MQ=255)
                                                                                   tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAGGACAAAACACGa                                                                >  1:94513/1‑149 (MQ=255)
                                                                                     gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCGGGATTTTCTGGTTTCGAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                              >  2:214620/1‑149 (MQ=255)
                                                                                           tatGTTTTCCTTTATTTTCCCAGAGGGCCGGTTTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                        <  1:121358/149‑1 (MQ=255)
                                                                                                   agtgtATTTTCCCAGGAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                                <  2:426023/145‑1 (MQ=255)
                                                                                                           ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                        <  1:118089/149‑1 (MQ=255)
                                                                                                                      ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                             <  1:540168/149‑1 (MQ=255)
                                                                                                                              tCTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAAAACGCACGATGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                     <  2:522934/149‑1 (MQ=255)
                                                                                                                               cTGGTTTCCAGTCG‑CGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                    <  1:215583/149‑1 (MQ=255)
                                                                                                                                                  aaaa‑aGATATCGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa   <  1:7109/148‑1 (MQ=255)
                                                                                                                                                   aaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  >  1:56297/1‑149 (MQ=255)
                                                                                                                                                 |||                                                                                                                                                     
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543879‑1544171

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