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

ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTCAAAGCAAAGATAATACCA  >  NC_000913_3_hsa_tpiA/1197360‑1197638
                                                                                                                                              |                                                                                                                                         
aTTAATGCATCAGAACACCCAAATCATCCGCAACCTAAATGAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCTCAAAGATATCGAGAGGGCATCGTATGAGGTTTTGCTGCATTCTGTGTGCTAATATTAACTCTGTGCata                                                                                                                                    <  2:563730/150‑1 (MQ=255)
    aTTGATCAGAATATGTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                                                                                                                                <  2:89345/150‑1 (MQ=255)
        atcaGACTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTGAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAAGTTTTGCTGCATTCGCTGAGCGCATATTAACTCTGTGCATATTATCGtt                                                                                                                            >  2:40090/1‑150 (MQ=255)
        atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATTCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                            <  2:329828/150‑1 (MQ=255)
            gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                        >  2:480233/1‑150 (MQ=255)
             aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                                                                                                       <  2:30631/150‑1 (MQ=255)
             aaTATCTAAGGTTTCCGCAGTCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGTGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                                                                                                       <  2:101369/150‑1 (MQ=255)
                atCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                    <  1:175701/150‑1 (MQ=255)
                 tCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAg                                                                                                                   >  2:503512/1‑150 (MQ=255)
                    aaGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                             <  2:507052/137‑1 (MQ=255)
                    aaGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                             >  1:507052/1‑137 (MQ=255)
                          tCCGTAGTCAAAATAAAGGGCGTGGCCATTTTTCTCAAGGGTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                                                                                                          <  2:318043/150‑1 (MQ=255)
                             gCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCAGTGCATCTTATAATACAAACCAAGTACCATTCaagt                                                                                                       >  2:656664/1‑148 (MQ=255)
                                gCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                    <  2:118844/150‑1 (MQ=255)
                                          aGGGCGGGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                                                                                          <  1:320604/150‑1 (MQ=255)
                                            ggCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTAAAAGCAAGTACCATTCAGTTAATGTTCTTTTTAGt                                                                                        >  1:198090/1‑150 (MQ=255)
                                                            tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATACTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatt                                                                        >  2:600187/1‑150 (MQ=255)
                                                                            aCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttt                                                        >  1:480763/1‑150 (MQ=255)
                                                                                     aatatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgactt                                               <  2:275225/149‑1 (MQ=255)
                                                                                       tatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttga                                             <  1:15290/150‑1 (MQ=255)
                                                                                          cTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATCTTTCCgccttgactt                                          >  1:629456/1‑150 (MQ=255)
                                                                                                  cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCt                                  <  2:480763/150‑1 (MQ=255)
                                                                                                   agcGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑g                                <  2:685520/148‑1 (MQ=255)
                                                                                                    tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑gg                               <  2:501774/150‑1 (MQ=255)
                                                                                                     cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGGTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑ggg                              <  1:303847/150‑1 (MQ=255)
                                                                                                          gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgactt                                               <  1:228362/129‑1 (MQ=255)
                                                                                                          gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgactt                                               >  2:228362/1‑129 (MQ=255)
                                                                                                             gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACtt                      >  2:360278/1‑149 (MQ=255)
                                                                                                               tttACTCCATTCTCTGAGCGAATATTAAATCTGTGCATATTATAGTTCGAAGCCAGTATCACTCAAGTAATGTTGTTTTTAGTGAAAAAGTGTTTTATTCTATCCTCGGCAAAATTTCAGAATTAACTTTAATTTCTCGGCCTTAACTGt                     >  2:80672/1‑150 (MQ=255)
                                                                                                                  gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTcaaa                 <  1:480233/150‑1 (MQ=255)
                                                                                                                   ctgcATTCTCTGAGCGACTAATAACTCAGTTCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTcaaag                >  2:235713/1‑150 (MQ=255)
                                                                                                                    tgcATTCTTTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTcaaagc               >  1:179975/1‑150 (MQ=255)
                                                                                                                             tGAGCTAATATTAACTCTGAGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCTGGATTTTTCAGAGTTAAAGCCCGATACT‑GGGCTTAACTGTTAAAGCAAAGataat      >  2:317133/1‑150 (MQ=255)
                                                                                                                               aGCGAATATTAACTCGGTGGATATTATGGTGCGAATGAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCGGACTTGACTTTCAGTTCT‑GGGCTTAACTGTCAAAGCAAAGATAacac    <  2:329734/150‑1 (MQ=255)
                                                                                                                                 cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTCAAAGCAAAGATAATACca  <  1:682410/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                         
ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCT‑GGGCTTAACTGTCAAAGCAAAGATAATACCA  >  NC_000913_3_hsa_tpiA/1197360‑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: ‑