Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,239 2 bp→AA 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,2390GA100.0% 139.1 / NA 43intergenic (‑71/‑458)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/18);  total (25/18)
*CP0068813,725,2400GA97.6% 111.6 / NA 43intergenic (‑72/‑457)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/16);  minor base C (0/1);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAG  >  CP006881/3725119‑3725363
                                                                                                                             ||                                                                                                                           
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAACCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAttt                                                                                                            >  5:94324/1‑139 (MQ=255)
     caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttcttt                                                                                                        >  7:155857/1‑139 (MQ=255)
               aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCaaa                                                                                               >  7:140267/1‑138 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctg                                                                                     >  5:157389/1‑136 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGctc                                                                                     >  1:276090/1‑137 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGctc                                                                                     >  7:274320/1‑137 (MQ=255)
                          tAAAGATACAACCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTATGGCAGGCTGCCTTTTTAGCAGCCCGATATCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCGAACTAGCGGCTCTa                                                                                   >  6:234321/1‑139 (MQ=255)
                             acaTACAACAAAACGTTTAATATAAAAACATCCCTCCGTTTTTTTTGGGCCGTCTGCTTTTTTCGACGCCGTTTTTTTCCGCCCCCGCCCCAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                <  4:220454/137‑1 (MQ=255)
                             aGATACAAGCAAAGTTTTCATTGAAAAACTGCCCACCGCTTCTTTTGGCCTGGCAGCCTTTTTAGCTGCGGTTTTTCACAGCCGCCGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                <  8:54389/139‑1 (MQ=255)
                              gaaaaaaGCACACGTTTAATACAAAAACAGCCCTCCGCTTCTTTTTGCCCGCCTGCCTTTTTCGCCGCCGTTTTTCTCCGCCCCACCCCCCAAGAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                               <  6:117983/134‑1 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTGGCAAACTAGCGTCTCTAGCAAAAc                                                                            >  5:54903/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGGGCTAGCaaaa                                                                             >  2:290331/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTATCGGCTCAAGCaaaa                                                                             >  3:220454/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCATTTCTCGCAAACTAGCGGCGCTAGCaaaa                                                                             >  5:318709/1‑138 (MQ=255)
                                    aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCAGCCCCAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                         <  4:284000/139‑1 (MQ=255)
                                     gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCACCCCGCATGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                        <  7:59620/139‑1 (MQ=255)
                                     ccAAACGTTTAATATAAAAACATCCCTCCGCTTATTTTTGGCCGGAAGACTTTTTCGCAGCCCGTTTTCTACGCCGCAGCCCCAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                        <  6:157389/138‑1 (MQ=255)
                                       aaaGGTTTAATACAAAAACAGCCCCCCGCTTCTTTTGGGCCGGCCCCCTTTTTACCCCCCGGTTTTCTCCGCCCCCCCCCGAAACAACCCACGGCCGATTTCTTTTCTCGCAAACTCGCGGCTCTAGCAAAACTCATGa                                                                      <  1:316627/139‑1 (MQ=255)
                                       aaaCGTTTCATCTAAAAACAACCCTCGGTTTCTTTTTGGCCGCCTGCTTTTTTAGCCGCCCGTTTTCTCTGCCGCTGCCCGCAAGAATCGACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  4:153597/139‑1 (MQ=255)
                                       aaaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCCCAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  7:206937/139‑1 (MQ=255)
                                       aaaCGTTTAATATAAAAACATCCCACCGATGCTTTTTGGCCGCCAGCCTTTTTAGCACCCCGTTTGCTCCGCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  3:157817/139‑1 (MQ=255)
                                       aaaCGTTCAATATAAAAACAACCCTCCGATTTTTTTGGGCTGCCAGCATTTTTCCCAGCCGGTTTACACCGCCGCCGCCCGAAATAAAACACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  5:291930/139‑1 (MQ=255)
                                         ccGTTTAAAAAAAAAACATCCCTCCGTTTCTTTTGGGCTGCCCGCATTTTTAGACGCCGTTTTTACCACCCGCCGCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    <  6:195497/138‑1 (MQ=255)
                                         ccGTTACAAACAAAAACATCCCCCGGGTGCTTTTGGGCGGGCAGCGTTTTTCGCCGCCCGTTTGCCCCGCCGCCCCCCCAAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    <  1:226990/138‑1 (MQ=255)
                                         acttttCATACAAAACCACCCCGCCGCTTATTTTGGCCCGGCCGCCTTTTTCGCCGCCGGTTTCCCCCGCCGCCCCCCGCAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    <  1:42571/136‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCGAGCCAAACTCATGAc                                                                     >  7:123259/1‑138 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTCTCGCAAACTAGCCGCGCAAGCAAAACACATgcc                                                                     >  5:177339/1‑136 (MQ=255)
                                         aCGTTTCATACAAAAACATCCCTCCGATTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCCCAGCCCGAACTAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    <  8:257323/139‑1 (MQ=255)
                                              caatccaaaacCATCCGCCCGATTATTTTTGCCAGCCCTCTTTTTACGCAGCCCTTTTGCTCCGCCGACCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATg                                                               <  4:290724/133‑1 (MQ=255)
                                               cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTATCGCAAACTAGCCGCGCTAGCAAAACTCATGACGCTATGa                                                              >  4:122532/1‑139 (MQ=255)
                                                atatAAAAAAATCCGTCCGTTTTTTTTTGGATGGACGATTTTTTAGCCCCCGTTTTTCTCACCCCCAGCCCCAAAAAAACCAGGGCCGATTTCTTTTCTCTCAAAGTAGCGGCTATAGCAAAACTCATGACTCTATGAt                                                             <  8:140267/139‑1 (MQ=255)
                                                     aaaaCATCCCTCCGCTTATTTTTGGCCGGCAGCCGTTTTAGCATCCCGTTTTCTCAGCCCCACCCCGCAACAAAAGCCCGCCGATTTCTTATCGCCCAAACTAGCGGCGCTACCAAAAACCATCACACCATTATAGAcc                                                        >  4:183984/1‑139 (MQ=255)
                                                      aaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCACCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAATTCAAGTCTCGCAAACTACCGGCTCTAGCAAAACTCATGACTCTATGATAGACCg                                                       >  5:130828/1‑139 (MQ=255)
                                                                                               agcagcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaaaaaa                 >  6:294405/1‑139 (MQ=255)
                                                                                                   agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTTGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc               >  2:89381/1‑138 (MQ=255)
                                                                                                   agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGTAAGAGAATTCTTTGTCACAAAAAGAGTGAAAAAAcc               >  3:131201/1‑138 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGATTTATCGCAAACTTGCGGCTGGAGCAAATCTCATGGCTCTATGATAGCCCCGGGAAAGAGATTTCTTTGGGGGCATTTGAGTGaaaaaagcta             >  4:63901/1‑135 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                 >  7:29698/1‑55 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                 <  8:29698/55‑1 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACGCATGTATCGATGATAGTCCGTGGTAAGAGAATTCTATGTCACAATATGTGTGAAAAAACCaa             >  8:166851/1‑139 (MQ=255)
                                                                                                     cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAATAGAATTCTTTTTCACAATATGATTTAAAAAACCaa             >  1:147010/1‑138 (MQ=255)
                                                                                                               cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtttt     >  3:225632/1‑136 (MQ=255)
                                                                                                                aGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGTAAAGAGAATTCTTTGTCTAAATATGAGTGAAAAAACAAACTAGTTTTTAg  >  6:99217/1‑138 (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 < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: