Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,237 A→G intergenic (‑69/‑460) 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,2370AG87.1% 72.8 / ‑3.9 32intergenic (‑69/‑460)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/5);  minor base C (1/3);  total (23/9)
Fisher's exact test for biased strand distribution p-value = 4.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725119‑3725369
                                                                                                                         |                                                                                                                                    
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                   >  3:160165/1‑139 (MQ=255)
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTGTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGGTTTCACACCCCCAGCCCCCAAGAAAAGACGGCCGCtt                                                                                                                   >  3:227820/1‑139 (MQ=255)
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTTGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACCGCCGAtg                                                                                                                   >  1:131064/1‑138 (MQ=255)
     caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCATCCGCAGCCCGCAAGAAAAGACGGCCGAgttcgtt                                                                                                              >  7:305466/1‑139 (MQ=255)
                      cccGTAACTATACAACCAAACGTTACATCAAAAAACATCCCTCGGTTTATTTTTGGCCGGAAGCTTTTTTAGCAGCCTTGTTTCTCTGCCGCCGCCCCCATAAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                             <  4:324146/139‑1 (MQ=255)
                        cGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAAATAGCGGctc                                                                                           >  7:290085/1‑139 (MQ=255)
                        cGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                           >  5:61976/1‑139 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGctc                                                                                           >  4:232437/1‑137 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTCGCAAACTAGCGGCGCTGGCAAAAc                                                                                  >  1:8230/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTACTTTTCTCGCAAACTAGCGGCTCTACCaaaa                                                                                   >  2:207806/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  7:82277/1‑139 (MQ=255)
                                   aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGTCCGATTTCTTTTGTCGCAACCTAACGTCTCGAGCAAAACt                                                                                 >  8:254030/1‑138 (MQ=255)
                                     gaaaaCTTTTAATAGAAAAAAATCCCTCCGCTTCTTTTTGGATGGCAGCTTTTTTAGCCGCCGGTTTTCTCACCCGCAGCCCTAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  8:305466/137‑1 (MQ=255)
                                     ccacaCGTTTAATATAAAAACATCCCGCCGATGATTTTTGGCCGCCCGCTTTTTTACCCCCCGGTTTTCTCCCCCCCCGCCCCCAACAAACGACGGCCGATTTCTTTTCTTGAAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  8:3009/135‑1 (MQ=255)
                                       aaaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                            <  7:188931/139‑1 (MQ=255)
                                         cggTTTCATACAAAACAGACCGTCCGCTGTGTTTTGGCCGGCCGCCTGTTTCGCCGCGCGTTTCCTCAGCCGCCGCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTCGCGGCTCTAGCAAAACTCATGAct                                                                          <  8:34469/137‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  3:338233/139‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTTGCAAACAAGCGGCTCGAGCAAAACTCATGAct                                                                          >  7:3009/1‑139 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          >  4:148668/1‑139 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTGCTTTTGGGCAGGCACCCTTTTTCGCAGCCCGTTTTCTCAGCCGCAGCCCCCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  8:197917/139‑1 (MQ=255)
                                                    aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAAAGAc                                                               >  5:139164/1‑139 (MQ=255)
                                                           tCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGa                                                         >  5:84248/1‑138 (MQ=255)
                                                              cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAATCTATGAAAGACCGTTGaaa                                                       >  3:62698/1‑137 (MQ=255)
                                                                 cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCACCCGCCCCCCGAAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa                                                  <  7:73015/139‑1 (MQ=255)
                                                                      cttTTTGCCCGGCAGCTTTTTTAGCAGCCGGTTTTCACAGCCCCACCCCGCAACAAACCACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                              <  6:61976/138‑1 (MQ=255)
                                                                        ttttGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTg                                           >  7:250034/1‑139 (MQ=255)
                                                                                                  gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCATTGTCACAATATGAGTGAAAAAACCaa                   >  5:60398/1‑139 (MQ=255)
                                                                                                  gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTCGCAAAATATCGGCTCTAGCCAAACTCAAGACTCTATGAAAGACCGTGGAAAGAGAATTATTTGTCACAATATTAGTGAAAAAACCaa                   >  2:267940/1‑139 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTTAAAAAACCaa                   >  5:92890/1‑138 (MQ=255)
                                                                                                       ttttCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGAGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAACAACAAGt              >  1:23518/1‑139 (MQ=255)
                                                                                                                ccGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTGTCGCAAACTAGCGGCTCTAGCCAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATTAGTGCAAAAACCCACAAGTTTTTATAAg     >  6:186156/1‑139 (MQ=255)
                                                                                                                   cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGATTGAAAAAACCAACAAGTTTTTTGaatgtt  >  4:3541/1‑135 (MQ=255)
                                                                                                                         |                                                                                                                                    
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725119‑3725369

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