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,7700TG92.9% 129.9 / 1.1 42D866E (GAT→GAGrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (14/25);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173709‑4173832
                                                             |                                                              
gTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATcc                                                          <  1:3059392/68‑1 (MQ=255)
     gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                     <  1:1473965/68‑1 (MQ=255)
     gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt                                                     <  1:964228/68‑1 (MQ=255)
     gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt                                                     >  1:2640675/1‑68 (MQ=255)
      aaGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATc                                                    >  1:3512319/1‑68 (MQ=255)
       agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:3031939/68‑1 (MQ=255)
       agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:1261374/68‑1 (MQ=255)
       agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:4147609/68‑1 (MQ=255)
       agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:4106016/68‑1 (MQ=255)
       agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   >  1:2536701/1‑68 (MQ=255)
        gagaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGt                                                  <  1:4243035/68‑1 (MQ=255)
                cGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACCTTg                                          >  1:3867035/1‑68 (MQ=255)
                 gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAttct                                         >  1:696049/1‑66 (MQ=255)
                 gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg                                         >  1:1289935/1‑68 (MQ=255)
                 gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg                                         >  1:1494148/1‑68 (MQ=255)
                  cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt                                        <  1:582275/68‑1 (MQ=255)
                  cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt                                        <  1:4196903/68‑1 (MQ=255)
                   tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg                                       >  1:1848962/1‑68 (MQ=255)
                   tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg                                       >  1:1348744/1‑68 (MQ=255)
                      cATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCgg                                    <  1:3738442/68‑1 (MQ=255)
                       aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGa                                   <  1:3093803/68‑1 (MQ=255)
                               cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:1331528/1‑68 (MQ=255)
                                gTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                          <  1:634501/68‑1 (MQ=255)
                                    gTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc                      <  1:3299794/68‑1 (MQ=255)
                                       aaGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc                   <  1:3362168/68‑1 (MQ=255)
                                         gCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGTCCGGTGGCGACAt                 >  1:1769664/1‑68 (MQ=255)
                                         gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 <  1:1249760/68‑1 (MQ=255)
                                         gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 >  1:3478653/1‑68 (MQ=255)
                                         gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 >  1:2298173/1‑68 (MQ=255)
                                          cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                >  1:3880043/1‑68 (MQ=255)
                                              gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            <  1:3960427/68‑1 (MQ=255)
                                                 ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:4115686/68‑1 (MQ=255)
                                                 ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:204712/68‑1 (MQ=255)
                                                 ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:82150/68‑1 (MQ=255)
                                                 ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:2673227/68‑1 (MQ=255)
                                                 ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:1068958/68‑1 (MQ=255)
                                                  tcCAATCTGGAGGAATCTGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttgg        <  1:215965/68‑1 (MQ=255)
                                                      aaCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag    <  1:3580310/68‑1 (MQ=255)
                                                      aaCTGGAGGAATCCGGTATCGTATACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag    <  1:2509106/68‑1 (MQ=255)
                                                       aCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaagg   <  1:1237484/68‑1 (MQ=255)
                                                        cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt  >  1:2561427/1‑68 (MQ=255)
                                                        cTGGAGGAATCCGGTATCGTTTACAATGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtacggt  >  1:3292607/1‑68 (MQ=255)
                                                             |                                                              
GTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173709‑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: