Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,197,502 C→T 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% 126.2 / NA 39G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

TCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCA  >  NC_000913_3_hsa_tpiA/1197369‑1197638
                                                                                                                                     |                                                                                                                                        
tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                           >  1:624440/1‑149 (MQ=255)
    aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                       <  2:130010/149‑1 (MQ=255)
            aGTTTTCCGCAGGCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAc                                                                                                               <  2:849440/149‑1 (MQ=255)
              ttttCCGCAGCCGAAATCAAGAGCGTGGTCATTTTTCTCAATGGTTCTAATATATTCGCAAAGATAGCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTActa                                                                                                             <  2:150745/149‑3 (MQ=255)
                 tCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAAAATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTc                                                                                                          <  1:595496/149‑1 (MQ=255)
                       gCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTaa                                                                                                    <  2:282488/149‑1 (MQ=255)
                          cataTCAAGTGCGTGGCCATTTTTCTCAATGCTTCTAATATATTCTCAAAGATGTCTAGAAGGCCGCGTTTGAGGTTTTGCTGCATTTTGTTAGCGAATATTTCATTTGTCGCTATTTTAGTTAAAAACAAGTACCATTCACTTTATGt                                                                                                 >  2:680193/4‑149 (MQ=255)
                                  gggCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                                                                                         <  1:110746/149‑1 (MQ=255)
                                     cGTGGCCATTTTTCTCAATGCTTCTAATATACTAGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTg                                                                                      >  1:708164/1‑149 (MQ=255)
                                            aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAtt                                                                               >  1:213304/1‑149 (MQ=255)
                                                   ttcaTGCTTCTAATATAACATCAAAGATGTCTAGAGGCGATCGTATGAGGTTTTGCTGCATTCTCTGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGGACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                        <  2:153530/146‑1 (MQ=255)
                                                   ccAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGTGTTGGTTGCCTTCTCTGAGCGAATATTAACTCTGTGCATATTATCGTTCAAAGCAAGTACAATTCTATTAATGTTATTGTTGTTGAAAAATTTttttat                                                                        >  2:655862/2‑149 (MQ=255)
                                                     aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                           <  2:569362/112‑1 (MQ=255)
                                                     aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                           >  1:569362/1‑112 (MQ=255)
                                                             tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTa                                                                                                     <  1:191435/110‑1 (MQ=255)
                                                             tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTa                                                                                                     >  2:191435/1‑110 (MQ=255)
                                                                 ataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTCTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAACAATTCTTTTATTTTATAGTACGGGa                                                          >  2:738412/1‑149 (MQ=255)
                                                                   aCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCgtgttt                                                        >  2:231598/1‑144 (MQ=255)
                                                                       gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                    >  1:686645/1‑149 (MQ=255)
                                                                       gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                    <  2:213304/149‑1 (MQ=255)
                                                                         aaaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCg                                                  >  1:113229/1‑149 (MQ=255)
                                                                                 cTAGAGGGCATCGTATGGGGTTTTGCTGCATTCGCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCCAGACTTTTCCgactggact                                          <  2:323376/149‑1 (MQ=255)
                                                                                    gagGGCATCGTATGATGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTc                                       <  1:52057/149‑1 (MQ=255)
                                                                                      gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAg                                     >  2:340683/1‑149 (MQ=255)
                                                                                      gggCATCGTATGAGGTTTTGCTGCATTCTATGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCGGGGATTTTTTCGACTTGACTTTcgg                                     >  2:155277/1‑147 (MQ=255)
                                                                                         cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTc                                  >  1:681609/1‑149 (MQ=255)
                                                                                          aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTCGTGAACAGTTCTTTTATTTTGTCGTCAGGGATTTTTCCGACTTGACTTTCAGTTca                                 >  1:335478/1‑148 (MQ=255)
                                                                                          aTCGTATGAGGGTTTGATGCATTCTCTGAGCGAATAATAACTCTCTGCATATTATAGTTCAAAGCAAGAACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATAGTCCGGGATCTTTCCGACGTGCCTTTCAGTTca                                 >  2:539602/1‑148 (MQ=255)
                                                                                           tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                                <  1:340683/149‑1 (MQ=255)
                                                                                            cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                               <  1:805286/149‑1 (MQ=255)
                                                                                               aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                            >  1:44732/1‑149 (MQ=255)
                                                                                                 gAGGATTTGCTGCATTCTTTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTAAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTATGTTCTGGCCTTa                          >  2:836578/1‑149 (MQ=255)
                                                                                                     ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTg                      <  2:276017/149‑1 (MQ=255)
                                                                                                      tttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGt                     <  2:686645/149‑1 (MQ=255)
                                                                                                         gctgcATTCTTTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaa                  <  2:862258/149‑1 (MQ=255)
                                                                                                                  tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtt                                                        <  2:275028/102‑1 (MQ=255)
                                                                                                                  tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtt                                                        >  1:275028/1‑102 (MQ=255)
                                                                                                                  tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat         >  2:704454/1‑149 (MQ=255)
                                                                                                                         gAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGGTTTTTCCGACGTGACTTTCAGTTCTGGGCTTAACTGCGAAGGCAAAGATAGTAcac  >  1:70881/1‑147 (MQ=255)
                                                                                                                                     |                                                                                                                                        
TCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCA  >  NC_000913_3_hsa_tpiA/1197369‑1197638

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