Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,209 +T intergenic (‑41/‑488) 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,2051.T83.0% 124.3 / 21.7 51intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base T (28/11);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT  >  CP006881/3725068‑3725337
                                                                                                                                             |                                                                                                                                      
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCAGGCCGCCt                                                                                                                                       >  1:183024/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCGGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCtt                                                                                                                                      >  8:282206/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTTGCAGGCAGCCtt                                                                                                                                      >  1:366287/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACACAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACTTTTCATATAAACAAAGCCCTCCGCGTCTTTTTGTCAGGCAGCCtt                                                                                                                                      >  4:359149/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCAGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGAGACAAGCAAACTTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAGCCta                                                                                                                                      >  7:430080/1‑138 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                          >  8:46607/4‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTGTTAGCAGCCCGtgt                                                                                                                      >  3:8433/1‑136 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                                     >  8:275588/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                                     >  3:218627/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                                     >  3:172803/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCAACAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAGACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGtttt                                                                                                                     >  2:156575/1‑139 (MQ=255)
                           gTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCTTTTTATCAGCCCGTTTTCTCAGCCGCa                                                                                                           >  8:307621/1‑139 (MQ=255)
                               ttCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcc                                                                                                        <  7:370443/138‑1 (MQ=255)
                                 cTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCTCGTTTTCTCAGCCGCAGcccg                                                                                                      <  7:65576/138‑2 (MQ=255)
                                 cTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                                      <  2:43838/138‑2 (MQ=255)
                                    tCCGGCGTCCCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAACCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCATTTTAGCAGCCCGTTTTCTCAGCCGCAGccagcaa                                                                                                   >  5:231880/1‑133 (MQ=255)
                                      cGGCGTCCCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTTGCCGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccgcaaga                                                                                                 >  5:100043/1‑132 (MQ=255)
                                      cGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGcccgcaaga                                                                                                 >  8:39020/1‑132 (MQ=255)
                                              cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGGAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                        >  5:192803/1‑139 (MQ=255)
                                              cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                        >  5:199749/1‑139 (MQ=255)
                                                  gTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCATTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAt                                                                                    >  1:299404/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                   >  8:19855/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGCCGGCCGAgttctgt                                                                              >  4:26443/1‑137 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttcttt                                                                              >  1:405840/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttattt                                                                              >  1:383664/1‑139 (MQ=255)
                                                                  aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc                                                                    >  4:422918/1‑139 (MQ=255)
                                                                         cccGTCAAGAAACACCCAACCTTTTCATACAAAAACATCCCTCGGCTTGTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTACTCAGCCGCAGCCCTCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                             <  7:39020/139‑1 (MQ=255)
                                                                         cccGAAAACAAAAAACCAAACGTTCCAAACAAAAACATCCCGCGGCTTCTTTTTGGCCGCCAGCATTTTACGCCCCCCTTTTTCTCAGCCGCCCCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                             <  7:402937/139‑1 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACAAGCGGctc                                                           >  2:355377/1‑137 (MQ=255)
                                                                                aGATACATCAAAACGTTTAAAAAAAAAACACCCCTCCGCTTCTTTTTGGCCGGCCGCCTTTTTCGCAGCCCTTTTTCTCCCCCCCCGCCCGCACAAAATCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                      <  4:218627/139‑1 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTTGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTGTTTTGTCGCCAACTCGCGGCTCTATCaaa                                                    >  6:33227/1‑137 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTCTTTTCTCGCAACCGAGCGGCTCTAGCAAAAc                                                  >  5:373337/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                  >  8:103549/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTATCATAGAAAAACATCCCTCCGCTTGTTTTAGGCAGGCCGCCTTTTTTGCAGCCCGTTTTCTCCCCCGCAGCCCGCAAGAAAAGACAGCCCAATAGGTGTCTCGCAAACTAGCGGCTGGACCaaaa                                                   >  4:189834/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTTGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCCAGTTCGTTTCTCGCAAACAAGCGGCTCGACCaaacag                                                 >  1:85066/1‑135 (MQ=255)
                                                                                        gaaaaCGTTTCATATAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCATTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCCCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  7:307621/137‑1 (MQ=255)
                                                                                        gCAACCGTTTCATAGAAAAACCTCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCACCCCGTTTGCTCAGCCGCCGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  8:386712/139‑1 (MQ=255)
                                                                                        gCAAACTTTTAATACAAAAACACCCCTCCGCTTTTTTTGGGCCGGCAGCCTTTTACCCCCCCCGTTTTCTCCCCCCCCGCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  1:355377/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATATAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCGGTTTTCTCACCCGCAGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  3:422918/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  7:275588/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTATTTTCTCGCAAACTAGCGGCTCTAGCCAAACTCAt                                              >  6:68866/1‑139 (MQ=255)
                                                                                        ccAAACGTTTCATACAAAAACATCCCGCCGCTTCTTTTTGGCCGCCAGCATTTTTCCCACCCCGTTTTCTCAGCCCCCCCCCCAAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  7:46607/138‑1 (MQ=255)
                                                                                        ccAAACGTGTCAAACAAAACCATCCCTCGGGTACTTTTGGCGCGGCAGCCTTTTTAGCACCCGGTTTTCTCCGCCGCTGCCCGCCAGAAAAGACGCCCGATTTATTTTCTCGCAAACAATCGGCTCTAGCAAAACTCAt                                              <  6:159651/138‑1 (MQ=255)
                                                                                        caaaaCGTTAAATATAAAAACATCCCTCCGCTTCTTTTTGGCAGCCAGCCTTTTTCCCAGCCCGTTTTCACCGCCGCAGCCCGAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  2:299404/137‑1 (MQ=255)
                                                                                         ccaaCGTTTACAAACAAAAACACCCCTCCGATTTTTTTTGGCGCGCCCCCTTTTTCCCCGCCTGTTTACTCCGCCCCACCCCCCAACAAAACACGGCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAAACTCAt                                              <  6:429304/137‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAGAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAcg                                          >  2:381336/1‑138 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          >  5:429304/1‑139 (MQ=255)
                                                                                            aCGTTTAATACAAAAACACCCCCCCGCTTCTTTTTGCCTGCCAGCTTTTTTACCACCCGGTTTTCACCGCCCCAGCCCCAAACAAACGCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  7:282206/139‑1 (MQ=255)
                                                                                             cGTTTCATAGAAAAACATCCATCCGCTTATTTTTGGCAGGCAGACTTTTTAGCAGACCGTTTTCTCAGCCGCAGCCCGCAAGAAAATACGGCCGATTTGTTTTCTCCCAAAATAGAGGGTCGAGCAAAACTCATTActc                                         >  1:139961/1‑139 (MQ=255)
                                                                                                 tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATTACTCTAt                                      >  2:286095/1‑138 (MQ=255)
                                                                                                   cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTAATGACTCTATGa                                    >  3:98079/1‑139 (MQ=255)
                                                                                                          aaaCATCCCTCCGCTTCTTTTTGGCAGGCAGACTTTTTAGCAGCCCGATTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTAGCAAACTAGCGGCTCTAGCAAAACTCATGACTCAAAGATAGTCCg                             >  8:446484/1‑139 (MQ=255)
                                                                                                                            ttttGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTATCGCAAACTAGCGGCGCGCGCAAAACTCATGACTCGATGATAGACCGTGGAAAGAGAATTGtttt           >  2:229296/1‑138 (MQ=255)
                                                                                                                             tttGGCAGGCAGCCTTTTTAGCCGCCCGTTTTATCAGCCGCAGCCCGCCAGAAAAGACGGCCGATTTAGTTTATCGCCAACAAACGGCGCTAGCAAAACACAATCCTCTATTAAAGACCCTGGAAAGACAAttaatttt          >  4:421995/1‑133 (MQ=255)
                                                                                                                                     gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat  >  5:170270/1‑139 (MQ=255)
                                                                                                                                     gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat  >  3:107167/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                      
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT  >  CP006881/3725068‑3725337

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: