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,7700TG97.8% 144.1 / ‑4.1 45D866E (GAT→GAGrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (18/26);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

TGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173706‑4173832
                                                                |                                                              
tGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGaa                                                             >  1:3773983/1‑68 (MQ=255)
 gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:103901/68‑1 (MQ=255)
 gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:2769837/68‑1 (MQ=255)
     ccGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCgg                                                        >  1:1117613/1‑68 (MQ=255)
     ccGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTCGAGGAATCCgg                                                        <  1:3264586/68‑1 (MQ=255)
         aaGAGATCACCGCTGAGATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATc                                                    <  1:2976632/68‑1 (MQ=255)
          agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAGTGGAGGAATCCGGTATCg                                                   <  1:2200259/68‑1 (MQ=255)
          agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:248086/68‑1 (MQ=255)
          agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:3487082/68‑1 (MQ=255)
          agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAACCGGTATCg                                                   >  1:5618898/1‑68 (MQ=255)
           gagaTCACCGCTGACATCCCTAACGTGTGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGt                                                  <  1:5495891/68‑1 (MQ=255)
             gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt                                                <  1:319848/68‑1 (MQ=255)
             gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt                                                <  1:2953839/68‑1 (MQ=255)
              aTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTa                                               >  1:3318245/1‑68 (MQ=255)
               tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc                                              <  1:2051163/68‑1 (MQ=255)
               tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc                                              >  1:1514885/1‑68 (MQ=255)
                cACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACa                                             <  1:5194898/68‑1 (MQ=255)
                 aCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACAt                                            >  1:3930812/1‑68 (MQ=255)
                  ccGCTGACATCCCGAACGTGGGTGAAGCTGCGCTGTCCAAACTGGAGGAATCCTGTCTCTTATACaca                                           >  1:829288/1‑66 (MQ=255)
                    gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg                                         >  1:2762522/1‑68 (MQ=255)
                      tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg                                       >  1:666624/1‑68 (MQ=255)
                          aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGCATCGTTTACATTGGTGCGGa                                   <  1:1031040/68‑1 (MQ=255)
                            cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 >  1:2133322/1‑68 (MQ=255)
                            cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 >  1:1413315/1‑68 (MQ=255)
                            cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 <  1:1175522/68‑1 (MQ=255)
                             ccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                <  1:830413/68‑1 (MQ=255)
                               gAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGa                              <  1:4439/68‑1 (MQ=255)
                                  cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:89875/1‑68 (MQ=255)
                                  cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:5633471/1‑68 (MQ=255)
                                  cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:2296559/1‑68 (MQ=255)
                                   gTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                          <  1:5625353/68‑1 (MQ=255)
                                       gTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc                      <  1:5272949/68‑1 (MQ=255)
                                         gAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATAGTTTACTTTGGTGCTGATGTGCCCGGTGGCGa                    >  1:4188391/1‑68 (MQ=255)
                                          aaGCTGCGCTCTCCAAACTGGAGGAAACCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc                   <  1:1889990/68‑1 (MQ=255)
                                            gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 <  1:4973068/68‑1 (MQ=255)
                                             cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                >  1:5418139/1‑68 (MQ=255)
                                             cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                <  1:3318420/68‑1 (MQ=255)
                                                    ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:662239/68‑1 (MQ=255)
                                                    ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:2501615/68‑1 (MQ=255)
                                                    ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:3111780/68‑1 (MQ=255)
                                                    ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGGGACCGGTGGCGACATTCtggttg         <  1:5118337/68‑1 (MQ=255)
                                                           cTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt  >  1:5104871/1‑68 (MQ=255)
                                                           cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggttaggt  >  1:1614745/1‑68 (MQ=255)
                                                           cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt  >  1:1351355/1‑68 (MQ=255)
                                                           cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt  <  1:102865/68‑1 (MQ=255)
                                                                |                                                              
TGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173706‑4173832

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