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.6% 238.9 / ‑4.8 71D866E (GAT→GAGrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (22/48);  total (23/48)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173707‑4173832
                                                               |                                                              
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGTGTGAAGCTGCGCTCTCCATACTGGAGGAAt                                                            <  1:3323028/68‑1 (MQ=255)
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:94591/68‑1 (MQ=255)
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:3093488/68‑1 (MQ=255)
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:4786755/68‑1 (MQ=255)
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:2404707/68‑1 (MQ=255)
gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt                                                            <  1:1931244/68‑1 (MQ=255)
 ggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATc                                                           >  1:2686446/1‑68 (MQ=255)
   tCCGGAAGAGATCACCGCTGACATCACGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCg                                                         >  1:3271963/1‑68 (MQ=255)
    ccGGATGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCgg                                                        >  1:759941/1‑68 (MQ=255)
       gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt                                                     <  1:2161968/68‑1 (MQ=255)
       gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt                                                     <  1:4082576/68‑1 (MQ=255)
        aaGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATc                                                    <  1:4559298/68‑1 (MQ=255)
         tgagaTCACCGCTGACATCCCGAACGTGTGGGACGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:1889708/67‑1 (MQ=255)
         agagagCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:4879867/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCTAACTGGAGGAATCCGGTATCg                                                   <  1:3395348/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:5026691/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:3629094/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:965464/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   >  1:3230234/1‑68 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:144319/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:4248527/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:5341786/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:3676159/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:4881260/68‑1 (MQ=255)
         agagATCACCGCTGACATCCCGAACGTGGGTGAACCTGCGCTCTCCAAACTGGAGGAATCCGGTAACg                                                   >  1:4423875/1‑68 (MQ=255)
         agagACCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg                                                   <  1:2870542/68‑1 (MQ=255)
          gagaTCACCGCTGACATCCCGAACGTGGGGGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGt                                                  <  1:1821172/68‑1 (MQ=255)
            gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt                                                <  1:395714/68‑1 (MQ=255)
            gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt                                                <  1:4798794/68‑1 (MQ=255)
              tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc                                              >  1:2712590/1‑68 (MQ=255)
              tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc                                              <  1:1114053/68‑1 (MQ=255)
                aCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAt                                            >  1:497687/1‑68 (MQ=255)
                 ccGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACAtt                                           >  1:995782/1‑68 (MQ=255)
                    cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt                                        <  1:6004216/68‑1 (MQ=255)
                    cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt                                        <  1:314540/68‑1 (MQ=255)
                    cTGACATCCCGAACGTGGGTGAAGCTGCGCGCTCCAAACTGGAGGAATCCGGAATCGTTTACATTGGt                                        <  1:5170596/68‑1 (MQ=255)
                        cATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGTTATCGTTTACAATGGTGCgg                                    >  1:2965980/1‑68 (MQ=255)
                         aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGa                                   <  1:3731601/68‑1 (MQ=255)
                           cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 >  1:347482/1‑68 (MQ=255)
                           cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 >  1:1334136/1‑68 (MQ=255)
                           cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg                                 <  1:2510171/68‑1 (MQ=255)
                                 cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:3434169/1‑68 (MQ=255)
                                 cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                           >  1:4139423/1‑68 (MQ=255)
                                 cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTAACATTGGAGCGGCAGTCACCg                           >  1:5546637/1‑68 (MQ=255)
                                  gTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                          <  1:102100/68‑1 (MQ=255)
                                      gTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc                      <  1:5950480/68‑1 (MQ=255)
                                         aaGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc                   <  1:4018937/68‑1 (MQ=255)
                                           gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 >  1:820048/1‑68 (MQ=255)
                                           gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                 >  1:1980102/1‑68 (MQ=255)
                                            cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                >  1:2880835/1‑68 (MQ=255)
                                            cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                >  1:5199083/1‑68 (MQ=255)
                                             tGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACTTTc               >  1:4300894/1‑68 (MQ=255)
                                              gcgcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCt              <  1:5116580/68‑1 (MQ=255)
                                                gcTCTCCATACTGGAGGAATCCGGGATCGTTTACATTGGGGCGGAACTGCCCGGTGGCGCCATTCtgg            >  1:3861923/1‑68 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            >  1:1512342/1‑68 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            <  1:3171602/68‑1 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            >  1:5335292/1‑68 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            <  1:1532723/68‑1 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg            <  1:4799727/68‑1 (MQ=255)
                                                gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGGGACCGGTGGCGACATTCtgg            <  1:2932905/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:5745520/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:1312714/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:5946908/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:2009909/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:6157438/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg         <  1:1339954/68‑1 (MQ=255)
                                                   ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGCGCGGAAGAGACCGGTGGCGACATTCtggttg         <  1:6155997/68‑1 (MQ=255)
                                                     ccAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt       <  1:5612814/68‑1 (MQ=255)
                                                     ccAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt       <  1:5776533/68‑1 (MQ=255)
                                                         aCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaagg   <  1:1980452/68‑1 (MQ=255)
                                                          cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt  >  1:349201/1‑68 (MQ=255)
                                                               |                                                              
GGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT  >  CP009273/4173707‑4173832

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