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.T88.1% 124.1 / 10.2 47intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (29/8);  total (36/11)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATG  >  CP006881/3725068‑3725338
                                                                                                                                            |                                                                                                                                        
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                         >  8:424242/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                         >  2:105878/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                         >  8:201361/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                         >  8:413058/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGgccgccgt                                                                                                                                         >  5:125388/1‑137 (MQ=255)
         tAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTagcag                                                                                                                              <  7:413058/139‑1 (MQ=255)
            accgCGGCAGCCGCAGTTTGTCTGGCCGGCGTCACCCAGATCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                           <  7:23586/136‑1 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCCCCCAGGTCATCAAAAAGACCCTCCCCACCCGTAAAGATACAAGCAAACGTGTCATAGAAAAACCTCCCTCCGCTTCTTTGTGGCAGGCGGCCTTTTTAGCAGCCCGt                                                                                                                         >  5:1058/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCga                                                                                                                         >  3:345704/2‑138 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCAGTTTAGGGGCCCGgttt                                                                                                                      >  1:444573/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTGTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCACCCTTTTTATCAGCCCGtttt                                                                                                                      >  1:89448/1‑139 (MQ=255)
                                        gCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGGCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGaaaa                                                                                               >  5:352306/1‑139 (MQ=255)
                                              cccAGTTCACCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGACGCAGCCCGCAAGAAAAGACgcc                                                                                         >  3:123241/1‑137 (MQ=255)
                                                   ttCATCTAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGGCTGTTTTGCAGGCAGCCGTTTTATCAGCCCCTTTTCTCAGCCGCAGCCCGCAAGCACAGACGGCAGagt                                                                                    >  5:490599/1‑137 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCg                                                                          >  4:26887/1‑139 (MQ=255)
                                                               aCCATCCACACCCTAAAAGATACAACCAAACGTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCa                                                                        <  7:145964/139‑1 (MQ=255)
                                                                  aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGCAAAGACGGCCGAGTTCTTTTGTCGCAAAc                                                                     >  2:192137/1‑139 (MQ=255)
                                                                   tCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc                                                                     >  3:349710/1‑138 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCg                                                                                                              <  2:438637/87‑1 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCg                                                                                                              >  1:438637/1‑87 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                            >  3:200315/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCGGGCGGCCTTGTTTGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGCAAAGCCGGCCGAGGTGTTTTGTCGCCTACTAggggggc                                                            >  6:509578/1‑134 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTATCGCAAACTAGCGGctctt                                                          >  8:303161/1‑138 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGGTTTAGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGCTTTCTTTTCTCGCAAACTAGCGGctc                                                            >  7:499840/1‑137 (MQ=255)
                                                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                       <  7:424242/139‑1 (MQ=255)
                                                                                 gATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGCa                                                       >  1:206543/1‑138 (MQ=255)
                                                                                  aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                      >  3:201986/1‑138 (MQ=255)
                                                                                  aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTGCTTTTTGGCCGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTGCTTTTCTCGGAAACTAGCGGCGCGAGcga                                                      >  4:212046/1‑136 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGATAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGAGGGCCGATTTCTTGTGTCGCAAGCGAGCGGCTCTAGCAAAAc                                                   >  8:387818/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCGGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGGTCTGTTGTCGCAAACAAGCGGCGCGAGcccaa                                                    >  2:184380/1‑134 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTCTCGCAAACTAGCGGCTCGAGCaaaa                                                    >  3:422769/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                    >  6:390499/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTGTAGCAAACAAGCGGCGCTAGCAAAACt                                                  >  5:421838/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACt                                                  >  8:242052/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGCAAAACt                                                  >  7:128358/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAATTGTTTTCACACAAAATAGCGGCGCAAGCAAAACt                                                  >  7:196187/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTTGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCAGCCCCCAAGAAAAGACGGCCGATTTAGTTTCTCGCACAAGAGCGGCACTAGCAAAACt                                                  >  6:4979/1‑138 (MQ=255)
                                                                                          aaaCGTTAAATCCAAAAACACCCCCCCGTTTATTTTTGCCAGGCAGCCTTTTTACACCCCCGTTTTCTCCGCCCCCGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                             <  1:105878/139‑1 (MQ=255)
                                                                                            agttttCATAGAAAAAATTCCCTCGGCTTATTTTTGGCAGGCAGCATTTTACGCAGCCCGTTTTCTCCGCGGCAGCCCGAAATAAATCACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                           <  8:128358/136‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAATACGGCCGATTTCTTTTATCGCAAAATAGCGGCTCTAGCAAAACTCATGAct                                           >  4:76016/1‑139 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTATCCCAAACTAGCGGCTCTAGCAAAACTCCAGAct                                           >  4:234163/1‑139 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCAGCAGCACGCAAGAAAAGACGGCAGATTACTTTTCTCGCAAACTATCGGCTCAAGCAAAACACATGAc                                            >  3:501468/1‑138 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCCCAAACTAGCGGCTCTAGCAAAACAAATGAc                                            >  3:64783/1‑138 (MQ=255)
                                                                                             cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                           >  6:243099/1‑138 (MQ=255)
                                                                                                       aaaaaCATCCCTCCGTTTTTTTTTGGCTGGCCGCTTTTTTCGCCGCCGGTTTTCTCCGCCGCCGCCCCAAATAAAACACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                                <  4:201986/139‑1 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                        >  4:79697/1‑137 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                        >  3:424544/1‑137 (MQ=255)
                                                                                                                    cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa                   >  6:199373/1‑139 (MQ=255)
                                                                                                                    cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATAGCAAAAAATCGGGTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa                   >  4:312448/1‑139 (MQ=255)
                                                                                                                         cttTTTGGCCGCCACCCTTTTTACCCGCCCTTTTTCTCCGCCGCCGCCCCAAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt               <  4:185280/138‑1 (MQ=255)
                                                                                                                               ggcaggcaGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTca         <  5:199373/138‑1 (MQ=255)
                                                                                                                                     cagCCTTTTTAGCAGCCCGTTTACACAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATg  <  1:11347/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATG  >  CP006881/3725068‑3725338

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: