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,5020CT98.2% 174.0 / NA 55G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (30/23);  minor base G (0/1);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACT  >  NC_000913_3_hsa_tpiA/1197356‑1197640
                                                                                                                                                  |                                                                                                                                          
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGTGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                          <  1:404643/149‑1 (MQ=255)
    aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACGCTGTGCat                                                                                                                                      >  1:38523/1‑149 (MQ=255)
    aTTAATTCATCAGAACATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                      >  1:401872/1‑149 (MQ=255)
        aTTCATCAGTATAACTAACTCTTCCGCAGCCAAAGTCAAGGGCGTGGCCATTTTTGTTAATTTTTTGAATATACTCGCAAAGATATCTAGAGGGCATCGTATAGAGCTATGATGCATTCTCTGAGCGAATATTAACAGTGTGCAtatta                                                                                                                                  >  2:717188/1‑149 (MQ=255)
             tcaGAATATCTAAGTTTTCCGCAGCCAAAATCACGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAAATTTAACCCCGTGCAtgttctgcat                                                                                                                             >  2:638029/1‑143 (MQ=255)
             tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                             >  1:122887/1‑149 (MQ=255)
                gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaa                                                                                                                          >  2:807370/1‑149 (MQ=255)
                   tatCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGc                                                                                                                       <  1:686413/149‑1 (MQ=255)
                      cTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAg                                                                                                                    >  1:311823/1‑149 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                         >  2:466875/1‑143 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                         <  1:466875/143‑1 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                                                                                                   >  1:11088/1‑149 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAAGATGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                                                                                                   >  2:104874/1‑149 (MQ=255)
                           ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                               <  1:22833/149‑1 (MQ=255)
                            tttCCCCAGCCAAGATCATGGGTGTGGCCATTTTTCTCGATGCTTCTAATATACTCGCAAAGATTTCTAGAGGGTATCGTATTTGGTTTTGCTGTATTCTCTTAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAt                                                                                                              <  2:398870/149‑1 (MQ=255)
                             ttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                             <  2:779257/149‑1 (MQ=255)
                                 gCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAAt                                                                                                         <  1:255582/149‑1 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                     <  1:176244/149‑1 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                     <  2:804266/149‑1 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAACCAAGTACCATTCAATTAAt                                                                                                     <  2:765317/149‑1 (MQ=255)
                                       aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCAGTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                                                                                                   >  2:838376/1‑149 (MQ=255)
                                               gggCGTGGCTATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                                                                                           <  2:900862/149‑1 (MQ=255)
                                                   gTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGa                                                                                       <  1:29994/149‑1 (MQ=255)
                                                           ttttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTACTGCATGCGCTGAGCGAATATTAACTCTGTGCATATTATAGTTCTGAGTAAGTGGCATTCGATTAGTGTTATTTTTGGTGAAAAAttgt                                                                               >  1:79089/1‑147 (MQ=255)
                                                               ctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATTTTATTTTTAGTGAAAAATTCtttta                                                                           >  1:179600/1‑149 (MQ=255)
                                                                tcAGTGCTTCGAATATACTCGGATAGATAGCTGGAGGGCATCGTATGAGGATTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAACATTCttttat                                                                          >  2:348573/1‑149 (MQ=255)
                                                                       ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCg                                                                   >  1:466785/1‑149 (MQ=255)
                                                                       ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTATTTTATTTTATCg                                                                   >  1:467034/1‑149 (MQ=255)
                                                                          tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGGGGTTTTTCTGCATACAATGAGAGAATATTAAATCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATAAATGTTATTTTTAGTGCAAAATTATTTTATTTTATCGtac                                                                >  2:123420/1‑147 (MQ=255)
                                                                          tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGTATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAATTATAATGAAAAATTCTTTTATTTTATCGTcc                                                                <  2:802291/149‑1 (MQ=255)
                                                                           aaTATACTCGCAAAGATATCTAGAGTGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTAAAATCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                               <  1:402404/149‑1 (MQ=255)
                                                                           aaTATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                               <  2:541321/149‑1 (MQ=255)
                                                                             tataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                >  2:67693/1‑146 (MQ=255)
                                                                             tataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                <  1:67693/146‑1 (MQ=255)
                                                                              ataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGa                                                            >  2:571747/1‑149 (MQ=255)
                                                                                  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                        >  1:555275/1‑149 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCCTATTATATTTCAAAGCAAGTAAAATTCAAATAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                      >  2:683452/1‑149 (MQ=255)
                                                                                     cAAAGATATCTAGAGGGCAACGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTcc                                                     <  1:131790/149‑1 (MQ=255)
                                                                                            atCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttga                                              >  2:338141/1‑149 (MQ=255)
                                                                                                  agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCa                                        <  2:401872/149‑1 (MQ=255)
                                                                                                    ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGt                                      >  1:426980/1‑149 (MQ=255)
                                                                                                       aTCGTATGAGGTTTTGCTGCATTCTATGAGCGAATATGGACTCTGTGCATATTATAGTTCAAAGCAAGTGGCATTCAATTAATGTTATTTTAAGTGCAAAGTTCTTTTATTTATTGGTAAGGTATTTTTGTACCTTGACtttagggtct                                   >  2:189180/1‑142 (MQ=255)
                                                                                                         cgaaTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTAGGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                                 <  1:388351/146‑1 (MQ=255)
                                                                                                          gTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGGTTTTTCCGCCTTGACTTTAAATTCTggc                                >  1:119306/1‑148 (MQ=255)
                                                                                                           tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGc                               >  2:239661/1‑149 (MQ=255)
                                                                                                              gAGGTTTTTGTGGATTCTCTGAGCGAATATTAATTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                            <  2:121539/149‑1 (MQ=255)
                                                                                                                      gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGTGATTTTTCTGACTTGACTTTCAGTTCTGGTCTTACCTGtaac                    >  2:36167/1‑146 (MQ=255)
                                                                                                                          cATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagca                >  1:676080/1‑149 (MQ=255)
                                                                                                                            ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAATACCATAAAATTAATGTTATTCTTAGTGAACAATTCTCCTATGTTATCGTCCTGTATTTTTCCGACTTGCATTTCATTTCTGGCCTTAACTGTcatagcaaa              >  2:99884/1‑149 (MQ=255)
                                                                                                                                ctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGata          <  2:676080/149‑1 (MQ=255)
                                                                                                                                     cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTAAACTGTCAAAGCAAAGATAATAcc     >  1:184473/1‑149 (MQ=255)
                                                                                                                                      gAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAg                                       <  2:165542/114‑1 (MQ=255)
                                                                                                                                      gAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAg                                       >  1:165542/1‑114 (MQ=255)
                                                                                                                                        atatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcact  <  1:571747/149‑1 (MQ=255)
                                                                                                                                        atatTAACTCTGTGCATAGTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcact  <  2:384739/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACT  >  NC_000913_3_hsa_tpiA/1197356‑1197640

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