Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,239 2 bp→AA 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% 75.6 / NA 27intergenic (‑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 (14/13);  total (14/13)
*CP0068813,725,2400GA96.2% 80.1 / NA 27intergenic (‑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 (14/11);  minor base C (0/1);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725114‑3725369
                                                                                                                                ||                                                                                                                                 
cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                                                        >  1:95199/1‑139 (MQ=255)
     ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGaat                                                                                                                   >  3:292157/1‑137 (MQ=255)
        atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttct                                                                                                                >  8:165188/1‑139 (MQ=255)
                       cacaCCCGTAAATATACAAGAAAACGTTAAATACAAAAACACCCCTCCGTTTCTTTTGGGCCGGCAGCCTTTTTCGCAGCCGGTTTTCTCCGCCCCCGCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                                 <  2:95199/139‑1 (MQ=255)
                             cGTAACTATACAACCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTCGCCCCCCGTTTTCTCAGCCGCACCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                           <  6:217877/139‑1 (MQ=255)
                               tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                           >  4:286810/1‑137 (MQ=255)
                                      aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  8:124825/1‑139 (MQ=255)
                                      aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCACCCCGTTTTCTCAGCCGCATCCAGCAAGATAAGACGGCCGAGTTATTTTATAGACAACATACGGCTCTAGCaaaa                                                                                   >  6:191661/1‑138 (MQ=255)
                                      aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCGTTTCTCGCAAAGTAGCGGCGCTAGCAAAAc                                                                                  >  2:50534/1‑139 (MQ=255)
                                        aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACt                                                                                 >  8:87047/1‑138 (MQ=255)
                                          gCAACCGTTTCAAACAAAAACACCCCTCCGCTTCTTTTTGCCAGCCAGCCTTTTTCGCAGCCCGTTTACACAGCCGCAGCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  6:99270/139‑1 (MQ=255)
                                          gCAAACGTTTCATAGAAAAACATCCCGCCGCTTCTTTTTGGCAGGCTGCCTTTTTAGCAGCCGGTTTTGTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  7:40176/139‑1 (MQ=255)
                                          gCAAACGTTTCATACAAAAACAACCCTCCGCTTCTTTTTGGCCGCCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  1:44230/139‑1 (MQ=255)
                                            aaccttttCATACAAAAACATCCCGCCGCTTTTTTTTGGCTGGCCGCTTTTTTCGCCGCCCGTTTTCACACCCCCCCCCCCCACCAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                            <  6:286813/134‑1 (MQ=255)
                                              ccttttAATAAAAAAACATCCCTCCGATTTTTTTTGGCCGGCCGCTTTTTACGCCGCCGTTTTTCCCCCCGGCCCCCCTCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  4:292157/136‑1 (MQ=255)
                                              acttttCAAACAAAAACACCCCACCGCTTCTTTTTGGCCGGCAGCATTTTTACCAGCCGGTTTTCACCGCCGCACCCCCAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  1:220456/136‑1 (MQ=255)
                                              acttttAATAGAAAAACACCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCTTTTTCTCAGCGCCAGCCCGAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  6:123105/136‑1 (MQ=255)
                                              aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCATTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  5:285168/139‑1 (MQ=255)
                                              aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTACCGGCACTAGCAAAACAAATGAct                                                                          >  8:40176/1‑139 (MQ=255)
                                              aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTATCAAAACTCAAGAc                                                                           >  3:189558/1‑138 (MQ=255)
                                              aCGTTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGGCAGCCTTTTTAGCACCCCGTTTTCTCAGCCGCAGCCCCCAAGAAAAGACCGCCGAATTCATTTCTCGCCAGCTAGAGGCTCATGCAAAACTCaacac                                                                           >  5:4356/1‑134 (MQ=255)
                                              aCGTTTAATACAAAAACAGCCCACCGTTTCTTTTGGGCAGGCAGCCTTTTTCGCAGCCCTTTTTCTCAGCCCCCGCCCCCAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  6:118738/139‑1 (MQ=255)
                                                                               ttggacggcagCCTTTTTAGCCGCGGGTTTGCGCAGCCGCTGCACGCATCAAAAGCCGGCCGATTTCTTTTCTCTCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                          <  1:6234/132‑1 (MQ=255)
                                                                                                      agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTAGTTTCTCGCAAACAAGCGGCTCAAGCAAAACGCATGACTGTATGAGAACCCCTGGAAAGAGAATACATTGTCAAAATATGATTGAAAAAAca                     >  2:205629/1‑137 (MQ=255)
                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                  >  4:144480/1‑139 (MQ=255)
                                                                                                                  cAGCCGCCGCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAg        <  7:190804/139‑1 (MQ=255)
                                                                                                                        cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGCCCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGGTTTTAGAAGttt  >  8:178891/1‑139 (MQ=255)
                                                                                                                                ||                                                                                                                                 
CCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725114‑3725369

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: