Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 269,605 A→T 9.4% V193E (GTG→GAG)  perR ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA269,6050AT 9.4% 101.1 / 3.4 43V193E (GTG→GAG) perRCP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base T (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TCAGCCAGTAGTTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGATGGGCGGGGTTCACCGT  >  NC_000913_3_hsa_tpiA/269458‑269743
                                                                                                                                                   |                                                                                                                                          
tCAGCCAGTAGTTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGcac                                                                                                                                           <  1:679252/149‑1 (MQ=255)
    ccAGTAGTTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCGCAtgg                                                                                                                                       >  1:649649/1‑147 (MQ=255)
     cAGTAGTTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGc                                                                                                                                      >  1:623526/1‑149 (MQ=255)
          gtTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCTCATCGCACTGa                                                                                                                                 >  1:674056/1‑149 (MQ=255)
          gtTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCTCATCGCACTGa                                                                                                                                 >  1:674041/1‑149 (MQ=255)
           ttttGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGaa                                                                                                                                >  1:223262/1‑149 (MQ=255)
                aTTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAAt                                                                                                                           <  1:493131/149‑1 (MQ=255)
                             tCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATaga                                                                                                              >  1:239816/1‑149 (MQ=255)
                               ccAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAgaga                                                                                                            >  1:470851/1‑149 (MQ=255)
                                       tCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGACCCTCTgg                                                                                                    >  2:519797/1‑149 (MQ=255)
                                                gcagcaATCGACATAAAGCTACGGTCAAATCGAAGGCCATAATTATTTGGTGGCGTCATCTTATTCTCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAg                                                                                           <  2:525946/149‑1 (MQ=255)
                                                 cagcaATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAGTTATTGGGTGGCGTCATCTTATTCGCGGCAAACCATCCTTTCCCCTGGTACACCTGCACATAGCGCTGAATCGATGTCAACGCGTCGAGATACGCTGATTTCTTTatc                                                                                          >  1:291964/1‑147 (MQ=255)
                                                 cagcaATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCTCATCGCACTGAATCAATGTCAACGCAGAGAGATCCTCTGGTTTCTTTAGc                                                                                          <  2:674056/149‑1 (MQ=255)
                                                 cagcaATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCTCATCGCACTGAATCAATGTCAACGCACAGAGATCCTCTGGTTTCTTTAGc                                                                                          <  2:674041/149‑1 (MQ=255)
                                                 cagcaAACGCCATAAAGCTACGGTCAAATCGAACGCCATAATTATTGGGTGGCGTCATCTTATTCGTCTCAAACTATCCTTTCCTCTGTGACAACTGTACGTCGCACTGACTCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGc                                                                                          <  2:300002/149‑1 (MQ=255)
                                                         gCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGc                                                                                  <  2:681319/149‑1 (MQ=255)
                                                         gCCATAAAGCCACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGc                                                                                  <  1:549246/149‑1 (MQ=255)
                                                           cATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGACATCCTCTGGTTTCTTTAGCCGCTCAGCCa                                                                                <  2:721371/149‑1 (MQ=255)
                                                           cATAAAGCTACGGTCAAATCGAAGGCCATAATAATTGGGTGGAGTCATCTTATTCGCCTCAAACCATCCTTTCCAGTGGTACAACTGCACATCGCACAGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTcgacca                                                                                >  2:690532/1‑144 (MQ=255)
                                                           cATAAAACTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAAAGCATAAAGAAACAATGGTTTCTTTAGCAGCTCAGccc                                                                                >  2:661439/1‑148 (MQ=255)
                                                                atcTACGGTCAAATTGAAGTCCATATTTATTGGGTGGCGTCATCTTATTCGCTTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGCCTATGCGTAGAGATCCTCTTGTTTCTTTAGCCGCTCAGCCAGTTgg                                                                           <  1:691112/147‑1 (MQ=255)
                                                                   tACTGTCAAATCGAACGCCAGAATTATTGGGTGGCGTCATCTTATTCGCCTCAAAGCATCCTTTCCACTGGTACAACTGCATATCGCACTGAATCAATGTCAACGTATAGAGATCCTCTGGTTTCTTTAGCTGGTCAGCCAGTTGGGga                                                                        <  2:597027/149‑1 (MQ=255)
                                                                   tACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGga                                                                        <  1:281582/149‑1 (MQ=255)
                                                                      ggTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGATAGCgagat                                                                     >  2:694582/1‑148 (MQ=255)
                                                                        tCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTGAAACCATCCTTTCCACTGGTACAACTGCACATCGCATTGAATCAATGTCAACGCATAGAGATGCTCGTGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAaca                                                                   <  2:6552/149‑1 (MQ=255)
                                                                           aaTCGAAGGCCATAATTATTTGGTTGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAg                                                                <  1:277877/149‑1 (MQ=255)
                                                                              cGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGGACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTgg                                                             <  1:436589/149‑1 (MQ=255)
                                                                                     cATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGt                                                      >  2:320994/1‑149 (MQ=255)
                                                                                         attattGGGTGGCGTCATCTTATTCGCCTCAAACCAGCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGttctt                                                  <  2:94238/149‑1 (MQ=255)
                                                                                                      gTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCgg                                     >  2:351232/1‑149 (MQ=255)
                                                                                                         aTCTTATTCGCCTTAAACCATCCCTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGAt                                  <  2:419287/149‑1 (MQ=255)
                                                                                                          tCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCt                                                                                               <  2:680131/87‑1 (MQ=255)
                                                                                                          tCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCt                                                                                               >  1:680131/1‑87 (MQ=255)
                                                                                                            ttATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCGGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGAtctt                               <  2:709471/149‑1 (MQ=255)
                                                                                                                    cTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAgg                       <  1:562891/149‑1 (MQ=255)
                                                                                                                        aaCCGTACTTGCCACTGGTACAACTGCACATCGCACTTAATCAATGTCAAGGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGTCGTAATTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGa                   <  2:453632/149‑1 (MQ=255)
                                                                                                                        aaCCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGa                   <  2:623526/149‑1 (MQ=255)
                                                                                                                        aaCCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGCCAACGCAGAGAGATCCTCTGGATTCTTAAGCCGCTCGGGCAGTGGGGGAAAACAAACTGGCCTACGTGCCACACCATCAAGCGGGACCTTACCATAgggcac                   >  2:157669/1‑146 (MQ=255)
                                                                                                                         aCCATCCTTGCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGCTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGAt                  <  2:628257/149‑1 (MQ=255)
                                                                                                                               cTTTCCACTGGTACAACTGCCCATCGCAGTGAATTAATTTCAACGCATAAAGATCCGCTGGTTTCTTTCGCTGCTCAACCAGTTAGAGGGAACACATGGGCGCAAGTCCTTCAACCGCAAGCTGGATCTTAACCTAGGGTAATGGCCgg            >  2:355063/1‑149 (MQ=255)
                                                                                                                                    cACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGATGGGCGGGGTTc       >  1:99666/1‑149 (MQ=255)
                                                                                                                                     aCTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGATGGGCGGGGTTCa      >  1:221225/1‑149 (MQ=255)
                                                                                                                                       tGGTACAACTGCACATCGCACTGAATCAATGTCAGCCCATAGAGATCCCCTGGCTTCTTTAGCCGCTCAGCCAGTTGGGGAGCACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATACGGTGATGGCCGCGGTTCAcc    >  2:389875/1‑149 (MQ=255)
                                                                                                                                         gTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGATGGGCGGGGTTTACCGt  >  2:606245/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
TCAGCCAGTAGTTTTGATTCAAGCACAACTCCCAAGCCGTCAACTGCGGCAGCAATCGCCATAAAGCTACGGTCAAATCGAAGGCCATAATTATTGGGTGGCGTCATCTTATTCGCCTCAAACCATCCTTTCCACTGGTACAACTGCACATCGCACTGAATCAATGTCAACGCATAGAGATCCTCTGGTTTCTTTAGCCGCTCAGCCAGTTGGGGAGAACACAGTGGCGTAAGTTCTTCAACAGCAAGCGGGATCTTCTCATAGGGTGATGGGCGGGGTTCACCGT  >  NC_000913_3_hsa_tpiA/269458‑269743

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