Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,433,439 C→A R317L (CGC→CTC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,433,4390CA95.2% 123.8 / ‑1.3 42R317L (CGC→CTC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (18/22);  total (19/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 = 8.59e-01

AAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACC  >  CP009273/3433377‑3433503
                                                              |                                                                
aaaaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATg                                                             >  1:4653563/1‑68 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:3447597/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:3961128/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:2185077/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:2648274/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:2390560/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGcc                                                           <  1:2455634/68‑1 (MQ=255)
  aaCCTTAACCTGTGATCCGGTTACTAGTCAGCGATTCTTGCCGGTGGCCAGTTTTCAAGGAGCAttac                                                           <  1:4353982/68‑4 (MQ=255)
   aCCTTAACCTGTGATCCGGTTACTCGTCAGCGATCCTTGCCTGTGGCCAGTTTTCCAGGAGTATGccc                                                          >  1:3720930/1‑68 (MQ=255)
   aCCTTAACCTGTGATCCGGTGACGCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGccc                                                          >  1:3431603/1‑68 (MQ=255)
          ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGCCa                                                   >  1:1550670/1‑68 (MQ=255)
          ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACa                                                   >  1:4167235/1‑68 (MQ=255)
          ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACa                                                   >  1:4540949/1‑68 (MQ=255)
          ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACa                                                   >  1:3942488/1‑68 (MQ=255)
           cTGTGATCCGGTTACTCGTCAGCGATGCTTGCTGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAg                                                  >  1:1764907/1‑68 (MQ=255)
            tctgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGt                                                 <  1:4380334/66‑1 (MQ=255)
                 tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACg                                            <  1:2660489/68‑1 (MQ=255)
                 tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACg                                            <  1:2939412/68‑1 (MQ=255)
                 tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCAACg                                            >  1:3610547/1‑68 (MQ=255)
                           cGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGc                                  >  1:1611234/1‑68 (MQ=255)
                               aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGt                              <  1:3910800/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:432823/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:959124/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:1599893/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:2360470/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:2577457/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:2604116/68‑1 (MQ=255)
                                gCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTc                             <  1:3064495/68‑1 (MQ=255)
                                  gATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCtt                           >  1:1081433/1‑68 (MQ=255)
                                  gATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCtt                           >  1:3349733/1‑68 (MQ=255)
                                          ccGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                   >  1:3968786/1‑68 (MQ=255)
                                          ccGGTGGCCAGTTTACCAGGAGCATGCCCAGAGACAGTCCACCGGAAGCCAGCACGTCTTTAATCTCa                   >  1:3878405/1‑68 (MQ=255)
                                            ggTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGt                 <  1:1882578/68‑1 (MQ=255)
                                            ggTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGt                 <  1:937350/68‑1 (MQ=255)
                                            ggTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGt                 <  1:3748173/68‑1 (MQ=255)
                                                   aGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAt          <  1:4528295/68‑1 (MQ=255)
                                                    gTTTTCCAGGAGCATGCCCAGAGACAGTCCATGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtt         <  1:4170524/68‑1 (MQ=255)
                                                           aGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc  >  1:4022776/1‑68 (MQ=255)
                                                           aGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc  >  1:3632049/1‑68 (MQ=255)
                                                           aGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc  >  1:3638493/1‑68 (MQ=255)
                                                           aGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc  >  1:4467430/1‑68 (MQ=255)
                                                           aGGAGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc  >  1:763812/1‑68 (MQ=255)
                                                              |                                                                
AAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACC  >  CP009273/3433377‑3433503

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: