Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_pgi 3,440,084:1 +G 65.2% coding (946/990 nt) rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi3,440,0841.G65.2% 71.5 / 57.6 46coding (946/990 nt)rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base . (6/10);  new base G (15/15);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAA  >  NC_000913_3_hsa_pgi/3439936‑3440221
                                                                                                                                                      |                                                                                                                                         
ggAACACAGCCTGGCGAGGGATGCTGTTGCGGTTCAGTGGTCGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCAt                                                                                                                                            <  2:331411/149‑1 (MQ=255)
 gAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGCTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATg                                                                                                                                          >  2:226420/1‑149 (MQ=255)
      tggCCGGGGGATGAGTGCTGTTGCTGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGATTACTCGTCAGGGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGgccca                                                                                                                                     <  2:290736/147‑3 (MQ=255)
      tAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGTCGGTGGCCAGTTTTCCAGGCGCATggccca                                                                                                                                     <  1:308947/149‑6 (MQ=255)
           tGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCCTGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATACTTGCCGGTGCCCAGTTTTCCAGGCGCATGGCCCCGAGAc                                                                                                                                 >  2:98807/1‑149 (MQ=255)
           cggggacggCTGCTGTCGCGGTTGAGTTGACGACCACTCTTATGATGGCGCATGTCCTTATCCTTCTCAGTGAAACGTTAACCTGTCATCCGGTCACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACa                                                                                                                                <  2:959/142‑1 (MQ=255)
            ggCGATGGCTGCTGTTGCGGTTCAGTTGACGTCCACTCTTACGATGGCGCATGACCTTCTCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACa                                                                                                                                <  2:434868/149‑1 (MQ=255)
            gcggATGGCGGCTGTTGCGCGTCAGCTGACGCGCACTCTTGCGATGGCGCATGTGCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTTCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACa                                                                                                                                <  2:260183/146‑1 (MQ=255)
                   gctgctGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACg                                                                                                                         >  1:247855/1‑149 (MQ=255)
                        tgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCCTGGCCCAGAGCCAGTCCACGGGAAg                                                                                                                    >  2:60973/1‑149 (MQ=255)
                          ttGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGAACTTATCCTTCTCAGTAAACCCTTAACCTGTGATCCGCTTACTCGTCAGCGACGCTTGTCGGTGGCCAGTTTTCCACGCGCCTGGCCTGGAGCCAGTCCACGGTCAGcc                                                                                                                  >  2:395646/1‑149 (MQ=255)
                           tGCGGTTCAGTTGACGACCACTCTTACGATGGCGCCTGTCCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTGAGCGATGCTTGCCGGTGGCAAGTTTTCCATGCGCATGCCCAGAGACAGTCCACGGGAAGCCAg                                                                                                                <  1:186384/149‑1 (MQ=255)
                                   aGTTGACGACCACTCTTACGATGGCGCATGACTTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCt                                                                                                         <  2:649252/149‑1 (MQ=255)
                                    gTTGACGACCACTCGTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACTCGAACCTGTCATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCtt                                                                                                        <  2:455214/149‑1 (MQ=255)
                                          gacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCGGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACGGGCGCGGCGAAGCGGGCAGGTCTTTaagccc                                                                                                 >  1:2016/1‑145 (MQ=255)
                                          gacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCAACGGGAAGCCAGCACGTCTTTAAtct                                                                                                  >  2:457298/1‑149 (MQ=255)
                                                tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAag                                                                                           >  2:333860/1‑149 (MQ=255)
                                                tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGACCCGGTTACTCGTCAGCGATGCTTGCCGATGGCCAGTTTTCCAGGACCATCCCCAGCGAGCGACCACCGGAAACCACCAGCACTTTAATCTCacgcaa                                                                                           >  2:524251/1‑144 (MQ=255)
                                                  ttACGATAGCGCATGACCCTATCGTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGGGGCCAGTTTTCCAGGCGAATGGCCCAGAGACGATCCACGGGAAGCGAGCACGTCTTTAATCTCAGTAaga                                                                                          >  1:615777/1‑149 (MQ=255)
                                                                  ccTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGt                                                                          >  2:377651/1‑149 (MQ=255)
                                                                  ccTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGtt                                                                         >  1:418147/1‑149 (MQ=255)
                                                                      atgcTTCTCAGTAAAACCTTAACCTGTGATACGGTCACTCGGTATCGATGCTTGCCGGTGGCCAGTTGTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGc                                                                     <  2:534463/146‑1 (MQ=255)
                                                                            ctcAGTAACACCTCAACCGGGGGTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGCCTTTAATCTAAGTAAGAGATTTTTTACAATGGTTAGGAGgttt                                                                >  2:412567/1‑149 (MQ=255)
                                                                               aGTAAAACCTTAACCTTGTCTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGAAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                             <  2:541147/149‑1 (MQ=255)
                                                                               aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                             <  2:413457/149‑1 (MQ=255)
                                                                                gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGa                                                           <  2:10493/149‑1 (MQ=255)
                                                                                 tcgaaCCTTAACCCGTGATCCGGTGACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCAGGGCCCAGAGACAGTCCCCCGGACGCTAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGa                                                           <  1:416776/146‑1 (MQ=255)
                                                                                   aaaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTACGGAggt                                                        >  2:586229/1‑147 (MQ=255)
                                                                                      ccTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCCCAGTAAGACATTTTTTACCAAGGGTAGGCGTTTTAAGGAGCTCa                                                      >  2:247816/1‑149 (MQ=255)
                                                                                            ccggtgATCCGCTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCg                                                <  1:402627/146‑1 (MQ=255)
                                                                                             cTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACCGTCCACGGGAAGCCCGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTCGGCGTTTTAAGGAGCTCAACCTCgc                                               >  2:319185/1‑148 (MQ=255)
                                                                                              tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACGAAGGTTAGGGGTTTTGAGGCGCTCAACCTCGGTa                                             >  1:501627/1‑149 (MQ=255)
                                                                                                  aTCCGGTTCCTCCTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCt                                         <  1:247008/149‑1 (MQ=255)
                                                                                                     cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGt                                       >  1:608928/1‑149 (MQ=255)
                                                                                                      cgTTACTCGGGTGCGATGCTGGCCGTTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGGACGTCTTTAATCTCAGTAAGAGTTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                      <  2:19940/148‑1 (MQ=255)
                                                                                                       gTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTAc                                     >  1:524720/1‑149 (MQ=255)
                                                                                                         tACTCGTGAGCGATGGTTGCCGCGGGCCAGGTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCa                                   <  2:63841/149‑1 (MQ=255)
                                                                                                           cTCGTCGGCGATGCTTGCCGGTGGCCAGTTCTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACCGCTTTCATCTCAGTAAGTGATTTTTTACCGAGGTTAGGCGTTTTAAGTAGCTCAACCTCGGTACGCTGTACCAGAt                                <  2:633124/149‑1 (MQ=255)
                                                                                                               tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCCGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTATTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                             >  1:352018/1‑149 (MQ=255)
                                                                                                               tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCCCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGGGCTCAACCTCGGTACGCTGTACACGATCAc                             >  1:357166/1‑149 (MQ=255)
                                                                                                                catagATGCTTGCCGGTCGCACGTTTTCCAGGCGCATGGCCCAGGGACAGTCCACGGTAAGCCAGCACGTCTTTACTCTCATTAAGGTATTTTTTACAAAGGTTAGGCGTTTTAAGGGGCTCAACCTCGTTAAGCTGTAGCAGATCAcc                            <  2:212963/145‑1 (MQ=255)
                                                                                                                cAGCGATGCTTGCCGGTGGCGAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCAGGGGAAGCCAGCACGTCTTTAATCTCAGTAATAGATTTTTTACCAAGGTTAGGTGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAcc                            <  2:21229/149‑1 (MQ=255)
                                                                                                                    gATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGGCCCAGAGAAAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGCAAGCGATTTGTTACCAAGGTTAGGCGTTTTAAGGAGCTCACCCTCGGTACGCTGTACCAGATCCCCGata                        <  2:693/149‑1 (MQ=255)
                                                                                                                                 ggCCAGTTTGCCAGGCGCATCCCCAGAGTCGGTCCACGGGAAGCCCGAACGTCTTTAATCTCAGTAAGCGATTTTTTACCACGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTc          <  1:255085/149‑1 (MQ=255)
                                                                                                                                  gCCAGTTTTCCAGGCGGATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAATGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTc          <  2:512349/149‑1 (MQ=255)
                                                                                                                                   ccAGTTTTCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCt         >  1:603693/1‑149 (MQ=255)
                                                                                                                                         ttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGCCTCACCTTTCGAGGAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa  <  1:519937/149‑1 (MQ=255)
                                                                                                                                          tCCAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa  <  1:480447/149‑1 (MQ=255)
                                                                                                                                          gccAGGCGCATGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa  <  1:60973/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                         
GGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAA  >  NC_000913_3_hsa_pgi/3439936‑3440221

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