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,2390GA93.5% 77.8 / ‑5.7 31intergenic (‑71/‑458)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/6);  minor base C (2/0);  total (25/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*CP0068813,725,2400GA96.8% 87.8 / NA 31intergenic (‑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 (25/5);  minor base C (0/1);  total (25/6)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAG  >  CP006881/3725119‑3725363
                                                                                                                           ||                                                                                                                           
ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGAtt                                                                                                             >  2:239757/1‑139 (MQ=255)
     caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttattt                                                                                                        >  4:189652/1‑139 (MQ=255)
            acactCCACACCCGAAAAGAAACACCCAACCGTTTAAAACAAAACCACCCCACCGCTTTTTTTTGGCCGGCCGCCTTTTTCGCAGCCGGTTTGCTCCGCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCa                                                                                                 <  8:264524/135‑1 (MQ=255)
             ccATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTATTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                 >  6:323195/1‑138 (MQ=255)
                     acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGAGTTCTTTTCTGGCAAACGAGCgg                                                                                        >  4:236367/1‑139 (MQ=255)
                     acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTCGCAAACTAGCgg                                                                                        >  5:109392/1‑139 (MQ=255)
                         gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGAGTTGTTTTGGGGCCCACTAGCGGcgcg                                                                                    >  1:162328/1‑136 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGctc                                                                                     >  4:159767/1‑137 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTGTCTCGCAAACTAGCGGctc                                                                                     >  7:130844/1‑137 (MQ=255)
                          tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCCTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGAGGTCTTTTCGACCCAACTATCGGcgg                                                                                     >  1:178624/1‑135 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                            >  1:51115/1‑139 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                             >  7:49466/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCAGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTATAGCaaaa                                                                             >  8:249861/1‑138 (MQ=255)
                                 aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTGTTTGGCAGGCAGCCGTTTTATCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTGAGTTTAGACCAAACTAGCGGCACTAGCAAAAc                                                                            >  5:130502/1‑139 (MQ=255)
                                     ccAACCGTTACATAAAAAAACACCCCCCCGCTTTTTTTGGCCCGGCAGCTTTTTTCGCCCCCCTTTTTCCCCCCCGCCGCCCGCAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                        <  1:226258/138‑1 (MQ=255)
                                       aacgttttCATAGAACAACCACCCGCCGCTTCTTTTTGGCGGCCCGCCTTTTTAGCCGCCCGTTTGCTCAGCCCCCCCCCCCACCAAACGCGGGGCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                      <  6:130502/134‑1 (MQ=255)
                                         ccGTTTCAAAGAAAAACACCCGGCGCCTGCTTTTGGGCAGCCCCCTTTTTAAGCAGCCGTTTTGCACCGCCGCTCCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    <  2:37583/138‑1 (MQ=255)
                                         aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGAct                                                                    >  6:324444/1‑139 (MQ=255)
                                          cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAACACTCATGActc                                                                   >  4:155112/1‑139 (MQ=255)
                                          cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActc                                                                   >  2:226258/1‑139 (MQ=255)
                                                atatAAAAACATCCCTCCGCTTCTTTTTGGCCGGCCGCCTTTTTACCCGCCCTTTTTCTCCGCCGCCCCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                                             <  1:37736/139‑1 (MQ=255)
                                                      aaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAGTTGGTTTATCGCCAACTAGCGGCGGGAGCAAAACTCCTGACTCTAAGTTAGGCCg                                                       >  1:133170/1‑139 (MQ=255)
                                                                      cttTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCGTTTCTCGCAAACTAGCGGCGCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCtt                                       >  1:225884/1‑139 (MQ=255)
                                                                      cttTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGACGCGATGTTAGACCGTGGAAAGAGAATTCtt                                       >  1:102616/1‑139 (MQ=255)
                                                                                ggcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGCCCGTGGAACGAGCATTCTTTGTGACAata                             >  1:143703/1‑139 (MQ=255)
                                                                                 gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat                            >  8:240950/1‑139 (MQ=255)
                                                                                  cagCCTTTTTCGCAGCCCGTTTTCCCAGCCGCACCCCTCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATg                           <  3:159767/139‑1 (MQ=255)
                                                                                                  gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa             >  8:291403/1‑139 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTTTCACAATATGATTGAAAAAACCaa             >  7:69634/1‑138 (MQ=255)
                                                                                                   cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTAACAATATGAGTGAAAAAACCaa             >  1:101619/1‑138 (MQ=255)
                                                                                                             cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGCCCGTGGAAAGAGAATTCTTTGTCCCAATATGTGTGAAAAAACCAACTAGgtttttg  >  1:280239/1‑137 (MQ=255)
                                                                                                                           ||                                                                                                                           
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAG  >  CP006881/3725119‑3725363

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: