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,2370AG84.5% 91.8 / ‑1.4 41intergenic (‑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 (23/11);  minor base C (1/6);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.41e-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.

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAA  >  CP006881/3725119‑3725352
                                                                                                                           |                                                                                                                   
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                  >  6:290660/1‑139 (MQ=255)
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCACCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                  >  5:248937/1‑139 (MQ=255)
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGGGTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCCAGAAAACCCGGCCGCtt                                                                                                  >  5:38436/1‑139 (MQ=255)
  catcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttc                                                                                                >  3:108914/1‑139 (MQ=255)
     caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCCCAAGAAAAGACGGCCGAtttcttt                                                                                             >  7:216765/1‑139 (MQ=255)
            accaccCACACCCGTAAAGATACCACCAACCTTTTCATACAAAAACACCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTCGCCGCCCGTTTTCTCCGCCGCACCCCCCAACAAACGACGGCCGATTTCTTTTCTCGCa                                                                                      <  8:126348/139‑1 (MQ=255)
               aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGCTTTCTTTTCTCggcaac                                                                                   >  4:49135/1‑134 (MQ=255)
                  cacaCCCGAAAAGATACACGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCACCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                <  4:87601/139‑1 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                          >  5:209771/1‑137 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTGTCGCAAACTAGCGGcgc                                                                          >  4:174587/1‑135 (MQ=255)
                             aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCCGCCCGTTTGCCCAGCCCCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                     <  3:267661/139‑1 (MQ=255)
                              gATACAAGCAAAGTTTTAATATAAAAACTTCCCTCCGCTTGTTTTTGCCTGGCAGCTTTTTTAGCAGCCGTTTTTCTCCGCCCCAGCCCGCAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                    <  8:216765/139‑1 (MQ=255)
                               aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATATCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                    >  7:59037/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                  >  1:153725/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                     >  8:209964/1‑135 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCGAGCaaaa                                                                  >  2:206170/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCCAACTAGCGGCTCTATCAAAAc                                                                 >  3:80157/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                     <  7:209964/135‑1 (MQ=255)
                                   aaGCAAACGTTTAATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGACGCAGCCAGCAAGAAAAGACGGCCGATTTCTTTTCTGGCAAACTAGCGGCTCGAGCaaaaat                                                                >  7:239397/1‑136 (MQ=255)
                                    aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCACCCTGTTTAGCAGCCCGTTTTCTCACCCCCAGCCCGCAAGAAAAGACGGCCGATATCTATTCTCGCAAACTATCGGCTCAAGCAAAACTCa                                                              >  8:98261/1‑139 (MQ=255)
                                     gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATAGCAAACAACCGGCTCTAGCAAAACTCAt                                                             >  2:24214/1‑139 (MQ=255)
                                     gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAAGACGGCCCAGTTCGTTTATAGCAAACTAGGCGCTGGAGCAAAACGCAt                                                             >  1:123458/1‑139 (MQ=255)
                                     ccAAACGTTACATACAAAAACATCCCTCCGCTGATTTTTGGCAGCCAGCCTTTTTAGCAGCCCTTTTTCTCAGCCCCAGCCCCCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                             <  5:141936/138‑1 (MQ=255)
                                       aaccGTTCCATATAAAAAAATCCCCCCGCTTTTTTTTGGATGCCACCCTTTTTAGACGCCCTTTTAAACCGCCGCCCCCCGCACTAAAACACGGCCGATTTCTTTTCTCGAAACCTAGCGGCTCTAGCAAAACTCATGa                                                           <  8:239397/136‑1 (MQ=255)
                                       aaaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                           <  5:140468/139‑1 (MQ=255)
                                         ccGTTTCATCCAAAAACATCCCGCCGCTTCTTTTTGGCTGGCAGCCTTTTTAGCCCCCCGTTTTCTCCGCCGCACCCCCCATTAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                         <  6:280972/138‑1 (MQ=255)
                                         acttttCAAATAAAAAAATCCCTCCGCTTATTTTTGGCCGCCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCCGCCCGAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                         <  8:160178/136‑1 (MQ=255)
                                         acttttAATAGAAAAACATCCCTCCGCTGCTTTTTGGCAGGCTGCCTTTTTCCCAGCCGGTTTTCACCGCCGCCGCCCGCAAGAACCGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                         <  3:141423/136‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACAGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                         >  3:275209/1‑139 (MQ=255)
                                         aCGTTTAATATAAAAACATCCCTCCGATTATTTTTGGCAGGATGACTTTTTAGCAGCCCTTTTGCTCAGCCGCACCCCGCAAAAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                         <  7:98261/139‑1 (MQ=255)
                                              tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTAt                                                     >  7:220237/1‑138 (MQ=255)
                                              tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCCGCTCTAGCAAAACTCAAGACTCTAt                                                     >  1:23684/1‑138 (MQ=255)
                                               cATAGAAAACCAACCCTCCGCTTCTTTTTGCCAGCCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCCCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGa                                                   <  2:9319/139‑1 (MQ=255)
                                                    aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAAACTCAGTACTCgatgatgtgt                                              >  7:251909/1‑135 (MQ=255)
                                                    aaaaaCACCCCTCCTCTTCTTTTGGGCCGGCCCCATTTTAACCCCCCCTTTTCCTCCCCCCCCCCCCCCAAAAAACCCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                                              <  4:108914/139‑1 (MQ=255)
                                                            cccTCCGCTTTTTTTTGGCCGGCCGCTTTTTTCGCAGCCGTTTTTCACCCCCGCCCCCCCAAACAACACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                                      <  2:23684/139‑1 (MQ=255)
                                                            cccTCCGCTTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCGGTTTTCTCCGCCGCAGCCCGAAAGAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                                      <  8:220237/139‑1 (MQ=255)
                                                                                 gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                             >  8:290324/1‑95 (MQ=255)
                                                                                 gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGAAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                             <  7:290324/95‑1 (MQ=255)
                                                                                                  agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACGCATGACGCTATGATAGACCGTGGAAAGAGAATTCGTTGTAACAATATGACTTaaaaaaac    >  1:255426/1‑136 (MQ=255)
                                                                                                    cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa  >  1:104886/1‑138 (MQ=255)
                                                                                                                           |                                                                                                                   
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAA  >  CP006881/3725119‑3725352

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: