Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_pgi 4,182,791 A→G 54.3% D516G (GAC→GGC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi4,182,7910AG54.3% 11.4 / 46.6 46D516G (GAC→GGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (15/6);  new base G (9/16);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

GGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAA  >  NC_000913_3_hsa_pgi/4182648‑4182939
                                                                                                                                               |                                                                                                                                                     
ggTACGTGTAGAGCGTGGGGTGAAAGGGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGa                                                                                                                                                  >  2:867217/1‑149 (MQ=255)
 gTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCGGCCAGCTGTCTCAGTTTATGGaaatgat                                                                                                                                                 >  2:905518/1‑143 (MQ=255)
       gtAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCAAGCTGTCTCAGATTATGGACCAGAACAAccc                                                                                                                                           >  1:532869/1‑149 (MQ=255)
              gTGCGGGTAAAGAGCGTCTGTTTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCATTGAGAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCt                                                                                                                                    <  1:191305/149‑1 (MQ=255)
                gccgTGAAAGCGCCTCTGTCTCTTGGCGAGCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTga                                                                                                                                  <  2:757961/146‑1 (MQ=255)
                 cGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTgag                                                                                                                                 <  1:659778/149‑1 (MQ=255)
                  tgggAAAGGGCGTCCGTCTCTGGGCGACCTGGCTACCCTACTGCCACAGTATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTgaga                                                                                                                                <  2:178348/146‑1 (MQ=255)
                      aaaGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTAc                                                                                                                            >  1:516631/1‑149 (MQ=255)
                         gagCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCAGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGca                                                                                                                         >  2:648611/1‑149 (MQ=255)
                               ctgtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAAcg                                                                                                                   >  2:649978/1‑149 (MQ=255)
                                 gtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgtc                                                                                                                 <  1:685033/149‑1 (MQ=255)
                                    tctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAGGAGTTCTTCGGTTGCAGCCAGCTGTATCAGCCTATTGACCAGAACAAGCCGATATCTGAGATTACGCACAAACGTCGTa                                                                                                              >  1:943275/1‑149 (MQ=255)
                                      tGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCCCTGTCTGAGATTCCGCCCAAACGTCGTAtc                                                                                                            >  2:108152/1‑149 (MQ=255)
                                          cGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAAAACAACCCGCTGTCTGAAATTACGCACAAACTTCGTATCTCCg                                                                                                       >  1:826955/1‑149 (MQ=255)
                                            aTCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCCGCCCGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCa                                                                                                     <  1:642206/149‑1 (MQ=255)
                                             tCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAATTTATGGAACAGAACAACGAGGTGTCTGAGATTACGCACAAACGTCTTATCTCCGCAc                                                                                                    >  1:276406/1‑149 (MQ=255)
                                                 gATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCATCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCgg                                                                                                >  1:488841/1‑149 (MQ=255)
                                                  aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGc                                                                                               <  2:775453/149‑1 (MQ=255)
                                                  aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGc                                                                                               <  1:746969/149‑1 (MQ=255)
                                                    aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACTCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGccc                                                                                             <  1:905518/149‑1 (MQ=255)
                                                    aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGATTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGAACAAACGTCGTATATACCCACTCGGccc                                                                                             >  2:536715/1‑149 (MQ=255)
                                                       cTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAgg                                                                                          >  1:924181/1‑149 (MQ=255)
                                                              cacaGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACGCGGCCCAGGCGGTCTg                                                                                   >  2:767263/1‑149 (MQ=255)
                                                                   gATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACaaa                                                                                                                     >  1:918398/1‑111 (MQ=255)
                                                                   gATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACaaa                                                                                                                     <  2:918398/111‑1 (MQ=255)
                                                                       tGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGaa                                                                          <  2:266930/149‑1 (MQ=255)
                                                                         aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAt                                                                                                             <  1:854155/113‑1 (MQ=255)
                                                                         aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAt                                                                                                             >  2:854155/1‑113 (MQ=255)
                                                                          tCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAACCGAAGAGCACACGCATGCGGCCAAGGAGGGCTGACCCGGGAACGt                                                                       >  1:276078/1‑149 (MQ=255)
                                                                           cAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTg                                                                      <  1:899659/149‑1 (MQ=255)
                                                                                  aaGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCGAGCCAGCTGTCTCAGTGTAGGGGCCAGGACAACCCTCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACGCGGCCCGGGGGTGGTGGCCCGTGAACGTGCAAGcgt                                                               >  2:690079/1‑147 (MQ=255)
                                                                                            ccGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGGCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                     <  1:130465/149‑1 (MQ=255)
                                                                                                  cagTGAAAGATTTCTTCCGTTCCACCCAGCTGTCTCAGTTTATGTGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                               <  1:802401/149‑1 (MQ=255)
                                                                                                  cagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTCCGCACAAACGTCGTATCTCCGCACTCGCCCCAGGCGGTCTAACCCGTAAACGTGCAGGCTTCGAAGTTCGAGATGTa                                               >  1:51724/1‑149 (MQ=255)
                                                                                                       aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTAAAccc                                          >  1:844794/1‑149 (MQ=255)
                                                                                                         agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGCCATCCAAATTCGACACCTGCACCCGa                                        >  2:299139/1‑149 (MQ=255)
                                                                                                                  cGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACTCCACTCACTACg                               >  2:275074/1‑149 (MQ=255)
                                                                                                                        cagccagcTGTCTCAGTTTATGGGCCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                         >  2:879423/1‑149 (MQ=255)
                                                                                                                        cagccagcGGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGGCGTATCTCCGCACTCGGCCCAGGCGGTCAGATCCGGCAACGTGTAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                         <  1:195899/149‑1 (MQ=255)
                                                                                                                              gcggTCCCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc                   <  1:921236/146‑1 (MQ=255)
                                                                                                                                   ttcAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc              <  1:979042/148‑1 (MQ=255)
                                                                                                                                   cttaGTTTATGTACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCTAGACGTACACCCGACTCACTTCGGTCGCGTATGTCCAATc              <  2:651930/146‑1 (MQ=255)
                                                                                                                                   ctcAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc              <  1:946353/149‑1 (MQ=255)
                                                                                                                                    tcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCTACTCACTACCGTCGCGTATGTCCAATCg             >  2:53709/1‑149 (MQ=255)
                                                                                                                                       gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa          <  1:649978/149‑1 (MQ=255)
                                                                                                                                       atttATGGGCCAGAACAACCCGCGGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa          <  2:948242/148‑1 (MQ=255)
                                                                                                                                          tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAACTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAACGAAAccc       >  2:709652/1‑149 (MQ=255)
                                                                                                                                               aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCACGGGATCCGTGCCGTGACCGTCCGAGCTTCGAAGTACGAGTCATACATCCCACTCACTACGGACGCAAATGTCCAATCGAAACCCCCGaa  >  1:487155/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAA  >  NC_000913_3_hsa_pgi/4182648‑4182939

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: