Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,147,332 A→T intergenic (+584/‑437) RPPX_18425 → / → RPPX_18430 enoyl‑CoA hydratase/conjugal transfer protein TraR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,147,3320AT88.2% 16.2 / 0.9 17intergenic (+584/‑437)RPPX_18425/RPPX_18430enoyl‑CoA hydratase/conjugal transfer protein TraR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (7/8);  total (8/9)
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.68e-01

CGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAG  >  CP009974/4147230‑4147480
                                                                                                      |                                                                                                                                                    
cGCTCTTATGGTGGGTTACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTgggg                                                                                                      <  1:1744222/151‑1 (MQ=9)
                                                  aaaCGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCCTCCTGGGGCCTCACCGCTTGACGGGCATCCCCGCCCGTCACCTCCCCCCCCAACGCCt                                                    >  2:4387382/1‑151 (MQ=11)
                                                                                     cgctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatca                <  1:3074906/151‑1 (MQ=11)
                                                                                     cgctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatca                <  1:4979131/151‑1 (MQ=11)
                                                                                      gctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaa               >  2:1417383/1‑151 (MQ=11)
                                                                                        tCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaaga             <  1:3899788/151‑1 (MQ=11)
                                                                                         ccGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagat            <  1:3102762/151‑1 (MQ=11)
                                                                                         ccGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagat            >  1:330475/1‑151 (MQ=11)
                                                                                          cGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc           <  1:1312301/151‑1 (MQ=11)
                                                                                           ggggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatca          <  2:4085320/151‑1 (MQ=11)
                                                                                            gggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaa         <  1:235622/151‑1 (MQ=11)
                                                                                            gggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCCGAAGTCGCGaacatcaagatcaagatcaa         >  2:4069032/1‑151 (MQ=11)
                                                                                             ggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaag        >  2:4057512/1‑151 (MQ=11)
                                                                                              gTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaaga       <  1:3052413/151‑1 (MQ=11)
                                                                                               tCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagat      >  1:4139638/1‑151 (MQ=11)
                                                                                                  ccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCCTGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCCCCTGAAGTCGCGaagatcaagatcaagatcaagatcaa   >  1:4640418/1‑151 (MQ=11)
                                                                                                  ccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGcagatcaagatcaagctcaagatcaa   >  2:716293/1‑151 (MQ=11)
                                                                                                   ctctctCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaag  <  2:3300254/151‑1 (MQ=14)
                                                                                                   ctctctCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaag  <  1:2675362/151‑1 (MQ=14)
                                                                                                      |                                                                                                                                                    
CGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAG  >  CP009974/4147230‑4147480

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