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,3320AT81.0% 13.8 / 6.4 22intergenic (+584/‑437)RPPX_18425/RPPX_18430enoyl‑CoA hydratase/conjugal transfer protein TraR
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (5/12);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAG  >  CP009974/4147203‑4147480
                                                                                                                                 |                                                                                                                                                    
cgcGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTACCCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTACGCTCCGCTCCGGGGTCCCCTCCCTACCGTCTTGCTCCCGGCAg                                                                                                                                 >  2:4330472/3‑151 (MQ=1)
                                 tgcggctgGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAg                                                                                                <  1:1145241/149‑1 (MQ=14)
                                   cggcggGCCACTTTCTGTCAACCGCTACAAAGAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAgcg                                                                                              <  1:3792369/151‑1 (MQ=255)
                                                                       aCCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCAGCGTGGGGCCTCAGCGCTTGACGGGCAGCCATGCCCGTCACCTCCTTCCGCa                                                          >  1:1722551/1‑151 (MQ=9)
                                                                          aaaagaCGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACCGGCATCCATGCCCGTCACCTCCGTCCGCAAgg                                                       >  2:5984348/1‑151 (MQ=9)
                                                                                                               gcgctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                 <  2:5022751/151‑1 (MQ=11)
                                                                                                               gcgctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                 <  1:6372036/151‑1 (MQ=11)
                                                                                                                  ctCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaag              <  2:978535/151‑1 (MQ=11)
                                                                                                                   tCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaaga             >  1:4724072/1‑151 (MQ=11)
                                                                                                                   tCCGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaaga             <  1:2924348/151‑1 (MQ=11)
                                                                                                                     cGGGGTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc           <  1:2198158/151‑1 (MQ=11)
                                                                                                                      ggggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatca          <  2:3150480/151‑1 (MQ=11)
                                                                                                                      ggggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatca          <  2:4748457/151‑1 (MQ=11)
                                                                                                                       gggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaa         >  1:509570/1‑151 (MQ=11)
                                                                                                                       gggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCCCCTGAAGTCGCGaagatcaagatcaagatcaa         >  1:3368341/1‑151 (MQ=11)
                                                                                                                        ggTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaag        <  1:6215650/151‑1 (MQ=11)
                                                                                                                         gTCCCCTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGacgatcaagatcaagatcaaga       >  1:6024456/1‑151 (MQ=11)
                                                                                                                           ccccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGTAGTCGCGaagatcaagatcaacatcacgatc     >  2:5677355/1‑151 (MQ=11)
                                                                                                                           ccccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatc     <  1:2746657/151‑1 (MQ=14)
                                                                                                                            cccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatca    <  2:2662294/151‑1 (MQ=14)
                                                                                                                             ccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaa   <  1:4511196/151‑1 (MQ=11)
                                                                                                                             ccTCTCTCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaa   <  2:1799094/151‑1 (MQ=14)
                                                                                                                              ctctctCCGGCCTCGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaag  <  2:5230052/151‑1 (MQ=11)
                                                                                                                                 |                                                                                                                                                    
CACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGGCCTACGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAG  >  CP009974/4147203‑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: