Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,173,770 T→G D866E (GAT→GAG rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,173,7700TG98.0% 195.0 / ‑5.0 51D866E (GAT→GAGrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (27/23);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCT  >  CP009273/4173635‑4173905
                                                                                                                                       |                                                                                                                                       
tgttgtTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGGGGGTGAAGCTGCGCTCTCCAAACTGGAGGaa                                                                                                                                      >  1:4283/1‑139 (MQ=255)
          gAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTggagga                                                                                                                                       >  1:2692732/1‑128 (MQ=255)
              aCCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt                                                                                                                         >  1:1959292/1‑138 (MQ=255)
                 gTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACAt                                                                                                                     <  1:1908840/139‑1 (MQ=255)
                  tttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACAtt                                                                                                                    >  1:1937892/1‑139 (MQ=255)
                          ccaTCCACATTCAGGAACTGGCGTGGGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCCCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGGGTCGGTTACATTGGTGCGGa                                                                                                            >  1:18528/1‑139 (MQ=255)
                             tCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                                                                                         >  1:924496/1‑139 (MQ=255)
                              ccACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGt                                                                                                                                <  1:1862408/115‑1 (MQ=255)
                                    tCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTa                                                                                                                               >  1:1369093/1‑110 (MQ=255)
                                    tCAGGAACTGGCGGGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg                                                                                                                <  1:1004957/125‑1 (MQ=255)
                                        gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACa                                                                                                                      >  1:438141/1‑115 (MQ=255)
                                        gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg                                                                                              >  1:1114862/1‑139 (MQ=255)
                                             ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGACCGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                                                                                         <  1:2769623/139‑1 (MQ=255)
                                              gCGTGTGTGTCCCGTGACACCCAGCTGGGTCCGGGAGCGATCACCGCTGACATCCCGAACGTGGGTGACGCTGCGCTCTCCAACCTGGAGGAATCCGGCATCGTTTACATTGGTGCGGGAGTGAc                                                                                                      >  1:2428458/1‑125 (MQ=255)
                                              gCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                                                                                          >  1:488697/1‑137 (MQ=255)
                                                gtgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                                                      <  1:535129/139‑1 (MQ=255)
                                                  gtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggt                                                                                    <  1:1912914/139‑1 (MQ=255)
                                                   tgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                                                   >  1:123033/1‑139 (MQ=255)
                                                    gggtCCCGTGACACCAAGATGGGTCCGGAAGAGATCACCGCTGACCTCCCGAACGTGGGTGAAGCTGCGCGCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                                                   <  1:1443228/136‑1 (MQ=255)
                                                        cccGTGCCACCAAGCTGGATCCGGACGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGTATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGTTGGCGACATTCTGGGTgggaa                                                                              >  1:1980307/1‑139 (MQ=255)
                                                        cccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaa                                                                              >  1:2458886/1‑139 (MQ=255)
                                                               cacCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGt                                                                                                                                <  1:75448/82‑1 (MQ=255)
                                                                acCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACg                                                                       >  1:65628/1‑138 (MQ=255)
                                                                 ccAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGcc                                                                     >  1:1284206/1‑139 (MQ=255)
                                                                      cTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAg                                                                >  1:2815990/1‑139 (MQ=255)
                                                                         ggTCCGGAAGAGATCATCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTg                                                             >  1:1225735/1‑139 (MQ=255)
                                                                         ggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACa                                                                                                                      >  1:1500063/1‑82 (MQ=255)
                                                                          gTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTg                                                             <  1:1547607/138‑1 (MQ=255)
                                                                          gTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGa                                                            >  1:1886097/1‑139 (MQ=255)
                                                                               gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTc                                                       >  1:1253649/1‑139 (MQ=255)
                                                                                       cACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGCCATTCTGGTTGGTAAGGTAACGCCGaa                                                                  <  1:1240667/120‑1 (MQ=255)
                                                                                         ccGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                                             <  1:1534834/139‑1 (MQ=255)
                                                                                         ccGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                                             <  1:1490255/139‑1 (MQ=255)
                                                                                          cGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCa                                            >  1:110250/1‑139 (MQ=255)
                                                                                             tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCaga                                          <  1:2707221/138‑1 (MQ=255)
                                                                                              gACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGGGGCGACATTCTGGGTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCggaag                                        >  1:2832406/1‑139 (MQ=255)
                                                                                                 aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACt                                                        >  1:2663775/1‑120 (MQ=255)
                                                                                                       aaCGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct                                 >  1:2239274/1‑137 (MQ=255)
                                                                                                               tGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggt                                                                                    <  1:1163044/78‑1 (MQ=255)
                                                                                                                  aGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGCAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGa                      >  1:1664845/1‑137 (MQ=255)
                                                                                                                     tGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCt                   <  1:43229/137‑1 (MQ=255)
                                                                                                                      gcgcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCg                >  1:324716/1‑139 (MQ=255)
                                                                                                                      gcgcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCg                <  1:2648864/139‑1 (MQ=255)
                                                                                                                      gcgcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCg                <  1:1847218/139‑1 (MQ=255)
                                                                                                                       cgcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCgg               >  1:771263/1‑139 (MQ=255)
                                                                                                                         ctctCCAAACTGGAGGAATCCGGTATCGTTTACATGGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggaa                                                                          >  1:2586281/1‑76 (MQ=255)
                                                                                                                            tcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTgaga          <  1:240425/139‑1 (MQ=255)
                                                                                                                            tcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTgaga          <  1:1749055/139‑1 (MQ=255)
                                                                                                                             ccAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGaa         <  1:1393103/139‑1 (MQ=255)
                                                                                                                             ccAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGaa         <  1:782674/139‑1 (MQ=255)
                                                                                                                                 aCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                                                      <  1:974969/58‑1 (MQ=255)
                                                                                                                                   tGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCctc   <  1:461852/139‑1 (MQ=255)
                                                                                                                                    ggaggaATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCctct  <  1:556169/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCT  >  CP009273/4173635‑4173905

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: