Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,239 2 bp→AA 100% intergenic (‑71/‑457) U712_19205 ← / → U712_19210 Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,2390GA100.0% 47.4 / NA 18intergenic (‑71/‑458)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/5);  total (13/5)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*CP0068813,725,2400GA92.0% 42.1 / ‑5.6 18intergenic (‑72/‑457)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/3);  minor base C (0/2);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 6.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCC‑‑TTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTT  >  CP006881/3725119‑3725360
                                                                                                                            ||                                                                                                                        
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                          >  7:180560/1‑139 (MQ=255)
                      cccGAAAAGAAACCACCACACGTTAAATCCAAACACATCCCTCCCCTGCTTTGTGGCCGCCAGCTT‑‑TTTTAGCAGCCCGTTTTCGCAGCCCCAGCCCGCACCAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                    <  2:317527/139‑1 (MQ=255)
                      agcctAAACATACAAGAAAATGTTAAAAACAAAAACATCCCCCCGCTTCTTTTTGGCTGGCAGCCTTTTTCGCCGCCGTTTTGCTCAGCCGCCGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                    <  1:366089/135‑1 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCCCAAGAAAAGACCGCCGATTTATTTTATCGCAAACTAGCGGcactt                                                                                >  5:136696/1‑138 (MQ=255)
                             ctaTACAAGCAAAGCTTCAATATAACAACACCCCACGGCTTTTTTTTGCCAGCCCGCC‑‑TTTTACGCACCCCTTTTTCCCCGCCCCCCCCCCCACCACACGCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                             <  4:85857/137‑1 (MQ=255)
                             aGACACAACCAAACGTTTCATACAAAAACACCCCCCCGCTGCTTTTTGGCAGCCAGCCTTTTTAGCAGCCCGTTTCCCCCGCCCCCCCCCGCAAGAACAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                             <  7:389399/139‑1 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAACACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTGTCTCGCCAACTAGCGGCGCGAGCAAAAc                                                                         >  4:157821/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACGAGCGGCTCTAGCAAAAc                                                                         >  7:388881/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTCCTTTTCTCGCAAACTAGCGGCTCTAGAAAAAc                                                                         >  3:126720/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACAACCGGCTCTAGCaaaa                                                                          >  5:279602/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCCGCCCGCAAAAAAAGACGGCCGATTTAGTTTCTCGCAAACTACCCGCACAAGCaaaa                                                                          >  5:128703/1‑138 (MQ=255)
                                   aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCACCCCGCAAGAAAACACGGCCGATTTATATTCACGCAAACAAGCGGCGCTAGCAAAAccg                                                                       >  2:366089/1‑137 (MQ=255)
                                      aaaaCTTTTCATCTAAAAACATCCCGCGGGTTCTTTTGGCACGCAGCCCTTTTTCGCGGCCGGTTTTCTCACCCGCTCCCCGAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                    <  7:75327/138‑1 (MQ=255)
                                                                                ggcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCGCGAGCAAAACTCATGCCGGGATGAAAGCCCGTGGAAAGAGAATTGATTTTCACAata                          >  2:402160/1‑139 (MQ=255)
                                                                                                   gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGATTGAAAAAACCa           >  2:300783/1‑138 (MQ=255)
                                                                                                   gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa          >  3:42858/1‑139 (MQ=255)
                                                                                                    cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTATCGCTAACTAGCGGCTCTAGCAAAACTCCTGACTCTATGTTAGACCGTGGAAAGAGAATTGTTTTTCACAATATGAGTGTAAAAACCaa          >  4:321692/1‑138 (MQ=255)
                                                                                                              cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtttt  >  3:140606/1‑136 (MQ=255)
                                                                                                                            ||                                                                                                                        
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCC‑‑TTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTT  >  CP006881/3725119‑3725360

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