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,0980GT100.0% 103.1 / NA 35L50L (CTC→CTA) ykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/21);  total (14/21)

TCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263957‑264243
                                                                                                                                             |                                                                                                                                                  
tCCCAGGCGCAGTGTCATGTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTATCCACTACTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                                                                                                             <  2:520109/149‑1 (MQ=255)
    aggccgaGTGAAATCACCCGGTGCAGGCCCCTGGCAGGGTTGCACATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                         <  2:157167/143‑1 (MQ=255)
     ggCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCg                                                                                                                                        >  1:67424/1‑149 (MQ=255)
          aGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                   <  1:4546/149‑1 (MQ=255)
           gtgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                                  <  1:574653/149‑1 (MQ=255)
              tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGa                                                                                                                               >  1:626921/1‑149 (MQ=255)
               aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                              <  2:497777/149‑1 (MQ=255)
                    cccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                                         <  2:67424/149‑1 (MQ=255)
                           cAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTCGGACTCCGGCGATATCCGGAACGTCCCTGCCgg                                                                                                                  >  1:212277/1‑149 (MQ=255)
                              gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCCTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGCCCATATCCGGAACGTCCCTGCCGgcgt                                                                                                               >  2:328037/1‑149 (MQ=255)
                                              cTCATTACAGTTCTCCACTATTGACAGTTCAGGGTGATGCTCGTCGGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTAcc                                                                                              >  1:539041/1‑149 (MQ=255)
                                                   tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 <  2:255003/126‑1 (MQ=255)
                                                   tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 >  1:255003/1‑126 (MQ=255)
                                                        ttCTCGACTTTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCAATTGCGGTCCTGGCGGCGTTCGGCGTATCCCGCTCGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTgg                                                                                    <  1:131496/149‑1 (MQ=21)
                                                           tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTAt                                                                                 >  1:510204/1‑149 (MQ=255)
                                                                 aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCg                                                                           <  1:380218/149‑1 (MQ=255)
                                                                     tcaggtcaggGTGATGCCCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGc                                                                       <  1:67205/149‑1 (MQ=255)
                                                                         ttgaGGGTCATGCTCATCATGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                   <  2:335292/146‑1 (MQ=255)
                                                                         ttCCGGGGGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                   <  2:335300/149‑1 (MQ=255)
                                                                               gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                             <  1:251643/149‑1 (MQ=255)
                                                                                    gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                        <  2:611087/149‑1 (MQ=255)
                                                                                      tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGTCGGCGGTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                      <  2:26196/149‑1 (MQ=255)
                                                                                             gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                         >  2:417144/1‑92 (MQ=255)
                                                                                             gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                         <  1:417144/92‑1 (MQ=255)
                                                                                             gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCCGTGTAGAGACAGGCG‑CTACCGTACCACTAGTAGTTGCCGCTGGCGCAGCGAAAAGTTCCGGACTccg                                               >  2:104156/1‑147 (MQ=255)
                                                                                               aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTCAGGACAGATGTATTTTTTATGGCT‑GTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                             <  1:501789/149‑1 (MQ=18)
                                                                                                  tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGGTGc                                          >  2:395351/1‑149 (MQ=255)
                                                                                                       ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 <  2:375125/74‑1 (MQ=255)
                                                                                                       ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 >  1:375125/1‑74 (MQ=255)
                                                                                                         cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                       >  2:83472/1‑129 (MQ=255)
                                                                                                         cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                       <  1:83472/129‑1 (MQ=255)
                                                                                                                     cTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCGGGATGCTGGCGGTATATCGTCATTg                       >  1:6705/1‑149 (MQ=255)
                                                                                                                              ttCGGCTTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt              <  2:510204/149‑1 (MQ=255)
                                                                                                                                cggcggAGACCGCTAGGACTGCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTa            <  2:626921/149‑1 (MQ=255)
                                                                                                                                  gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACt          >  1:491821/1‑149 (MQ=255)
                                                                                                                                          cgcgAGGACGCCGGCGATAGCGGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAttt  <  1:348485/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACG‑CTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263957‑264243

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