Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,665,957 Δ1 bp coding (523/1281 nt) csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,665,9560C.98.0% 221.1 / ‑3.8 49coding (522/1281 nt)csiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (29/19);  total (30/19)
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

GATACCAATCTACAAGCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCTCCAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGA  >  NC_000913_3_hsa_tpiA/2665810‑2666099
                                                                                                                                                  |                                                                                                                                               
gaTACCAATCTACAAGCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCTGCTcaa                                                                                                                                               >  1:229884/1‑147 (MQ=255)
  tACCAATCTACAAGCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCTcaacg                                                                                                                                             >  1:682694/1‑145 (MQ=255)
        tCTACAAGCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTCCGCCTCAGTTTAATCCGC‑CCACCGTCGCTg                                                                                                                                      >  2:543100/1‑149 (MQ=255)
            caggCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTACCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCCGGGCg                                                                                                                                  <  2:648980/146‑1 (MQ=255)
               gCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGaa                                                                                                                               >  1:619517/1‑149 (MQ=255)
                   tAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGAAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCt                                                                                                                           <  2:574209/149‑1 (MQ=255)
                     gTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAAACGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCTGCT‑CAACGTCGCTGGGCGGATTCCTGc                                                                                                                         >  2:814616/1‑149 (MQ=255)
                          tctcTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGa                                                                                                                    >  1:46921/1‑149 (MQ=255)
                             ctGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAAtt                                                                                                                 >  1:685305/1‑149 (MQ=255)
                              tGCTCCCGCCGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAAttt                                                                                                                >  2:102661/1‑149 (MQ=255)
                              tGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAAttt                                                                                                                >  1:593905/1‑149 (MQ=255)
                               gCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTc                                                                                                               <  2:619517/149‑1 (MQ=255)
                                cTCCCGACGGCTGCAANCNCNCNTNGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCa                                                                                                              <  1:567412/149‑1 (MQ=255)
                                cTCCCGACGGCTGCAAAAACGCTTTGCAACGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGGATCCTGCCTTGAATTTCa                                                                                                              >  1:837496/1‑149 (MQ=255)
                                cTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCCGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCa                                                                                                              <  1:754199/149‑1 (MQ=255)
                                cTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCa                                                                                                              <  2:685782/149‑1 (MQ=255)
                                cTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCa                                                                                                              <  2:685305/149‑1 (MQ=255)
                                cTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCa                                                                                                              <  2:300159/149‑1 (MQ=255)
                                                aaaCGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAAtt                                                                                              <  1:152685/149‑1 (MQ=255)
                                                 aaCGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTg                                                                                             >  1:316343/1‑149 (MQ=255)
                                                      tttGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCGAGCAATGAATTTCAGGTAGCGCAGGAAATTGGACGc                                                                                        >  1:788746/1‑149 (MQ=255)
                                                                     atccctttCCTGTGTCTGTATCTGCAATATTGTTTGCCGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCGGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCgg                                                                         <  2:82017/144‑1 (MQ=255)
                                                                     atTCGTTTCCTGTTTCTGTATCTGCAATATTGTTTCCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCgg                                                                         <  2:801881/149‑1 (MQ=255)
                                                                     atTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCGGGTAGCGCAGGAAATTGGACGCCACTGGCAACGACgg                                                                         >  1:404657/1‑149 (MQ=255)
                                                                              cTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAg                                                                >  2:443337/1‑149 (MQ=255)
                                                                               tgtgTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGc                                                               >  1:761464/1‑149 (MQ=255)
                                                                                    cTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTa                                                          >  1:309844/1‑149 (MQ=255)
                                                                                     tGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGa                                                                                           >  1:502935/1‑115 (MQ=255)
                                                                                     tGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGa                                                                                           <  2:502935/115‑1 (MQ=255)
                                                                                        aTCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTaccac                                                      >  2:22029/1‑149 (MQ=255)
                                                                                         tCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTaccacc                                                     <  1:606738/149‑1 (MQ=255)
                                                                                           tGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAAACCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCtg                                                   >  1:469374/1‑149 (MQ=255)
                                                                                            gCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCtga                                                  >  2:228506/1‑149 (MQ=255)
                                                                                             cAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCtgat                                                 <  1:701041/149‑1 (MQ=255)
                                                                                              aaTATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCtgatg                                                >  2:94823/1‑149 (MQ=255)
                                                                                               aaatTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCtgatga                                               <  2:714950/147‑1 (MQ=255)
                                                                                                tatTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAg                                              >  1:232432/1‑149 (MQ=255)
                                                                                                     tttGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGGGCCACt                                         >  1:224748/1‑149 (MQ=255)
                                                                                                     tttGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACt                                         >  1:496047/1‑149 (MQ=255)
                                                                                                     tttGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACt                                         >  1:241006/1‑149 (MQ=255)
                                                                                                     tttGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCCCCAGCCTGTACCGCCTGATGAGCCACt                                         >  1:18543/1‑149 (MQ=255)
                                                                                                         ctgcAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGttt                                     >  2:157827/1‑149 (MQ=255)
                                                                                                           gcagcaCCCCGCTGGAATTACGCCCCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTAt                                   <  1:516643/149‑1 (MQ=255)
                                                                                                           gcagcaCCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTAt                                   <  1:622711/149‑1 (MQ=255)
                                                                                                           gcagcaCCACGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTAt                                   <  2:593905/149‑1 (MQ=255)
                                                                                                                 ccacGCTGGAATTACGCCTCAGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                             >  1:710347/1‑149 (MQ=255)
                                                                                                                                     aGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCggg         >  2:420054/1‑149 (MQ=255)
                                                                                                                                     aGTTTAATCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCggg         >  1:357187/1‑149 (MQ=255)
                                                                                                                                        ttACTCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTCTTCCATGTTGCGGGTTc      <  1:500932/149‑1 (MQ=255)
                                                                                                                                           aTCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCg   >  1:52831/1‑149 (MQ=255)
                                                                                                                                            tCCGCT‑CAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGa  >  1:698376/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GATACCAATCTACAAGCGTTAGTGAATCTCTGCTCCCGACGGCTGCAAAAACGCTTTGAATCGCGCGATATTCATTTCCTGTGTCTGTATCTGCAATATTGTTTGCTGCAGCACCACGCTGGAATTACGCCTCAGTTTAATCCGCTCCAACGTCGCTGGGCGGAATCCTGCCTTGAATTTCAGGTAGCGCAGGAAATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGA  >  NC_000913_3_hsa_tpiA/2665810‑2666099

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