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.A78.8% 152.3 / 36.6 52coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (6/5);  new base A (20/21);  total (26/26)
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
*NC_000913_3_hsa_tpiA1,544,0220CA98.1% 182.7 / ‑3.3 52L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (26/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543878‑1544170
                                                                                                                                                 |||                                                                                                                                                    
aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaaaa                                                                                                                                                  <  1:539347/150‑2 (MQ=255)
    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATa                                                                                                                                               >  4:471623/1‑150 (MQ=255)
    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaaaagata                                                                                                                                              <  3:199534/150‑6 (MQ=255)
      aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAAc                                                                                                                                             >  2:385203/1‑150 (MQ=255)
                    cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCa                                                                                                                               >  1:376592/1‑150 (MQ=255)
                       tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAttt                                                                                                                            >  4:186667/1‑150 (MQ=255)
                       tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAttt                                                                                                                            <  2:376592/150‑1 (MQ=255)
                           aaGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  4:255900/149‑1 (MQ=255)
                               aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCCGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                    <  1:403524/150‑1 (MQ=255)
                                       acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                            >  3:394255/1‑150 (MQ=255)
                                         aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          >  4:217525/1‑150 (MQ=255)
                                         aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                         <  3:230000/150‑1 (MQ=255)
                                             ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                      <  1:532766/150‑1 (MQ=255)
                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTCTTTTTCTCATTTAGCAGTGCGATATTGTTCCTcac                                                                                                     >  1:277051/1‑150 (MQ=255)
                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaaaa                                                                                                                                                  >  1:606049/1‑103 (MQ=255)
                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaaaa                                                                                                                                                  <  2:606049/104‑2 (MQ=255)
                                               cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    <  2:371262/150‑1 (MQ=255)
                                               cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    >  4:384338/1‑150 (MQ=255)
                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                        >  2:744306/1‑150 (MQ=255)
                                                            actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGACGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                       <  3:299082/150‑1 (MQ=255)
                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                                <  2:277533/139‑1 (MQ=255)
                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                                >  1:277533/1‑139 (MQ=255)
                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     >  4:11271/1‑150 (MQ=255)
                                                                  tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                 >  1:525876/1‑150 (MQ=255)
                                                                   gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                   <  1:189947/115‑1 (MQ=255)
                                                                   gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                   >  2:189947/1‑115 (MQ=255)
                                                                   gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa                                                                                <  2:525876/150‑1 (MQ=255)
                                                                     tatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaa                                                                              <  1:744306/150‑1 (MQ=255)
                                                                       tCCAGAGATTGCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGGCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  4:307276/150‑1 (MQ=255)
                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCCCACATTGGTAATGAAAAAAAGa                                                                          >  4:341288/1‑150 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                <  3:337160/111‑1 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                         >  1:24353/1‑150 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                >  4:337160/1‑111 (MQ=255)
                                                                                                ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAGTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AAt                                                  <  3:384338/150‑1 (MQ=255)
                                                                                                  ttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAg                                                >  1:698572/1‑150 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                           >  1:638140/1‑88 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                           <  2:638140/88‑1 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAg                                            >  1:412351/1‑150 (MQ=255)
                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGaa                                   >  2:653038/1‑150 (MQ=255)
                                                                                                                        cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGc                                               >  2:733296/1‑129 (MQ=255)
                                                                                                                        cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGc                                               <  1:733296/129‑1 (MQ=255)
                                                                                                                           aTTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAAt                       <  1:634804/150‑1 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  3:11271/150‑1 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  2:412351/150‑1 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:671507/150‑1 (MQ=255)
                                                                                                                               tCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcgc                                                      >  1:548277/1‑116 (MQ=255)
                                                                                                                               tCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  2:698572/150‑1 (MQ=255)
                                                                                                                               tCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATCCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcgc                                                      <  2:548277/116‑1 (MQ=255)
                                                                                                                                 tGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAAAAGCAAGTAATTAAAATTAAGGTTGAATTgggg                  >  1:193471/1‑149 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                          >  4:170535/1‑44 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                          <  3:170535/44‑1 (MQ=255)
                                                                                                                                        cAGTCGCG‑AAAAAAGATAACGACTATTTTTTTTTTTTTTTAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  4:331332/1‑150 (MQ=255)
                                                                                                                                            tcgcgaaaaaaGATAACGACTATTTTTTTTATTTTTCTTTTCCATATTGTTCCTCACATGCACACTTTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGAATGAATTTGGTGAAGTGCTGAt       >  2:78139/8‑150 (MQ=255)
                                                                                                                                            tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       >  3:202160/8‑150 (MQ=255)
                                                                                                                                            tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       >  4:490305/8‑150 (MQ=255)
                                                                                                                                                 aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtagcaa  >  1:196/3‑145 (MQ=255)
                                                                                                                                                 aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  3:217525/148‑1 (MQ=255)
                                                                                                                                                  aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     <  4:220567/66‑1 (MQ=255)
                                                                                                                                                  aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     >  3:220567/2‑67 (MQ=255)
                                                                                                                                                 |||                                                                                                                                                    
ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGC‑AATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543878‑1544170

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