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.T84.2% 95.0 / 14.8 42intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base T (17/15);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAA  >  CP006881/3725068‑3725321
                                                                                                                                             |                                                                                                                     
tGAGCCGAATAGACCGCCGCAGCCGCAGTTATTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACAACCCTCCGCGGCATTAAGGCAGGCAGcca                                                                                                                       >  3:99306/1‑138 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                     >  7:321536/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                     >  7:75082/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCtt                                                                                                                     >  5:310882/1‑139 (MQ=255)
       aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtttttag                                                                                                                >  3:150519/1‑135 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCGGCCCGtttt                                                                                                     >  1:159320/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGttt                                                                                                      >  5:268720/1‑138 (MQ=255)
                               ttCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcc                                                                                        <  1:34706/138‑1 (MQ=255)
                               ttCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcc                                                                                        <  8:321536/138‑1 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGCAAAGACGGCCGAtt                                                                   >  4:121211/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                   >  6:59322/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                   >  2:59642/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTTGCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAATACGGCCGAtt                                                                   >  8:139905/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGAgttagtt                                                              >  6:26501/1‑135 (MQ=255)
                                                           aaaccccacccacACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTGCTTTTTGGCAGGCAGCCTTTTTCGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAACGACGGCCGAtttcttttct                                                           <  1:41696/134‑1 (MQ=255)
                                                           aaaGACCAACCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCTTTTTCTCAGCCGCAGCCCGCAACAAACGACGGCCGAtttcttttct                                                           <  3:34100/139‑1 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                         >  3:85162/1‑139 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                         >  4:300285/1‑139 (MQ=255)
                                                                  aTCCACCCCCGTAAAGATACAAGCAACCGGTTCATAGAACAACATCCCTCCGCTGCTGTTTGGCGGGCAGCCGTTTTCGCAGCCCGTTTTGTCAGCCGCGGCCCGCAAGAGAAGGCGGCCGATGTCTTTGCTCGctatc                                                    >  4:43261/1‑135 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGGAGCCTTTTTAGCAGCCCGTTTTATCCGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTGTCGCAAACTAGCGGctcaa                                         >  5:352098/1‑137 (MQ=255)
                                                                                acaTACACGCAAACGTTTCATACAAAAACACCCCTCCGCTTGTTTTTGGCAGGCAGCCTTTTTACCACCCCGTTTTCACCCCCCCCCCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                      <  5:59322/137‑1 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCTAGCAAAAc                                  >  4:171456/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTGTAGCaaaa                                   >  6:336577/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                   >  7:347413/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                  >  3:273990/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAAATAGCGGCTCGAGCaaaa                                   >  5:214196/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTGCGTTTTGGCAGGCAGCCGGTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTGTTGTGTCGCCAACTAGCGGCTCGAGCaaacc                                  >  6:277896/1‑137 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTAGCAAAAcg                                 >  3:332048/1‑137 (MQ=255)
                                                                                        gaaaaCGTTTCAAACAAAAACATCCCCCCGCTTTTTTTTGGCAGGCAGCCTTTTTCGCCGCCCGTTTTCTCCGCCGCCGCCCCCAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  8:347413/137‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  7:11250/139‑1 (MQ=255)
                                                                                        gCAAACGTTTAATATAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCAGCCCGTTTTCTCACCCGCAGCCCCCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  7:197592/139‑1 (MQ=255)
                                                                                        ccAAAGGTTTCATACAAAAACATCCCCCCGCTTCTTTTTGGCCGGCCGCCTTTTTCGCCGCCGGTTTTCTCAGCCGCAGCCCCAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  3:303777/138‑1 (MQ=255)
                                                                                        ccAAACTTTTAATACAAAAACACCCCTCCGATTTTTTTTGGCCGGCACCATTTTTAGCCCCCCTTTTACCCCCCCGCCCCCCCAAAAAAAACACGGCCGATTTCTTTTCTCGCAAAAAAGCGGCTCTAGCAAAACTCAt                              <  5:90459/138‑1 (MQ=255)
                                                                                          aaagttttCATAAAAAAACACCCCCCCGCTTATTTTGGCCCGGCCGCTTTTTTCCCCCCCCGTTTGCTCCCCCCCCGCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                            <  5:336577/134‑1 (MQ=255)
                                                                                          aaaCGTTTCAAACAAAAACAACCCTCCGCTTCTTTTTGGCCGGCACCCTTTTTAGCAGCCCGTTTTCTCCGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                            <  7:7458/139‑1 (MQ=255)
                                                                                           aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                            >  4:349867/1‑138 (MQ=255)
                                                                                            ccGTTTCATCCAAAAACACCCCACGGCTTTTTTTTGCCCGCCCGCCTTTTTCCCCGCCCGTTTTCCCCGCCGCCCCCCGACAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  8:75082/138‑1 (MQ=255)
                                                                                            ccGTTTCATACAAAAACACCCCCCCGCTTCTTTTGGGCCGCCCGCCTTTTACCCCGCCCGTTTTCCCAGCCCCCCCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  4:150519/138‑1 (MQ=255)
                                                                                            ccGTTACATACAAAAAAAGCCCTCCGCTTATTTGTGCCCGCCCGCCTTTTTCGCCGCCCGTTTTCCCCGCCGCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  8:350348/138‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                           >  3:164273/1‑138 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCCCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  7:220345/139‑1 (MQ=255)
                                                                                            aCGTTACAAAAAAAAACAACCCCCCCCTTATTTTTGGCCGCCCCCCTTTTACCCCCCCCGTTTCCCCCCCCCCCCCCCCCACCAACACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  4:273990/139‑1 (MQ=255)
                                                                                             agTTTCAATACAAAAAAATCCCTCCTCTTTTTTTTGGCTGGCCGCTTTTTTCGCACCCGTTTTTCTCCGCCCCCGCCCCAAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  4:164273/138‑1 (MQ=255)
                                                                                                     tAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGGCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGCTCGAGCAAAACTCATGACTCTATGa                    >  1:278554/1‑137 (MQ=255)
                                                                                                             cATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg          <  6:120603/139‑1 (MQ=255)
                                                                                                             cATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg          <  3:349867/139‑1 (MQ=255)
                                                                                                                     cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTGTGGCCAAATAGCGGCTCGAGCAACACTCATGCCTAGATTGTATACCGGGGaaatataa  >  4:27638/1‑134 (MQ=255)
                                                                                                                                             |                                                                                                                     
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAA  >  CP006881/3725068‑3725321

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: