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% 158.4 / NA 48intergenic (‑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 (35/13);  total (35/13)
*CP0068813,725,2400GA93.6% 139.9 / ‑6.4 48intergenic (‑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 (35/9);  minor base T (0/3);  total (35/13)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAG  >  CP006881/3725119‑3725366
                                                                                                                             ||                                                                                                                              
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCATCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                >  8:239090/1‑139 (MQ=255)
   atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                >  5:214303/1‑136 (MQ=255)
     caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttcttt                                                                                                           >  7:265628/1‑139 (MQ=255)
     caAAAAGACCATCCACACCCGAAAAGATACAACCAAACGTTGCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCCGCCCGTTTTCTCCGCCCCACCCCGCAAGAAAACACGGCCGAtttc                                                                                                              <  4:185941/136‑1 (MQ=255)
        aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttgttttct                                                                                                        >  1:129917/1‑139 (MQ=255)
        aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                           >  5:145792/1‑136 (MQ=255)
             ccATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGGTTTCTCAGCCGCATCCCGCAAGAAAAGACGTCCGAtttgttttggagca                                                                                                    >  8:384754/1‑132 (MQ=255)
                tCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                    >  5:308901/1‑135 (MQ=255)
                      cccGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                          >  8:317759/1‑139 (MQ=255)
                        accacggATATATGCAAAGGTTTATTATAAAAAATACCCTGCGTTTTTTTTTGGCCGGCAGGATTTTTAGCTGCCCGTTTTCTCTCCCGCTCCGCGCATCAAACCTCGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                        <  8:326442/133‑1 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAACCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGTCCGATGTCGTTTCTCGCAAACTAGCCGcgc                                                                                        >  7:27174/1‑135 (MQ=255)
                             aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                        >  3:178612/1‑134 (MQ=255)
                               aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAg                                                                                     >  5:413402/1‑135 (MQ=255)
                                 aCAACCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGTCAGGCAGCCTTTTTTGCAGCCCGTTTTCACAGCCGCAGCCCCCAAGAAAAGACGGCCGAATTATTTTCTCGGAAACTAGCGGCTCATGCaaaa                                                                                >  7:9713/1‑138 (MQ=255)
                                  cAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCat                                                                                  >  4:116548/1‑134 (MQ=255)
                                    aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCTATAGCAAAAc                                                                               >  4:211782/1‑136 (MQ=255)
                                     gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                           <  7:317759/139‑1 (MQ=255)
                                     ccAAACGTTTCAAACAAAAACATCCCCCCGCTTCTTTTGGGCAGGCCGCCTTTTTCGCCGCCCGTTTTCTCCCCCCCCCCCCGAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                           <  2:283620/138‑1 (MQ=255)
                                       acaCGTTGAATACAAAAACAGCCCCCCGCTTTTTTTGGGCCGGCAGCTTTTTTCGCCGCCGGTTTTCGCCGCCGCCGCCCCCACCAAACCACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCa                                                                            <  6:308901/134‑1 (MQ=255)
                                       aaaCTTTTCATAAAAAAACAGCCCTCCGCTTCTTTTTGGCAGGCCGCGTTTTTAGCCGCCCGTTTACTCAGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                            <  3:211782/136‑1 (MQ=255)
                                       aaaCGTTTCATACAAAACCATCCATCCGGTGATTTTTGGCCCGCAGCCTTTTTCGCAGCCGTTTTTCTCAGCCGCCCCCCGCAACAAACCACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                            <  6:413402/136‑1 (MQ=255)
                                         ccttttAATATAAAAACATCCCTCCGTTTGTTTTTGGCTGGCAGCTTTTTAAGCTCCCGTTTTTCTCAGCCGCCGCCCCAATTAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                          <  6:214303/133‑1 (MQ=255)
                                         ccGTTTAATAGAACAACACCCCTCCGGTTCTTTTGGCCAGGCAGCCTTTTACCCTCCCCTTTTTCTCCGCCGCCGCCCCAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAACCTCATGAct                                                                       <  8:119151/138‑1 (MQ=255)
                                         aggTTTAATACAAAAACATCCCTCGGCTTATTTTGGGCAGGCAGCTTTTTTCGCTGCCCGTTTTCTCAGCCGCAGCCCGCATGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       <  7:239090/137‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTCGCAAACAAGCGGCTCTAGCAAAACTccagtc                                                                        >  7:326442/1‑133 (MQ=255)
                                         aCGTTTAATAGAAAAACACCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTCGCAGCCCGTTTTCTCCCCCGCAGCCCCCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       <  7:229833/139‑1 (MQ=255)
                                         aCGTTACATACAAAACCCTCCCTCCGCTTCTTTTTGGCAGCCCGCATTTTTCGCCCCCCTTTTCCTCCGCCGCCCCCCGAAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       <  2:139459/139‑1 (MQ=255)
                                            tttCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       >  3:216040/1‑136 (MQ=255)
                                            tttCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTATCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       >  4:169773/1‑136 (MQ=255)
                                            tttCATAGAAAAACATCCCTCCGCTTCTTTTTGGCACGCCGCCTTTTTATCAGCCCGTTTTCTCAGCCGCAGCCCTCAACAAAAGACGGAGGAATTCTTTTATCGAAATCTGTCCGCTCGAGCAAAAATCATGCct                                                                       >  6:51714/1‑136 (MQ=255)
                                              tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCGAGCAAAACTGATGACGCTATg                                                                  >  2:54646/1‑139 (MQ=255)
                                                 tAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAcgagat                                                                   >  6:315501/1‑130 (MQ=255)
                                                    aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAGTTCGTTTGTCGCACACTAGCGGCTGTAGCAAAACTCATGGCTGTTTGTTAGAc                                                            >  2:301827/1‑139 (MQ=255)
                                                       aaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                                                            >  6:303633/1‑136 (MQ=255)
                                                            cccTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTCGCCAACTAGCGGGACTAGGAAAACTCATGACTCAATGATAGCCCGTGGaaa                                                    >  8:223844/1‑139 (MQ=255)
                                                                      cttTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATTGACCGTGGAAAGAGaaa                                              >  3:388116/1‑134 (MQ=255)
                                                                         tttGGCAGGCAGCCTTTTTAGAAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGGTCTAGCAAAACTCAAGACTCTATGATATAACCTGGAAAGAGAATTCCTTGt                                       >  8:192367/1‑139 (MQ=255)
                                                                          ttGGCAGGCTGCCTTTTTAGCACCCCGTTTTCTCTGCCGCAGCCCCAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCtt                                          <  4:388116/135‑1 (MQ=255)
                                                                                                   agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGCCTCTATGATAGTCCGTAGAAAGAGCATTGGTTGTGCCAATATGAGTGAAAAAAcc                  >  1:173660/1‑138 (MQ=255)
                                                                                                   agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCGCGCCAACTAGCGGCTCTAGCAAACCTCATGCCTCTATGATAGACCGTGGAAAGGGAATTAGTTGTCACAATATGAGTTaaaaac                    >  4:313/1‑135 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa                >  7:314573/1‑139 (MQ=255)
                                                                                                    gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCa                 >  1:141905/1‑138 (MQ=255)
                                                                                                      ccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTGTCGCAAACTAGCGGCTCTAGCAAAACTCCTGACTCTATGATAGACCGGGGAAAGAGAATTCTTTGGGACCATATTAGGGAAAAAAcc                  >  4:5245/1‑135 (MQ=255)
                                                                                                        gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCTTGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc               >  5:175426/1‑136 (MQ=255)
                                                                                                        gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGCCTCGATGATAGACCGGGGAAAGAGAATTCTTTGTCAGAATATGAGTGATAAAACCGAc               >  3:280245/1‑136 (MQ=255)
                                                                                                        gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc               >  3:425342/1‑136 (MQ=255)
                                                                                                                     cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTTTGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGGTTTTAGAAg  >  4:370704/1‑136 (MQ=255)
                                                                                                                     cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCCCAAAAGTGCGGCTCTAGCAAAACTGATGACTCGATGATAGTCCCTGGAAAGAGAATTATTTGTGAACAAATTACTGAAAAAAACCACGACTGTTtagtag  >  4:63259/1‑136 (MQ=255)
                                                                                                                             ||                                                                                                                              
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAG  >  CP006881/3725119‑3725366

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: