Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,205:1 +T 100% intergenic (‑37/‑492) 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.8% 60.0 / 11.9 33intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base T (18/7);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 5.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT  >  CP006881/3725068‑3725337
                                                                                                                                           |                                                                                                                                       
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                        >  7:188762/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                        >  3:293952/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCtt                                                                                                                                      >  1:258828/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                                      >  4:258373/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGGAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGGGGGTTTTTGGCAGGCGGGCtt                                                                                                                                      >  4:294838/1‑139 (MQ=255)
       aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCtttttagc                                                                                                                                >  7:323965/1‑135 (MQ=255)
       aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtttttagc                                                                                                                                >  1:27181/1‑135 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACCTCCCTCCGCTGCTTTTTGGCAGGCCGCCTTTTTAGCAGccc                                                                                                                           >  7:207266/4‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATCCAAGCAAACGTTTCATAGAAAACCATCCCTCCGCTTGTTTTTGGCAGGCAGCCTTTTTagcagc                                                                                                                             >  6:127958/4‑137 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                           >  3:270498/4‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                                         >  3:286278/2‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACTTCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtc                                                                                                                     >  8:96541/1‑139 (MQ=255)
                   cAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTATTTAGCAGCCCGTTTtct                                                                                                                    <  2:4181/139‑1 (MQ=255)
                                             aCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACgg                                                                                         >  6:74602/1‑139 (MQ=255)
                                              cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                        >  7:157968/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGCCGTTTTAGCCGCCCGGTTTCTCAGCCGCCGCCCGCAAGCAAAGACGGCCGAtt                                                                                   >  6:72144/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTATCAGCCCGTTTTCTCAGCCCCACCCCGCAAGAAAAGACGGCCGAtt                                                                                   >  1:245536/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                   >  8:141976/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                              <  3:258373/139‑1 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGACCGTTTTCTCAGCCGCCGCCCCCAAGAAAAGACGGCCGAtttcgtt                                                                              >  5:244875/1‑139 (MQ=255)
                                                        caAAAAGAACATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTCGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGACCGCAAGAAAATACGACCGAtatattt                                                                              >  7:42767/1‑139 (MQ=255)
                                                                        acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCgg                                                              >  2:135/1‑139 (MQ=255)
                                                                          ccGAAAAGATACACGCAAACTTTTAATACAAAAACAGCCCTCCGTTTTTTTTTGGCCGGCCTCCTTTTTACCCGCCCGTTTTCCCCGCCGCAGCCCCCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGct                                                            <  6:244875/139‑1 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCATCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTTTAGCCTACTATCggggc                                                           >  3:8176/1‑134 (MQ=255)
                                                                                aGATAAAACCAAACTTTACATACAAAAACTTCCCTCGGTTTTTTTTTGGCTGGCTGCTTTTTTAGCAGCCCGTTTTCCCCGCCGCAGCCCGAAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                      <  4:293952/139‑1 (MQ=255)
                                                                                  aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                     >  3:136817/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGAGTTCGTTTGTCGCCAACTAGCGGCTCGAGGaaaa                                                   >  2:134569/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTATCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTATCGGCTCTAGCAAAAc                                                  >  1:42814/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                   >  2:176775/1‑138 (MQ=255)
                                                                                        gaaaaCTTTTAATAAAAAAACTCCCCACGGATTCTTTTTGGCCGGCCGCTTTTTTCCCCGCCGGTTTGCCCCGCCCCCGCCGCAAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  4:320349/137‑1 (MQ=255)
                                                                                        ccAAACGTGCCATCCAAAACCATCCCTCCGCTGCTTTTTGGCCGGCAGCCTTTTTAGCCGCCGGTGTTCTCCGCCGCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  8:188762/138‑1 (MQ=255)
                                                                                         cAAACGTTTCAAACAAAAACACCCCACCGCGTCTTTTTGGCAGGCACCCTTTTTACCCCCCCGGTTTCCCCCCCCCCCCCCCCAACACAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  7:79951/139‑1 (MQ=255)
                                                                                          aaccttttCATATAAAAACCTCGCGCCGTTTTTTTTGGGCCGGCCTCATTTTTCCCAGCCCTTTTTCCCCCCCGCCCCCCGAAAAAAACCACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                            <  8:136109/134‑1 (MQ=255)
                                                                                            ccGTTTCATACAAAAACATCCCACCGTTTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCCGTTTTCTCACCCGCCGCCCGCAATAAAAGACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  7:21023/138‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          >  1:279150/1‑139 (MQ=255)
                                                                                            aCGTTTCATACAAAAACACCCCTCCGCTTCTTTTGGGCTGGCCCCATTTTACGCCGCCGTTTTTCTCCGCCGCCGCCCCCAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  7:141976/139‑1 (MQ=255)
                                                                                                      gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                       <  2:223499/100‑1 (MQ=255)
                                                                                                      gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                       >  1:223499/1‑100 (MQ=255)
                                                                                                            caccccTCCGCTTCTTTTTGGCCGGCCGCCTTTTTAGCAGCCGGTTTTCTCCCCCCCAGCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg                          <  1:135/136‑1 (MQ=255)
                                                                                                                                    gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat  >  1:247157/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: