Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,791 A→G 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_0009134,182,7910AG97.4% 124.2 / ‑5.2 38D516G (GAC→GGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (18/19);  total (19/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

TGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAA  >  NC_000913/4182654‑4182925
                                                                                                                                         |                                                                                                                                      
tgtAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAAcc                                                                                                                             >  5:39707/1‑149 (MQ=255)
  tAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                           >  8:133907/1‑149 (MQ=255)
             gTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCt                                                                                                                     >  8:75542/1‑144 (MQ=255)
             gTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCt                                                                                                                     <  7:75542/144‑1 (MQ=255)
              tGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTg                                                                                                                    >  2:123601/1‑144 (MQ=255)
              tGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCACCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTg                                                                                                                    <  1:123601/144‑1 (MQ=255)
                   gagCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGTGATTACGca                                                                                                          <  6:39707/149‑1 (MQ=255)
                   gagCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGca                                                                                                          <  2:56212/149‑1 (MQ=255)
                                       tCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTgag                                                                                                                  >  6:76406/1‑121 (MQ=255)
                                       tCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTgag                                                                                                                  <  5:76406/121‑1 (MQ=255)
                                           gATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCgg                                                                                  >  4:114478/1‑149 (MQ=255)
                                            aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGCTTCCAGCCAGCTGTCTCAGATTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGc                                                                                 <  7:133907/149‑1 (MQ=255)
                                             tACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGcc                                                                                >  2:118866/1‑149 (MQ=255)
                                                   gATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGGCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCg                                                                          <  1:131597/149‑1 (MQ=255)
                                                    aTGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGCCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                                                         >  2:81815/1‑149 (MQ=255)
                                                         acaGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGa                                                                    >  4:102442/1‑149 (MQ=255)
                                                                  gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGGTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCTTGAAc                                                           >  5:98601/1‑149 (MQ=255)
                                                                     caccGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTg                                                        >  7:123954/4‑149 (MQ=255)
                                                                      aaCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCt                                                                      >  1:13461/1‑134 (MQ=255)
                                                                      aaCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCt                                                                      <  2:13461/134‑1 (MQ=255)
                                                                          ccAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTAAGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGCTCTGACCCGTGAACGTGCAGGc                                                   >  5:43633/1‑149 (MQ=255)
                                                                            aaGCCGATATCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGCTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCtt                                                 >  2:144127/1‑149 (MQ=255)
                                                                              gCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg                                               <  6:43633/149‑1 (MQ=255)
                                                                                cGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGCCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCTGGCTTCGaa                                             <  1:81815/149‑1 (MQ=255)
                                                                                  aTTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                                                         >  8:105040/1‑119 (MQ=255)
                                                                                  aTTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                                                         <  7:105040/119‑1 (MQ=255)
                                                                                  aTTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGt                                           >  3:108517/1‑149 (MQ=255)
                                                                                       cGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                       <  3:114478/148‑1 (MQ=255)
                                                                                             agTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                 >  1:87941/1‑148 (MQ=255)
                                                                                                   agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACACCCCTCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAGCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                          <  3:30206/149‑1 (MQ=255)
                                                                                                    gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                         <  8:123954/149‑1 (MQ=255)
                                                                                                            cGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                   <  1:118866/147‑1 (MQ=255)
                                                                                                             ggTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                 <  5:106208/148‑1 (MQ=255)
                                                                                                               ttCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                 >  2:5382/1‑146 (MQ=255)
                                                                                                               ttCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                 <  1:5382/146‑1 (MQ=255)
                                                                                                               gtCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACACCCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc              <  2:104892/148‑1 (MQ=255)
                                                                                                                        gcTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc     <  1:144127/149‑1 (MQ=255)
                                                                                                                           gTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa  >  6:110227/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAA  >  NC_000913/4182654‑4182925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑