Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,237 A→G 100% 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,2370AG86.5% 70.4 / ‑2.5 34intergenic (‑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/7);  minor base C (0/5);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 2.85e-03
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.

CAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA  >  CP006881/3725124‑3725362
                                                                                                                      |                                                                                                                             
caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAtttcttt                                                                                                       >  4:36961/1‑139 (MQ=255)
caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATACCTCCGCTTCTTTTTGGCAGGCAG‑C‑CATTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGatttattt                                                                                                       >  8:27984/1‑139 (MQ=255)
     aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                                                                  >  4:77854/1‑139 (MQ=255)
       aCCATCCACACCCGAACAGATAAACGCAAACTTTTCAAACAAAACCAGCCCTCCGCTTCTTTTTGGCAGGCAG‑CGT‑TTTT‑A‑GCCGCCCGTTTCCGCCGCCGCCGCCCGAAAGAAACGACGGCCGATTTCTTTTCTCGCa                                                                                                <  7:243976/139‑1 (MQ=255)
           tCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAG‑C‑CGTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGACGAATTCTTATCTCGCaa                                                                                               >  3:33650/1‑136 (MQ=255)
                     tAAAGATACAAGCAAACGTTTCATAGAACAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGGCCGTTTTCTCCGCCGCAGCCCGCCAGAAAAGACGGCCGATTGCTTTTCTCGCAAACTAGCGGcttgt                                                                                  >  2:268618/1‑136 (MQ=255)
                     tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGATAAGACGGCCGAGTTCTTTTCTGGCAACCTAGCGGctctc                                                                                  >  2:138322/1‑138 (MQ=255)
                     tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                  >  4:330537/1‑139 (MQ=255)
                        aGATACAAGCAATTGTTTCATATAAAAAAATCCCTCCTGTGTTTTTTGTGAGTCAGAG‑T‑TTTT‑C‑GCTGCCTGGTTTCTCCCCCGCCCCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                               <  6:254091/139‑1 (MQ=255)
                            aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                           >  2:12730/1‑139 (MQ=255)
                            aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTCTTTTCTCGCAACATAGCGGCTCTAGCaaaa                                                                            >  5:286888/1‑138 (MQ=255)
                            aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCTCAAACTAGCGGCTCTAGCaaaa                                                                            >  8:361302/1‑138 (MQ=255)
                            aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTAGCaaaa                                                                            >  5:254091/1‑138 (MQ=255)
                            aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTATCGCAAACTAGCCGCGCTAGAAAAAc                                                                           >  1:187479/1‑139 (MQ=255)
                                gCACACGTTTCATACAAAAACATCCCTCCGCTTGTTTTTGGCTGGCAG‑C‑CTTTTT‑C‑GCAGCCGGTTTTCTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                       <  6:286888/139‑1 (MQ=255)
                                gCAACCGTTGCAAACAACAACATCCCGCCGCTGATTTTGGGCAGGCCG‑C‑CTTTTT‑C‑GCCGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                       <  7:72748/139‑1 (MQ=255)
                                caaaccGTTTCATAGAAAAACATCCATCCGCTTCTTTTTGCCAGGAAG‑CAT‑TTTT‑A‑GCACCCCTTTTTCTCAGCCGCCCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                       <  7:361302/134‑1 (MQ=255)
                                    ccGTTCCATAAAAAACCACCCAGCCGCTTCTTTTGGCCTGCCCG‑C‑CTTTTT‑A‑GCCGCCCGTTTTCCCCGCCGCAGCCCCCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  1:64467/138‑1 (MQ=255)
                                    atttttCATATAAAAACATCGCTGCGCTTTTTTTTGCCTGTCCG‑CAT‑TTTT‑A‑GCAGCCGGTTTTCGCAGCCGCACCCCCAAAAAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  8:69013/136‑1 (MQ=255)
                                    aCGTTTCATAGAAAAACATCCCTCCGGTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  3:157955/139‑1 (MQ=255)
                                    aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑CAT‑TTTT‑A‑GCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAAGACGGACGAATTATTTTATAGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   >  5:131165/1‑139 (MQ=255)
                                    aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAtcccg                                                                   >  8:256491/1‑135 (MQ=255)
                                    aCGTTTAATACAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑CAT‑TTTT‑A‑GCAGCCCGTTTTCCCCGCCGCAGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  3:330537/139‑1 (MQ=255)
                                    aCGTTAAATATAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCG‑C‑C‑TTTTAA‑GCACCCCGTTTTCTCCGCCCCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  3:77854/139‑1 (MQ=255)
                                      gTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActct                                                                 <  8:231783/139‑1 (MQ=255)
                                                 aaaCATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGAAGGCCGATTTAGTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAATCTATGATAGACCg                                                      >  4:220986/1‑139 (MQ=255)
                                                              cttcttTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTCGCAAACTAGCGGCTCGAACAAAACTCATTACTCTATGATATAACGTGGAAAGAGAAtt                                         >  3:33326/1‑139 (MQ=255)
                                                                 cttTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCtt                                      >  6:178733/1‑139 (MQ=255)
                                                                     ttGGCCGGCAG‑C‑C‑TTTT‑ATCCGGCCCTTTTGCTCAGCCGCAGCCCCCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                   <  7:27984/138‑1 (MQ=255)
                                                                                     tttt‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGa                     >  8:110719/1‑138 (MQ=255)
                                                                                               gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa            >  8:56151/1‑139 (MQ=255)
                                                                                                cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc           >  1:42027/1‑139 (MQ=255)
                                                                                                cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCGTTTGTCGCAAACAACCCGCGCGAGCAAAACTCCTGCCGCTATTATAGACCGTGGAAAGAGAAGTCGTTTGCCCAATATGAGTGTAAAAACCAAc           >  1:173964/1‑139 (MQ=255)
                                                                                                          cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACAACCGGCTCTAGCAAAACACATTACTCGATGATAGACCGTGGAAAGAGAATTCTTTGTCACAAAATGATTGTAAAAACCAAACAATTTTTa  >  2:135032/1‑138 (MQ=255)
                                                                                                                      |                                                                                                                             
CAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA  >  CP006881/3725124‑3725362

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: ‑