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,2370AG81.3% 90.2 / 0.7 45intergenic (‑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 (26/10);  minor base C (2/7);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA  >  CP006881/3725119‑3725362
                                                                                                                             |                                                                                                                             
ttCATCAAAAAGACCATCCCCACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCtt                                                                                                            >  1:191822/1‑139 (MQ=255)
ttCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGATAATCATCCCTCCGCGTATTTTTGGCATGCAGCCTTTTTAGCAGCCCGTTTTCTCCTACGCAGCCCGCAAGAAAAGCCGGCCGagt                                                                                                            >  6:201568/1‑137 (MQ=255)
   atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                            >  3:99100/1‑136 (MQ=255)
     caAAAAGACAATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCACGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                       >  1:141128/1‑139 (MQ=255)
        aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGGCGAtttcttt                                                                                                       >  4:151019/1‑136 (MQ=255)
            aCCATCCACACCCTTAAAGATACAAGCAAACGTTTCATATAAAAACAGCCCTCCGATTCTTTTTGGCCGGCAGCTTTTTTCGCAGCCGGTTTTCTCAGCCGCAGCCCGCAATAAAACACGGCCGATTTCTTTTCTCGCa                                                                                                <  5:198719/139‑1 (MQ=255)
               aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCaaa                                                                                              >  7:413370/1‑138 (MQ=255)
                  cacaCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                          >  8:29745/1‑139 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCTGctc                                                                                    >  2:243097/1‑137 (MQ=255)
                             acaTACAAGCAAACGTTCAAAACAAAAACAACCCACCGCTTATTTTTGGCCGGCACCCTTTTAACCACCCCGTTTTCGCCCCCCCAGCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                  <  3:409649/134‑1 (MQ=255)
                             acaTACAACCAAACGTTTAATATAAAAAAACCCCCCCGCTGATTTTTGGCTGGCAGCTTTTTTCGCCGCCGGTTTTCTCCGCCCCCCCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                  <  3:152516/134‑1 (MQ=255)
                             aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                               <  7:41452/139‑1 (MQ=255)
                             aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                               <  5:90137/139‑1 (MQ=255)
                             aGAAACAAGCAAACGTTTCATACAAAAACATCCCTCCGCTGCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCTTTTTCTCAGCCGCAGCCCGCATTAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                               <  6:349639/139‑1 (MQ=255)
                              ctacaaaGAAAACGTTTCATACAAAAACATCCCTCCGCTTATTTTTGCCCGGCAGCCTTTTTAGCAGCCCGTTTGCACAGCCGCCCCCCCCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                              <  1:243097/134‑1 (MQ=255)
                               aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTATCTCAGCCGCACCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                              >  8:13318/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCGGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                            >  8:317978/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                           >  8:332055/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGAAGCCCGTTTTCTCAGCCGCAGCCCCCAAGAAAAGACGGCCGCTTTCTTTTCACGCAAAAAAGCGGCTCAAGCaaaa                                                                            >  5:334611/1‑138 (MQ=255)
                                    aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTATTTGGCACGCAGCGTGTTTAGCACCCCGGTTTCTCATCCGCAGCCCCCAAGAAAAGACGGCCGCTTTCGTTTCACGCAAACAACGCGCTCTAGCaaacc                                                                           >  4:136995/1‑134 (MQ=255)
                                     gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCAGCCCGAAAGAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                       <  7:249352/139‑1 (MQ=255)
                                     caaaaCGTTCCATACAAAACAATCCCCCCGCTTTTTTTTGGCCGGCAGCCTTTTTAGCCGCCCGTTTTCCCCGCCGCCCCCCCCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                       <  1:65486/137‑1 (MQ=255)
                                      cAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTATCGGCACTAGCaaaaat                                                                          >  4:371284/1‑133 (MQ=255)
                                       aaaCGTTTCAAAAAAAAACTTCCCCCCGCTTATTTTTGCCTGGCAGCCTTTTTACCCGCCCTTTTTCTCCGCCCCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGAAAACTCGCGGCTCTAGCAAAACTCa                                                                        <  4:212341/136‑1 (MQ=255)
                                        aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                                                                    >  6:192692/1‑139 (MQ=255)
                                         cgttttCTTACAAAAACACCCCTCGCTTTTTTTTGGCCCGCCAGCTTTTTTCCCAGCCCGTTGTCCCCCCCCCCCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                      <  4:149152/133‑1 (MQ=255)
                                         acttttAATACAAAACCATCCCCCCGCTTTTTTTTGGCGGGCAGCCTTTTTCGCAGCCCGTTTTCTCCGCCCCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   <  7:317978/136‑1 (MQ=255)
                                         actttacATAAAAAAACACCCCTCGCTTTCTTTTTGGCCGCCCGGTTTTTTCGCCGCCCTTTTTCTCCCCCCCCCCCCCCAACAAACCCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                      <  4:145132/130‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                   >  1:19484/1‑139 (MQ=255)
                                           gTTTCATATAAAAACAGCCCGCCGTTGATTTTTGGCCGGCCGCGTTTATAGCCGCCCGTTTTCTCCGCCGCCGCCCGAACAAAAACACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                                                                    <  4:331593/136‑1 (MQ=255)
                                               cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACACTAtta                                                             >  8:181624/1‑137 (MQ=255)
                                                      aaaCATCCCTCCGCTTCTTTTTGGACGGCAGCCGTTTTAGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATAGCAAACTAGCGGCTATATCAAAACAAATTCCTAGATGATAGAACg                                                      >  6:225287/1‑139 (MQ=255)
                                                            ccccccGATTTTTTTTGCCCGGCAGCTTTTTTAGCAGCCGGTTTTCCCAGCCCCCGCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg                                                   <  3:371284/136‑1 (MQ=255)
                                                                 cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa                                           >  7:203556/1‑139 (MQ=255)
                                                                      atttttGGCCGGCAGCCTTTTTAGCAGCCGTTTTTCTCAGCCCCCGCCCGAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                       <  8:203556/137‑1 (MQ=255)
                                                                       tttttGGCAGGCAGCCTTTTTAGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCGCGAGCAAAACACATGACTAGTTGATTGACCGTGGAAATAGaaata                                        >  4:391721/1‑133 (MQ=255)
                                                                                      tgcctttttCGCAGCCCGTGCTCTCCCCCGCAGCCCCAAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGa                         <  7:13318/134‑1 (MQ=255)
                                                                                              ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACGCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATTTGAGTGaaa                   >  8:211762/1‑138 (MQ=255)
                                                                                                     agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc              >  8:11979/1‑138 (MQ=255)
                                                                                                         cGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTTTAACAATTTTAGTGAAAAAAcaca            >  4:39958/1‑133 (MQ=255)
                                                                                                          gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGATTGAAAAAACCaa            >  3:10242/1‑135 (MQ=255)
                                                                                                          gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAACaa            >  4:317559/1‑135 (MQ=255)
                                                                                                                 cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCCAACTAGCGGCTCGAGCAAAACTCATGACGCTATGATAGACCGTGGAAAAAGATTTGTTTGTAACCATATGAGTTAAAAAAACAACTAGattttt  >  2:251957/1‑137 (MQ=255)
                                                                                                                 cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtttt    >  6:296280/1‑136 (MQ=255)
                                                                                                                 cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTAGCAAACTAGCGGCTCGAGCAAAACTCATGACTCGATGATAGCCCGTGGAAATATAATTATTTGTCCAAATATGAGTTGAAAAAACAAATAGtttttt  >  2:381564/1‑137 (MQ=255)
                                                                                                                             |                                                                                                                             
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA  >  CP006881/3725119‑3725362

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: