Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,197,502 C→T 100% G11R (GGA→AGA)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,197,5020CT100.0% 72.7 / NA 25G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/12);  total (13/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACC  >  NC_000913_3_hsa_tpiA/1197358‑1197637
                                                                                                                                                |                                                                                                                                       
tGAATAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAAGGCGTGACCCTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGCTTTTGCTGCATTCTCTAAGCGAACATTAACTCTGTGc                                                                                                                                     >  1:112831/1‑149 (MQ=255)
      aTTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGACGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTATGAGCGAATATTAACTATGTGCAtatta                                                                                                                               >  2:425201/1‑149 (MQ=255)
       ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTACCTCTGTGCAtattat                                                                                                                              <  1:322860/149‑1 (MQ=255)
               aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                      <  2:577846/149‑1 (MQ=255)
                   tCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGAAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTTTAGTTCAAAGCaa                                                                                                                  >  1:125840/1‑149 (MQ=255)
                     tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATAGTATGAGGTTTTGCTGCATTCTCTTAGCGAATATTAACTCTGTGCATATTTTTGTTCAAAGCACGt                                                                                                                >  2:402136/1‑149 (MQ=255)
                         ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                            >  1:329421/1‑149 (MQ=255)
                                                          tttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt                                                                           <  1:517198/149‑1 (MQ=255)
                                                              ttaaTGCGTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGGGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                       <  1:282029/147‑1 (MQ=255)
                                                              tcAATGCTTCTTATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGAATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATGTTTAGTGAAAAATTCttttat                                                                       <  2:582427/149‑1 (MQ=255)
                                                                             taCTCGAAAAGATATCGAGAGGGCATCGTATGAGGTTTTGCTGTATTCTCTGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAt                                                        <  2:534398/149‑1 (MQ=255)
                                                                                tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                     >  1:474255/1‑149 (MQ=255)
                                                                                             tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGGAAGTACCATTCAATTAATGTTATTTTTAGTAAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCTgccttgactt                                        >  2:151490/1‑149 (MQ=255)
                                                                                             tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAGCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgactt                                        >  1:609229/1‑149 (MQ=255)
                                                                                                agGGCATGGTATTAGGTTTGGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAAGTAATGTTATTTTTGGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCa                                     <  2:271246/149‑1 (MQ=255)
                                                                                                agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTTTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCa                                     <  2:532968/149‑1 (MQ=255)
                                                                                                     aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCt                                <  1:308298/149‑1 (MQ=255)
                                                                                                      tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                               <  2:251527/149‑1 (MQ=255)
                                                                                                                     ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                >  1:391975/1‑149 (MQ=255)
                                                                                                                      tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaagg               >  1:278015/1‑147 (MQ=255)
                                                                                                                             tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTACTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTGTTGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCTAAGat        >  1:498947/1‑149 (MQ=255)
                                                                                                                               ttaGCGAATATTAACTCTGTGCATATTATATTTCAAATCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa      <  1:34173/147‑1 (MQ=255)
                                                                                                                               tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa      <  1:445467/149‑1 (MQ=255)
                                                                                                                               tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa      >  1:246594/1‑149 (MQ=255)
                                                                                                                                   cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAcc  >  1:474356/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                       
TGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACC  >  NC_000913_3_hsa_tpiA/1197358‑1197637

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