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,2390GA95.7% 82.0 / ‑5.9 29intergenic (‑71/‑458)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/5);  minor base C (2/0);  total (24/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*CP0068813,725,2400GA100.0% 82.2 / NA 29intergenic (‑72/‑457)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/5);  total (24/5)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAG  >  CP006881/3725119‑3725363
                                                                                                                           ||                                                                                                                           
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACTGCCGAtt                                                                                                             >  8:170490/1‑139 (MQ=255)
ttCATAAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGATTCATAGAAAAACATCCCACCCCGTCTATTTGGCAGGCAGCCTGTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGcccagt                                                                                                             >  8:400240/1‑135 (MQ=255)
 tcatcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCTGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAttt                                                                                                            >  7:160464/1‑139 (MQ=255)
               aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCACCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCAAGCaaa                                                                                               >  7:380341/1‑138 (MQ=255)
               aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCACGCAGCCATTTTAGCACCCCGGTTTCTCAGCCGCAGCCCCCCAGAAAAGCCCGCCGCATTCTAATCTCGCaaa                                                                                               >  6:242390/1‑138 (MQ=255)
                            aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGAAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTGTTTTCACGCAAACTAGCGGCTCTAg                                                                                  >  3:175736/1‑138 (MQ=255)
                             aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCACCCCGTTTTCTCAGCCGCCCCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                <  3:337730/139‑1 (MQ=255)
                               aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCCGCCCGCCAGAAAAGACGGCCGCTTTCTTTTCTCGCAACCTAGCGGCGCTATCaa                                                                               >  8:357591/1‑138 (MQ=255)
                               aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                               >  4:361695/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                             >  8:12548/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                             >  1:277184/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                            >  3:288943/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                            >  8:342772/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTAGTTTTGGCAGGCAGCCATTTTAGCAGCCCGGTTTATCAGCCGCAGCCCCCAAGAAAAGACGGCCGATGTATTTTATCGCAAAATAGCGCCACTAGCaaaa                                                                             >  4:119252/1‑138 (MQ=255)
                                    accAACCGTTTCATATAAAAACCTCCCTCCGTTTCTTTTTGGCAGCCCGCCTTTTTAGCCGCCCGTTTCCGCCGCCCCCCCCCTAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                         <  8:29728/137‑1 (MQ=255)
                                     ccAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCTTTTTTAGCAGCCCGTTTTCTCCCCCGCAGCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                        <  8:14399/138‑1 (MQ=255)
                                       aaaCGTTTCATACAAAAACTCCCCTCCGCTTCTTTTTGGCAGGCAGCTTTTTTAGCAGCCCGTTTTCTCAGCGGCAGCCCGAAAGAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  3:361695/139‑1 (MQ=255)
                                               cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTGATTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGGCACTaag                                                               >  3:275637/1‑136 (MQ=255)
                                               cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGa                                                              >  2:14886/1‑139 (MQ=255)
                                                                 cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCCAACTACCGGCTCGAGCAAAACAACAGACACTATGATAGCCCCGGGAAAGAGaa                                            >  6:400038/1‑139 (MQ=255)
                                                                                ggcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCAAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCATTGTCACAata                             >  1:285662/1‑139 (MQ=255)
                                                                                 gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATTACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatac                            >  7:323104/1‑138 (MQ=255)
                                                                                 gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGAGGCTCTAGCAACACTCATGACTCGATTCTAGCCCGTGGAAAGCGAATTCTTTGTCACAatat                            >  6:102206/1‑139 (MQ=255)
                                                                                         ttttAGCAGCCCGTTTTCTCAGCCGCAGCCCGAAAGAAAATAAGTCCGAATTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGTAAAGAGAATTATTTGTAAAAAAATGATTGa                      >  2:62351/1‑138 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATATACCGTGGAAAGAGAATTCTTTGTCACAATAAGAGTGAAAAAAACaa             >  6:165106/1‑138 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa             >  3:375442/1‑138 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc            >  3:362589/1‑139 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCTCTAGCAAAACTCATGACTGTATGTTAGACCGTCGAATGAGAATTGTTTGTCACAATATGAGTTAAAAAACCac             >  6:294307/1‑137 (MQ=255)
                                                                                                             cAGCCGCTGCCCGCAAGAAAGGACGTCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAg  <  8:42063/139‑1 (MQ=255)
                                                                                                                           ||                                                                                                                           
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAG  >  CP006881/3725119‑3725363

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: