Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 4,332,850 G→T 8.6% R135R (CGC→CGA)  basS ← sensory histidine kinase in two‑component regulatory system with BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,332,8500GT 8.6% 112.0 / 2.8 48R135R (CGC→CGA) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base T (2/2);  total (26/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 = 2.81e-01

ATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGC  >  NC_000913_3_hsa_tpiA/4332703‑4332996
                                                                                                                                                   |                                                                                                                                                  
aTCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGTTCAGGc                                                                                                                                                   >  1:542272/1‑149 (MQ=255)
     gCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                                                                                              >  1:482513/1‑149 (MQ=255)
         tGCCAGTAAGGGTGGTGCATCAATGTGATGCGTTTGCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAGACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                                                                                                          <  2:390845/149‑1 (MQ=255)
         tGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                          <  1:152315/149‑1 (MQ=255)
         tGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                                                                                                          >  1:589668/1‑149 (MQ=255)
            cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGCTGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTCCTGGTCAGGCGACTGAcccaat                                                                                                                                       >  2:87649/1‑145 (MQ=255)
                   ggTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCATATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTaa                                                                                                                                <  1:533535/149‑1 (MQ=255)
                           aTCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaaa                                                                                                                        >  2:588108/1‑149 (MQ=255)
                                tgtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                          >  2:38912/1‑126 (MQ=255)
                                tgtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                          <  1:38912/126‑1 (MQ=255)
                                tgtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTAGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGGCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                                                                                                   >  2:216496/1‑149 (MQ=255)
                                tgtgATGCGTTGTCGCCAGCAGTCCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                                                                                                   <  1:102905/149‑1 (MQ=255)
                                gggtaTGCGTTTACGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCATATTGAATAACACCGAGaccacc                                                                                                                   <  2:304856/145‑1 (MQ=255)
                                   gATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                                                                                                >  1:418389/1‑149 (MQ=255)
                                    aTGCGTTTACGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGGTAGCAGTTCGGGCGCCACGTCCGCGGTAAACAACCTTTCGTTATCCCGCGTGCTGGTCAGACGACTAACCACCTCGTTCAACCCCCGAACCACCGcgt                                                                                                               >  2:558042/1‑147 (MQ=255)
                                     ttcGTTTTCGCCAGCAGTTCAAGGTGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                                                                                                              <  1:60971/147‑1 (MQ=255)
                                         tcttCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTCCGTGCGCGTCGTCAGCGGTAAACAACCTTTCGTTATCCAGCTTACTGGTCAGGCGACTGACCAGCCCATTTCAGGTCGAAATCACCGCTCCGttt                                                                                                          >  2:395213/3‑149 (MQ=255)
                                               ccAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCACTTCGATTTCGAgg                                                                                                    >  1:351955/1‑149 (MQ=255)
                                                cagcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                                                                                                   >  1:675328/1‑149 (MQ=255)
                                                cagcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCACGGAATTTGTTCAGTCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTCAGATTTCGAggg                                                                                                   >  1:157917/1‑149 (MQ=255)
                                                  gcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCTGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                                                                                                 <  1:449308/149‑1 (MQ=255)
                                                  gcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGTTg                                                                                                 >  1:626966/1‑149 (MQ=255)
                                                      ttCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgggc                                                                                             >  1:292474/1‑149 (MQ=255)
                                                      ttCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                                                                                             >  1:257360/1‑149 (MQ=255)
                                                                   cGCAGCCCCGGCAGGGGCGTTCGCAGTTCGGGCGCGACGTCAGGGGTAAACAACCTTTCGTTATCGAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt                                                                                <  2:430035/149‑1 (MQ=255)
                                                                     taCCCCCGACAGTGACCTTTGCAGTTCGTGCGCGACGTCAGCGGGAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATgg                                                                              <  2:52858/148‑1 (MQ=255)
                                                                            gCCCGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTaa                                                                       <  1:597730/149‑1 (MQ=255)
                                                                                  ggCGTTCGCAGTGCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTc                                                                 <  1:199182/149‑1 (MQ=255)
                                                                                     gTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTGAAGTTGTcggc                                                              >  1:279318/1‑149 (MQ=255)
                                                                                            gTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc                                                       <  1:543046/149‑1 (MQ=255)
                                                                                                gcgcgcGACGGCAGCGGTAAACACGCTTTCGTTATCCAGCGTGCGGGGCAGGCGGCTGAGCAGGTCGGTTAACGGCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGGGAAGGGCAATGGGCGTGAAGTTGTCGGCGGTGCGCGCtg                                                   <  2:93111/147‑2 (MQ=255)
                                                                                                   cgcgAGCTCGGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAGTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCa                                                <  2:317500/149‑1 (MQ=255)
                                                                                                   cgagACGTCAGCGGTAAACAACCTTGTGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAAAGCCGAAAAAACAGACTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGGCGGCGGTGCGCGCTTCCa                                                <  2:316734/146‑1 (MQ=255)
                                                                                                     cgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTTACCTGAGCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGc                                              >  2:632247/1‑149 (MQ=255)
                                                                                                        cGTCAGCGGTAAACACCCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGctct                                           <  1:318248/149‑1 (MQ=255)
                                                                                                             gCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGGTTTCGAGGGTGGCGCCGTGAATCGCAACGGGCGTTACGCTGTCCGCGGTGCGCGCTTGCAGCGCTTTTTg                                      >  2:88634/1‑149 (MQ=255)
                                                                                                             gCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTg                                      <  2:508696/149‑1 (MQ=255)
                                                                                                              cGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTGAGTTGTGGGCGGTGGGCGCTTCAAGCTGTTTTTga                                     >  2:338939/1‑148 (MQ=255)
                                                                                                                ggaaaCAACCTTTCGTTATCCAGTGTGCTGGTTAGACGACTGTCCAGATCGGTTTACGCCGAAACCACCGCCTCGATTTCGAGGGTGGCGCTGTGAAGGGCAGTGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTGTGCAg                                   <  2:474467/147‑1 (MQ=255)
                                                                                                                   aacaacCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCGTCGATTTCGAGGGTGCCGCCGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCACCTc                                >  2:24416/1‑149 (MQ=255)
                                                                                                                    acaacCTTTCGTTATCCAGCGTTCTGCTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGCGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCCCTTTTTGCAGCTcc                               <  2:358191/149‑1 (MQ=255)
                                                                                                                    acaacCGTTCGGTATCCAGCGGGCGGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTcc                               <  2:287744/149‑1 (MQ=255)
                                                                                                                          tttCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgc                         >  1:668842/1‑149 (MQ=255)
                                                                                                                          tttCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAATTCGCGGGTCACGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCGCTTATTGTCGCTccgcccgc                         >  2:2699/1‑149 (MQ=255)
                                                                                                                           ttCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcg                        <  1:106168/149‑1 (MQ=255)
                                                                                                                                    cAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAGGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGa               >  1:180593/1‑149 (MQ=255)
                                                                                                                                             ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCGGTGAATGGCACTGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGCCGGGTGATGGGGCGtg      >  2:352220/1‑148 (MQ=255)
                                                                                                                                                 aGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGc  >  1:686854/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGC  >  NC_000913_3_hsa_tpiA/4332703‑4332996

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