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,5020CT97.3% 120.8 / NA 37G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/17);  minor base A (1/0);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA  >  NC_000913_3_hsa_tpiA/1197356‑1197642
                                                                                                                                                  |                                                                                                                                            
aCTGATTAGTTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                            <  1:602356/149‑1 (MQ=255)
 cTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                           <  1:473914/149‑1 (MQ=255)
        atgcatcaGAATATCTAAGTTTTCAGCAGCCAAAATCAAGGGCGTGGCCATGTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATCCTCCGAGCGAATATTAACTCTGTGCAtatta                                                                                                                                    <  2:392888/146‑1 (MQ=255)
          tcatcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTAAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                  <  2:275133/149‑1 (MQ=255)
          tcagcaGAATAACTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTGAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                  <  2:167435/149‑1 (MQ=255)
            atcgGAATCACTAAGTTGGCCTAAGACAAAATCAAGGGCGTGGCTCTTTTTCACAGTGCTTCTAATATCCTCGCACAGCTATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTTTCTGAGCGAATATTAACTTTGTGCATATTATAGt                                                                                                                                <  2:363681/149‑1 (MQ=255)
            atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                                >  1:569380/1‑149 (MQ=255)
                 aaTATCTAAGTCTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGTGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                           <  1:483157/149‑1 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                                                                                                     >  1:127852/1‑149 (MQ=255)
                                       aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                                                                                                     >  1:339917/1‑149 (MQ=255)
                                       aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCAGTGCATCTTATAGTTCAAAGCAAGTACCATTCAATTAATGt                                                                                                     >  2:704984/1‑149 (MQ=255)
                                        aaaTCGAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGtt                                                                                                    >  1:472449/1‑149 (MQ=255)
                                          aTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAt                                                                                                  >  2:604410/1‑149 (MQ=255)
                                             aaGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCCTATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtttt                                                                                               <  2:127852/149‑1 (MQ=255)
                                                 gCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGt                                                                                           >  2:98292/1‑149 (MQ=255)
                                                                 cAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTCTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatt                                                                           <  2:472449/149‑1 (MQ=255)
                                                                  aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttattt                                                                          >  1:227483/1‑149 (MQ=255)
                                                                  aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttattt                                                                          >  1:252437/1‑149 (MQ=255)
                                                                    tGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                            <  1:24194/145‑1 (MQ=255)
                                                                    tGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                            >  2:24194/1‑145 (MQ=255)
                                                                       ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCg                                                                     >  2:256336/1‑149 (MQ=255)
                                                                                        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                    <  1:616509/149‑1 (MQ=255)
                                                                                        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                    <  2:227483/149‑1 (MQ=255)
                                                                                            atCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATGCTTTTATGTTATCGTCCGGGATTTTTCCgacttga                                                <  2:339917/149‑1 (MQ=255)
                                                                                              cTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATCTTATTTATCGGTATTTTTCCgacttgact                                              <  2:577616/149‑1 (MQ=255)
                                                                                               tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgactt                                             >  1:643544/1‑149 (MQ=255)
                                                                                                agagGGCATCGTATGAGTTGTTGCTGCATTCTCCGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                            <  1:373588/149‑1 (MQ=255)
                                                                                                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGNCCGGNATTTTTCCGACTTGGCTTTCGGTTCTg                                    >  1:829/1‑149 (MQ=255)
                                                                                                             tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                                                                                                     >  1:671917/1‑63 (MQ=255)
                                                                                                             tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                                                                                                     <  2:671917/63‑1 (MQ=255)
                                                                                                              gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTTACTTTTAGTTCTGGGCTTa                              >  2:23741/1‑149 (MQ=255)
                                                                                                                  ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTg                          >  1:84520/1‑149 (MQ=255)
                                                                                                                       ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATGGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                     >  2:460091/1‑149 (MQ=255)
                                                                                                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                    >  2:28843/1‑149 (MQ=255)
                                                                                                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                    >  1:224899/1‑149 (MQ=255)
                                                                                                                            ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa                <  2:252437/149‑1 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGGCTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat             <  1:468899/149‑1 (MQ=255)
                                                                                                                                          atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca  >  1:362480/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA  >  NC_000913_3_hsa_tpiA/1197356‑1197642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑