Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,928 A→T E562V (GAA→GTA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9280AT97.1% 99.7 / NA 34E562V (GAA→GTA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (20/13);  minor base C (0/1);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTG  >  NC_000913/4182783‑4183065
                                                                                                                                                 |                                                                                                                                         
tttATGGACCAGAACAACCCGCTGTCTGAGCTTACGCACAAACGTCGTATCTCCGCACGCGGCCCAGGCCGACTCACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCCCTCCGGTCGCGTATGTCCCATCGTAAc                                                                                                                                        >  1:62274/1‑149 (MQ=255)
     ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTCAACCCGACGCACTACGGTCGCGCATGTCCAAGCGTAACCCCTg                                                                                                                                   >  2:147212/1‑149 (MQ=255)
      gaccagaccaCCCCGCTGTCTGAGATTACGCACCAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                  <  1:155119/149‑1 (MQ=255)
      gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                  >  2:379819/1‑149 (MQ=255)
         cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTTCGGTCGCGTATGTCCAATCGTAACCCCTGAAgg                                                                                                                               <  1:287559/149‑1 (MQ=255)
                  ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACAt                                                                                                                      >  1:202218/1‑149 (MQ=255)
                      tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCATAACCCCTGAAGGTCCGAACATCGGt                                                                                                                  <  1:92995/149‑1 (MQ=255)
                        tCTGAGATTACGCACAAACGTCGTATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCt                                                                                                                <  2:427602/149‑1 (MQ=255)
                           gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCGGAt                                                                                                             >  1:63002/1‑149 (MQ=255)
                                        aaagtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                                                                >  2:99157/4‑149 (MQ=255)
                                        aaCGTTGTATCCCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                                                                <  1:245887/149‑1 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCGGTCTGATCAACTCTCTGTcc                                                                                                <  1:331972/149‑1 (MQ=255)
                                            tcgtATCTCCGCACTCGGACCAGACAGTCTGACCCGTCGACGTACAGGCTTCGAATTTCTAGGCGTACACCCGACTCACGACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                            >  2:262063/1‑149 (MQ=255)
                                            tcgtATCTCAGCACTCGGCCCAGAAAGTCTGGCCAGTGAACGTACGGGCGTCGAGGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                            >  2:262049/1‑149 (MQ=255)
                                            tcgtATCGCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCTAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTCAAGGTCCGAACATCGGTCTGACCGACTCTCAGTCCgtgt                                                                                            >  2:142395/1‑149 (MQ=255)
                                                 tctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                                                       >  1:431876/1‑149 (MQ=255)
                                                   tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcaca                                                                                     >  1:390149/1‑149 (MQ=255)
                                                          tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGATCCGAACATCGGCCTGATCAACTCTCGGTCCGGGTACGCACAGACTAAc                                                                              >  2:352112/1‑149 (MQ=255)
                                                                  gggcgTCTGACCCGGGAACGTGCAGGCTTCGAAGCTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                                                      <  2:267987/145‑1 (MQ=255)
                                                                       tCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTGCGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGCACGCACAGCCTACCGAATACGGCTTcc                                                                 >  1:254748/1‑149 (MQ=255)
                                                                          gACCCGTGAACGGGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                              <  2:447207/149‑1 (MQ=255)
                                                                             ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                                                           >  2:422590/1‑149 (MQ=255)
                                                                                tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTGGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAGGAATACGGCTTCCTTGGGACTc                                                        >  2:157076/1‑149 (MQ=255)
                                                                                tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                                                        >  1:374476/1‑149 (MQ=255)
                                                                                 gAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGAACAAATCTCTGTCCGTGTACGCACAGACTAACGAATACGGCGTCCTTGAGACTcc                                                       >  1:393532/1‑149 (MQ=255)
                                                                                           tcgtCGAAGTTCGATACGTACACCCGTCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                             <  2:122776/146‑1 (MQ=255)
                                                                                                             tACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTa                           >  2:282884/1‑149 (MQ=255)
                                                                                                              acacCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                          <  1:363211/149‑1 (MQ=255)
                                                                                                                              ggTGGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCGGTCCGTGTAGGCACAGACTAACGAATACGGCTGCCTTGCGACTCCCGATCTTAAAGTGAACGACGCTATTTTTACTGAGCAAATACCCTa          >  2:47920/1‑149 (MQ=255)
                                                                                                                               gTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc         <  1:435511/149‑1 (MQ=255)
                                                                                                                               gTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc         <  1:435507/149‑1 (MQ=255)
                                                                                                                                 cgcgTATGTCCAACCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCt       <  1:319399/149‑1 (MQ=255)
                                                                                                                                     tATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATCCACTACCTGCCt   >  2:168556/1‑149 (MQ=255)
                                                                                                                                      aTGTCCAATCGCAACCCCTGAAGGTCCGAACATCGGTCTGACCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTg  <  1:2857/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTG  >  NC_000913/4182783‑4183065

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: