Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,197,717 C→A 20.7% C174F (TGC→TTC)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,197,7170CA20.7% 79.1 / 4.5 41C174F (TGC→TTC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/6);  new base A (2/7);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 1.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAT  >  NC_000913_3_hsa_tpiA/1197576‑1197864
                                                                                                                                             |                                                                                                                                                   
tCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCaaaaaaa                                                                                                                                              >  2:405627/1‑149 (MQ=255)
  cGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGa                                                                                                                                            >  1:183644/1‑149 (MQ=255)
   gggATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGat                                                                                                                                           >  1:701390/1‑149 (MQ=255)
                      tttCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAAt                                                                                                                        >  1:634425/1‑149 (MQ=255)
                      tttCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATAAACATATGTTAGCAAAAAAAGATTATCAGAAAGGAAATTact                                                                                                                        >  2:16095/1‑147 (MQ=255)
                        tCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATGAAGACAGCAAAGAACCTTCAACATGAAAAAAATCCATTTGTTTCAAAAAAAAGATTATTAGCAAGGACAtaacact                                                                                                                      >  2:2531/1‑143 (MQ=255)
                          aGTCCTGGGCTAAACTGTCAAAGAAAAGATAATAACACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGATCCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCaa                                                                                                                    <  1:633633/149‑1 (MQ=255)
                              cTGGGATTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAACAGAATAACGCCCGTCAAGTACCTTCAACATGAAAAATATCCATTTGTTTGAAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTAt                                                                                                                <  2:626138/149‑1 (MQ=255)
                                           tCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCAAAAATTAaaaa                                                                                                   >  1:184252/1‑149 (MQ=255)
                                             aaagcaaagATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCaaaaaa                                                                                                 >  1:30699/1‑149 (MQ=255)
                                                         ataCCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                                                     >  1:271511/1‑149 (MQ=255)
                                                          taCCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                                                    <  1:254491/149‑1 (MQ=255)
                                                            ccACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGta                                                                                  >  1:596306/1‑149 (MQ=255)
                                                                        aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTAATTCCa                                                                      >  1:312019/1‑149 (MQ=255)
                                                                         aTGATAAAATTACAAAAAAGAATAAAGAAAGCAAAGAATCTTAAAAATGAAAAATATCCAATAAAAAAAAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                                                     <  2:13845/149‑1 (MQ=255)
                                                                               aaaTTACAAAAAAGAATAAAGACAGGAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTCATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGt                                                               <  2:666390/149‑1 (MQ=255)
                                                                                      aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGACGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTACTTCCAGAAGAGTTCAatat                                                        >  1:493657/1‑149 (MQ=255)
                                                                                      aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                                                        >  1:597307/1‑149 (MQ=255)
                                                                                       aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATAACCATATGAATGAAAACACACCAACTTCAGAAGGAAATTCATGCAACTGTCGTGACTCTAGAAAATGGGCAAAAAGATTGAACTTTATGCCAGAACAATTAAAtata                                                       >  2:151941/1‑149 (MQ=255)
                                                                                        aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATaa                                                      >  1:544651/1‑149 (MQ=255)
                                                                                           aGAAAAAAGACACGAAAGATCTCTCAACATGAAAAATATACAAAAAAAAAAAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                                                   <  2:192222/149‑1 (MQ=255)
                                                                                               aaaaaACAGCAAAGAACCTTAAACATGAAAAAGATCCAATAGTATGAAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTCTACTTCATACCAGAAGAGTTCAATATAATGTTTGt                                               <  2:522161/148‑1 (MQ=255)
                                                                                               aaaaaACAGAAAAGAACCTTTAACATGAAAAATATCAATTTGATAGAAAAAAAAGATTCTTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGCAGATACCCATAGAATCTTTGt                                               <  2:632979/148‑1 (MQ=255)
                                                                                                  aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCATAAGAGTTCAATATAATGTTTGTCtt                                            >  1:428551/1‑149 (MQ=255)
                                                                                                   gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAACTTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCAAGAAGAGTTCAATATAATATTTGTCTTc                                           >  2:52579/1‑149 (MQ=255)
                                                                                                   gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATCTAAAATTTGTCTTc                                           >  2:283316/1‑149 (MQ=255)
                                                                                                   gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc                                           <  2:271511/149‑1 (MQ=255)
                                                                                                    aCAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCa                                          >  1:198388/1‑149 (MQ=255)
                                                                                                    aCAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTGTATACCAGAAGAGTTCAATATAATGTTTGTCTTCa                                          >  1:198304/1‑149 (MQ=255)
                                                                                                      tgCAAAGAGGCATGAACAAGAAAAATATGAAGTTATTTAAAAAAAAAGAATATTAAGAAGTAAATTAATGCAAATAAAAAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAt                                        <  2:694648/148‑1 (MQ=255)
                                                                                                       gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtt                                       >  1:622923/1‑149 (MQ=255)
                                                                                                          aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGCAGAGTTCAATATAATGTTTGTCTTCAAttttc                                    >  2:75388/1‑148 (MQ=255)
                                                                                                          aaGAACCTTCAACATGAAAAATATCAATATGAATGAAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt                                    <  2:345043/149‑1 (MQ=255)
                                                                                                             aaCCTTAAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt                                 <  1:314040/149‑1 (MQ=255)
                                                                                                               ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCAAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAc                               >  2:544235/1‑149 (MQ=255)
                                                                                                                cTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTAAATATAATGTTTGTCTTCAATTTTTATTACt                              >  1:466921/1‑149 (MQ=255)
                                                                                                                 ttCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAGTTTTTCTTCCtt                             >  2:645314/1‑149 (MQ=255)
                                                                                                                           aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAata                   >  1:85381/1‑149 (MQ=255)
                                                                                                                                 tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCACTATAATGTTTGTCTTAAATTTTTCTTACTTCAGGGTAATATAGAtt             >  2:259029/1‑149 (MQ=255)
                                                                                                                                    aTTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCGGGCTAATATAGATTGCt          >  2:719160/1‑149 (MQ=255)
                                                                                                                                            gCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGATTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAt  <  1:203159/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAT  >  NC_000913_3_hsa_tpiA/1197576‑1197864

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