Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 264,098 G→T L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA264,0980GT97.1% 100.5 / NA 35L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (18/16);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_hsa_tpiA/263952‑264242
                                                                                                                                                  |                                                                                                                                                
cGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                                <  2:350409/149‑1 (MQ=255)
               aGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                 >  1:773717/1‑149 (MQ=255)
                                catgCCCTTGGTAGGGGTGCTCATTATAGTGCTTCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGAGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGGCCCTGCCgg                                                                                                                <  2:802361/146‑1 (MQ=25)
                                  ggCCCCTGGTAGGGGTGCTCATTACAGTTCTCCACTATTGTTAGTTTAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                                                                                                              <  2:674591/149‑1 (MQ=255)
                                             gggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                   <  1:224033/149‑1 (MQ=255)
                                                      aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg                                                                                          >  2:633846/1‑149 (MQ=255)
                                                          cAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCAc                                                                                      <  1:418733/149‑1 (MQ=255)
                                                                tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATATGGAACGTCGCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                                >  1:273758/1‑149 (MQ=255)
                                                                   aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACAGTGCCACTGGTATTTg                                                                             >  1:21731/1‑149 (MQ=255)
                                                                    cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGCCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCATGCCACTGGAATTTGc                                                                            >  2:338468/1‑149 (MQ=255)
                                                                      tctgtgAGTTCGGGGGGATGCTTATATGGGGGGTAGAGGGCCCATGGGGGTGCTGGCGGAGTTCGGCGTATACCGCGAGGACTCGGGCGATATCCGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATGTGCCg                                                                          <  2:201019/147‑1 (MQ=18)
                                                                         gTCAGTTCAGGGGGATGCTCATCCGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAAAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                                       <  1:215516/149‑1 (MQ=255)
                                                                                                aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCGGTAGCGAAAGATTCGGGGCTc                                                >  1:224969/1‑149 (MQ=255)
                                                                                                 ggcaggTATAGGGCCAATTGCGGTCCTGGCGGCGTTCGGCGTATACAGCTAGGACTCCGGCGATATCCGCAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCCGTTCAGAAGCGAAACATTCGGGACTCa                                               >  2:679817/1‑149 (MQ=18)
                                                                                                       tataGTGCCCCTTGCGGCCCCGGCGGCGTTCTGCGTATACCGCTAGGACTCCGGTGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                         <  1:713321/149‑1 (MQ=17)
                                                                                                       gataGGGCCCAGTGGCGTCCTGGGGGCGTCCGGCGTATGCCGCTAGGACTCCTGTGATATACGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                         <  1:713329/148‑1 (MQ=14)
                                                                                                          aGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTgg                                      <  1:66017/149‑1 (MQ=255)
                                                                                                            ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                          <  1:529351/79‑1 (MQ=255)
                                                                                                            ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                          >  2:529351/1‑79 (MQ=255)
                                                                                                              cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGt                                  >  2:156579/1‑149 (MQ=255)
                                                                                                                 aTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtata                               >  2:819140/1‑149 (MQ=255)
                                                                                                                  ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAtt                                                                               >  2:171022/1‑100 (MQ=255)
                                                                                                                  ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAtt                                                                               <  1:171022/100‑1 (MQ=255)
                                                                                                                    gCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCCGGATGATAGCGGTATATCg                            >  2:806474/1‑149 (MQ=255)
                                                                                                                      ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTc                          <  1:300019/149‑1 (MQ=255)
                                                                                                                          cTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTg                      >  2:767271/1‑149 (MQ=255)
                                                                                                                                   ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt             >  1:25733/1‑149 (MQ=255)
                                                                                                                                     cGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCGGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTa           <  2:531194/149‑1 (MQ=255)
                                                                                                                                      gccttacaCCGCTGGGACTCCGGCGATTTCCGGAACGTCCCTTCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc          <  1:49302/142‑1 (MQ=14)
                                                                                                                                          tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGccc                   >  2:396029/1‑136 (MQ=255)
                                                                                                                                          tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGccc                   <  1:396029/136‑1 (MQ=255)
                                                                                                                                             aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                               <  1:734731/105‑1 (MQ=255)
                                                                                                                                             aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                               >  2:734731/1‑105 (MQ=255)
                                                                                                                                             aCCGCTAGGACTCCGGCGATACCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCCGTCGCGAAAGATGCGGGACTCAGGATGCTGGCGGCATACCGTACTTGCCGTGCGTCTTCTGATCAt   >  2:124322/1‑149 (MQ=14)
                                                                                                                                              ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAtt  >  2:458362/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_hsa_tpiA/263952‑264242

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