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.A75.0% 68.4 / 20.6 28coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (1/6);  new base A (13/8);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 7.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_hsa_tpiA1,544,0220CA100.0% 87.0 / NA 28L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/14);  total (14/14)

TTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543876‑1544171
                                                                                                                                                   |||                                                                                                                                                     
ttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGaaa                                                                                                                                                      >  1:749626/1‑149 (MQ=255)
      agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCTCGAaaaaagat                                                                                                                                                <  1:513343/149‑5 (MQ=255)
                     gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtct                                                                                                                                 <  2:776265/149‑1 (MQ=255)
                       gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCa                                                                                                                               <  2:597683/149‑1 (MQ=255)
                              aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        >  1:41461/1‑149 (MQ=255)
                                    cAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAATGTTAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                                  <  2:471234/149‑1 (MQ=255)
                                            gAAGGAAAGAAGTAGACAACTTACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          <  2:597165/149‑1 (MQ=255)
                                                  aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    >  2:296203/1‑149 (MQ=255)
                                                   ataagTAGGCAACTGACTGATATCCAGAGATTTCTTCGATCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                   <  2:711712/147‑1 (MQ=255)
                                                     aagTAGACAACTGACTGATATCCAGAGATCTCTTCGAGGAATATATTTTCGTTTCTTGTCCGAGAAGGCCCGAGGTTCTGGTTTCCAGTCGCGAAAAAAGATAGCGACTATTTTTCTCATTTTGCGGTGCCATATTGTTCCTCAGATGc                                                                                                 >  1:705294/1‑149 (MQ=255)
                                                                tgactgaGATCCAGGGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAATGTTCCTCGCATGCACCCATTGTTa                                                                                      <  2:184221/149‑1 (MQ=255)
                                                                 ttctgaTATACAGAGATTTCTTCGAGCAATATGTTTTTTTTTATTTTGAGGGGGTACTACATGTTCTGGTTTACAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     <  2:388129/147‑1 (MQ=255)
                                                                 gcctgaTATCCAGAGAGTTCTTTGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     <  2:9629/147‑1 (MQ=255)
                                                                 gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCATCACATGCAAACATTGGTaa                                                                                     >  1:552351/1‑149 (MQ=255)
                                                                    tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAtag                                                                                  >  1:388129/1‑147 (MQ=255)
                                                                          acAGAGATTGCTTCGAGCAATATGTTTTCATTTATTTTCCCAGAAGTCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  1:298347/148‑1 (MQ=255)
                                                                               gaTTACTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                       >  2:635066/1‑149 (MQ=255)
                                                                                  ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACAGATTGGTAATGAAAAAAAGACAAAac                                                                    <  1:635066/149‑1 (MQ=255)
                                                                                            aaTATTTGTTCCTTTATTTGCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                          <  1:261983/149‑1 (MQ=255)
                                                                                            aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                          <  1:83060/149‑1 (MQ=255)
                                                                                                  ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    >  1:260364/1‑149 (MQ=255)
                                                                                                  ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    >  1:570829/1‑149 (MQ=255)
                                                                                                        tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  1:665890/1‑149 (MQ=255)
                                                                                                        tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAAGTAAGGAGCAAAAGcgc                                              >  2:771923/1‑147 (MQ=255)
                                                                                                         ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTCTTTTTCTCATTTTGCAGTGCCATATTGTTCCTAACAACCAACCACTGGCAACGAAAAATAGACAAAACAAGAGTGAAGGCGCAATAGccgg                                             >  2:713699/1‑147 (MQ=255)
                                                                                                                          cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTCTTTTTCTCATTTTGCCGTGCCCTATTGTTCCTCCCATGCACCCATTGGTAATGAAAAAAAGACAAAACACGAGGAAAGGCGCAATCGCCCGTTATTAGAATTAAGGaa                            >  2:753933/1‑148 (MQ=255)
                                                                                                                             aTTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAACACACGAGGTAACGCGCAATAGCCAGTTATTAGAATTAAGGATGaa                         >  2:790769/1‑149 (MQ=255)
                                                                                                                                 tCTGGTTTCCAGTCGCGAAAAAAGATAACGCCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCCCACATTGGTAATGAAAAAAAGACAAAAGACGAGGTAAGGCGCAATAGCCAGTTGTTAGAATTAAGGGTGAATTgg                     >  1:224418/1‑149 (MQ=255)
                                                                                                                                   tGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                   <  2:260364/149‑1 (MQ=255)
                                                                                                                                              tcgagaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  1:673322/142‑1 (MQ=255)
                                                                                                                                              tcgagaaaaaaGATAACGACTAGTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGTTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  2:48987/142‑1 (MQ=255)
                                                                                                                                                   aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa    >  1:84496/3‑149 (MQ=255)
                                                                                                                                                    aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa   >  1:758257/2‑149 (MQ=255)
                                                                                                                                                     aaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  >  1:577693/1‑149 (MQ=255)
                                                                                                                                                   |||                                                                                                                                                     
TTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543876‑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: