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% 130.2 / NA 43G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/19);  total (24/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCAC  >  NC_000913_3_hsa_tpiA/1197356‑1197643
                                                                                                                                                  |                                                                                                                                             
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGTCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                             <  2:92507/149‑1 (MQ=255)
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                             <  1:678314/149‑1 (MQ=255)
aCTGATTAATTAATCAGAATATCTAATTTTTCCGTAACCTAATTCAAGAGCTTGTCAATGTTTCGCAAAGCTGCTAATATACTCGCAAAGATATCTAGAGGTCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                             <  2:609119/149‑1 (MQ=255)
  tGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtgt                                                                                                                                           >  1:157375/1‑148 (MQ=255)
          tcatcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                   <  2:726479/149‑1 (MQ=255)
              caGACTATCTAAGTTTTCCGCAGACAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATAATCGCAAAGATATCTAGAGGGAATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTc                                                                                                                               <  1:198155/149‑1 (MQ=255)
                 gatatCGAAGTTTTCCGCAGATGTAGTAAATGACTTGGCCATTTTTTCCAATGCTTCTAATATTCTCGCAAAGCTATCCAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGACCGAATATTAACTCTGTGCATACTATAGGTCaaa                                                                                                                            <  2:663388/148‑1 (MQ=255)
                           ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                  >  2:95642/1‑149 (MQ=255)
                           ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                  >  1:512467/1‑149 (MQ=255)
                             ttCCGCAGCCCAAATTAAGGGCGTGGCCATTTTTTTCAATGCTTCTAATATACTCGAAAAGATATCTAGAGGGCATCGTATGAGAATTTCCTTCAGTCTCTGAGAGAATATGAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                                <  2:335030/149‑1 (MQ=255)
                                    gCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTaa                                                                                                         <  1:353512/149‑1 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                        <  2:674938/149‑1 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                        >  1:626138/1‑149 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                        <  1:542521/149‑1 (MQ=255)
                                       aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCCTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                                                                                                      >  1:136167/1‑149 (MQ=255)
                                         aaTCAAGGGCGTGGCCATTTTTCTCAATGCTTATAATATACTCGCAAAGATATCTAGAGGACATCGTATGAGGTTTTACTCCATTCTCTGAGCGAATATTAACTCTGTGCACTTTATCGTTCAAAGCAAGTGCCGTTCAATTAATCttt                                                                                                    >  2:276507/1‑148 (MQ=255)
                                           tCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATCTTAACTCTATGCATATTATAGTTCAAAGCAAGTACCATACAATTCATGTTAtt                                                                                                  >  2:314040/1‑149 (MQ=255)
                                                ggCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                             >  1:7457/1‑149 (MQ=255)
                                                     ggCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTaaaa                                                                                        >  2:558599/1‑149 (MQ=255)
                                                        cATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                                                                                                           >  1:129301/1‑111 (MQ=255)
                                                        cATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                                                                                                           <  2:129301/111‑1 (MQ=255)
                                                        cATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAt                                                                                     >  1:186782/1‑149 (MQ=255)
                                                                       ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCg                                                                      >  1:180450/1‑149 (MQ=255)
                                                                                 cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttt                                                            >  1:446086/1‑149 (MQ=255)
                                                                                 cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGAATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttt                                                            <  1:239130/149‑1 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                         >  1:291949/1‑149 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                         >  1:666390/1‑149 (MQ=255)
                                                                                              cgagagGGCATCTTATGAGGTTTTGCTGCATTCTCTTAGCGAATATTAACTCTGTGCATATTATATTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgact                                               <  1:37951/147‑1 (MQ=255)
                                                                                               tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgactt                                              >  1:694061/1‑149 (MQ=255)
                                                                                                        tCGTGTGGGGTTTTGCTTCATTTGTTGGTCTAATATTCATTCTGTGCATATTATGGTTCAAAGCAATTACCATTCAATTAATGTTATTTCTAGTGAAAAATTCTTTGATTTTATCGTCCGGGCTGTTTCCGTCTCGAATTTCAGTTCTg                                     <  2:235923/149‑1 (MQ=255)
                                                                                                          gTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTggg                                   <  1:95642/149‑1 (MQ=255)
                                                                                                              gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                               >  2:347834/1‑149 (MQ=255)
                                                                                                               aGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                              <  1:569015/149‑1 (MQ=255)
                                                                                                               aGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCATCCGGGATATTTCCGACTTAGCTTTCAGTTTTGTGCTTaa                              >  2:409186/1‑149 (MQ=255)
                                                                                                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAAc                             >  2:422226/1‑149 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                            <  1:180716/149‑1 (MQ=255)
                                                                                                                         gcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGCCTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                    >  1:23816/1‑149 (MQ=255)
                                                                                                                            ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTTTGGGGTCCAGTGTtaaagcaaa                 >  1:275963/1‑149 (MQ=255)
                                                                                                                               tctGAGCGAATAGTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat              <  2:439002/149‑1 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAATataa            >  2:400794/1‑149 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa            >  1:425200/1‑149 (MQ=255)
                                                                                                                                  gAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTGATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAACataat           >  2:658210/1‑149 (MQ=255)
                                                                                                                                        atatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAACATAATACcact     >  2:203159/1‑149 (MQ=255)
                                                                                                                                           ttAACTCTGTTCATATTAAAGTTAAAAGCAATTACCATTCATCTACTCTTATTCTTACTGAAATATTCTTATATTTGCTCGTCCGCGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAACGATAATACcactcac  <  2:362977/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCAC  >  NC_000913_3_hsa_tpiA/1197356‑1197643

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