Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,172,812 T→G V547G (GTT→GGT)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,172,8120TG98.3% 220.1 / ‑3.6 59V547G (GTT→GGT) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (25/33);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCT  >  NZ_CP009273/4172726‑4172896
                                                                                      |                                                                                    
caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCg                                                                                   <  1:703297/90‑1 (MQ=255)
caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCg                                                                                   <  1:466957/90‑1 (MQ=255)
 aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCga                                                                                  <  2:955104/90‑1 (MQ=255)
 aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCga                                                                                  <  1:147015/90‑1 (MQ=255)
 aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCga                                                                                  <  2:133140/90‑1 (MQ=255)
 aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCga                                                                                  <  1:543619/90‑1 (MQ=255)
  acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCgag                                                                                 <  1:175254/90‑1 (MQ=255)
  acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCgag                                                                                 >  2:701380/1‑90 (MQ=255)
    ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGAc                                                                               >  2:171066/1‑90 (MQ=255)
       cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTa                                                                            >  1:118609/1‑90 (MQ=255)
       cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTa                                                                            >  2:532564/1‑90 (MQ=255)
       cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTa                                                                            >  2:274507/1‑90 (MQ=255)
       cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTa                                                                            >  1:1016184/1‑90 (MQ=255)
          tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTacac                                                                         <  2:607828/90‑1 (MQ=255)
           cTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACAcc                                                                        >  1:947497/1‑90 (MQ=255)
            tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACGccc                                                                       <  2:118609/90‑1 (MQ=255)
            tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACAccc                                                                       <  2:130337/90‑1 (MQ=255)
               gaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    <  1:831192/90‑1 (MQ=255)
               gaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    <  1:746234/90‑1 (MQ=255)
                   aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc                                                                >  2:749656/1‑90 (MQ=255)
                   aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCAc                                                                >  2:235933/1‑90 (MQ=255)
                   aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCAc                                                                >  2:227832/1‑90 (MQ=255)
                   aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCAc                                                                >  1:758224/1‑90 (MQ=255)
                   aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCAc                                                                >  1:426119/1‑90 (MQ=255)
                         aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACATGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGt                                                          >  2:839631/1‑90 (MQ=255)
                          aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCATTACGGTc                                                         <  2:779843/90‑1 (MQ=255)
                          aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTc                                                         <  2:689202/90‑1 (MQ=255)
                                gtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTAt                                                   >  2:283903/1‑90 (MQ=255)
                                  aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt                                                 >  1:49708/1‑90 (MQ=255)
                                    ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc                                               <  2:577793/90‑1 (MQ=255)
                                    ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc                                               <  2:197454/90‑1 (MQ=255)
                                          aCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    <  1:779787/63‑1 (MQ=255)
                                          aCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    >  2:779787/1‑63 (MQ=255)
                                           cTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa                                        <  1:532564/90‑1 (MQ=255)
                                             cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCgaga                                                                                <  2:336952/48‑1 (MQ=255)
                                             cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCgaga                                                                                >  1:336952/1‑48 (MQ=255)
                                             cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAc                                      >  2:380031/1‑90 (MQ=255)
                                                    ggCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                               >  2:706020/1‑90 (MQ=255)
                                                    ggCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                               <  2:639478/90‑1 (MQ=255)
                                                     gCGGTCTGACCCGTGAACGTGCGGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAg                              <  2:544756/90‑1 (MQ=255)
                                                         tCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    <  1:311034/48‑1 (MQ=255)
                                                         tCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGAc                                                                    >  2:311034/1‑48 (MQ=255)
                                                         tCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                          >  2:594813/1‑90 (MQ=255)
                                                         tCTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                          >  2:548529/1‑90 (MQ=255)
                                                          cTGCCCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg                         <  1:673767/90‑1 (MQ=255)
                                                          cTGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg                         <  1:322996/90‑1 (MQ=255)
                                                           tGACCCGTGAACGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGa                        <  2:1016184/90‑1 (MQ=255)
                                                                    aaCGTGCAGGCTTCGAAGGTCGAGACGTTCACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTc               <  2:62841/90‑1 (MQ=255)
                                                                    aaCGTGCAGGCTTCGAAGGTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTc               <  2:1016429/90‑1 (MQ=255)
                                                                    aaCGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTc               <  2:435687/90‑1 (MQ=255)
                                                                    aaCGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTc               <  1:186314/90‑1 (MQ=255)
                                                                    aaCGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTc               <  1:252413/90‑1 (MQ=255)
                                                                      cGTGCAGGCTTCTAAGGTTGAGTCGTACACCCGACTCACTACGGTCGCGTATTTCCAATCTAAACCCCTGAAAGGCCGAACATCGGTCTg             <  1:432464/90‑1 (MQ=255)
                                                                      cGTGCAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTg             >  1:334741/1‑90 (MQ=255)
                                                                          cAGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa         <  1:380031/90‑1 (MQ=255)
                                                                           aGGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCaa        >  1:194795/1‑90 (MQ=255)
                                                                              cTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActc     >  1:38609/1‑90 (MQ=255)
                                                                              cTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActc     <  1:701380/90‑1 (MQ=255)
                                                                                 cGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct  >  2:327799/1‑90 (MQ=255)
                                                                                      |                                                                                    
CAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCT  >  NZ_CP009273/4172726‑4172896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: