Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,104,736 C→T 47.4% intergenic (‑339/‑55) opgC ← / → opgG osmoregulated periplasmic glucan succinylation membrane protein/osmoregulated periplasmic glucan (OPG) biosynthesis periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,104,7360CT47.4% ‑4.8 / 18.4 19intergenic (‑339/‑55)opgC/opgGosmoregulated periplasmic glucan succinylation membrane protein/osmoregulated periplasmic glucan (OPG) biosynthesis periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base T (5/4);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GAAGTCGATGCTTCTGTCTTTAGGAGAAGCACGGAAAGTGAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTGTATACATCTTCAAGCTGGGCTTTCAGTATTGATGATGTCGCA  >  CP009273/1104595‑1104877
                                                                                                                                             |                                                                                                                                             
gAAGTCGATGCTTCTGTCTTTAGGAGAAGCACGGAAAGTGAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATaa                                                                                                                                        <  1:191875/149‑1 (MQ=255)
    tCGATGCTTCTGTCTTTAGGAGAAGCACGGAAAGTGAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGcac                                                                                                                                    <  1:341413/149‑1 (MQ=255)
                      ggAGAAGCACGGAAAGTGAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCtt                                                                                                                  >  2:127590/1‑149 (MQ=255)
                                 gAAAGTGAAGACGGTGGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACCCAAAGGGGGAAGTGCTTACTAATTATGa                                                                                                       <  1:345619/149‑1 (MQ=255)
                                       gAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATa                                                                                                 <  2:259747/149‑1 (MQ=255)
                                                 gCAATCAGGTGCTTTATCCATGAGCCAGCGTGCTGAACGATACAGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGGTTACTAATTATGAAACATAAACTACAAat                                                                                       <  2:224588/149‑1 (MQ=255)
                                                   aaTCAGGTGCTTAATCCATGAGCCTGCGTGACGACAGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTATAATTGATCGTATCGATATAAGCACACAAAGGGGGAAGGGCTTACTAATTATGAAACATAAACTACAAatga                                                                                     <  2:120668/149‑1 (MQ=255)
                                                         gTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGaaaa                                                                               <  2:288301/149‑1 (MQ=255)
                                                            cTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATTGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAAGATGAGGTACAGATGATGaaaccgc                                                                            >  2:280401/1‑145 (MQ=255)
                                                                          cAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGtgct                                                              <  1:330913/149‑1 (MQ=255)
                                                                             cGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCa                                                           <  1:127590/149‑1 (MQ=255)
                                                                                  tGAACGATACCGGGATTCTGTTGGCGGAATGGCTGGTTATCCATTAAAGTAGATCGGGACGATATAAGCACACAAAGAGGAAAGTGCATACTAATTACAGAAGATAAACTACACACGATAACGAAGCGTTGGTGGATtggtgcgtgaag                                                      >  1:102447/1‑143 (MQ=255)
                                                                                           ccGGGATTCGTGTGTCGCAATGGTTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACGACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCt                                             <  2:52594/149‑1 (MQ=255)
                                                                                           ccGGGAGTCTGTTGTCGGAATGGCTGGTTATCCATTAAAAGAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCt                                             <  2:152042/149‑1 (MQ=255)
                                                                                            cGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTg                                            >  2:30333/1‑149 (MQ=255)
                                                                                                     gttgtCGGAATGGCTGGTTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCTTACTACTTATGAAACATAAACTACAAATGAGGAAAATGCGTTGGTTGCGTGCTGCAGTAATGTTAACCCTGTATACATCt                                   >  2:138062/1‑149 (MQ=255)
                                                                                                        gtCGGAATGGCTGATTATCCATTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAGACTACAAATGATGAAAATGCGTTGGTTGAGTTCTGCAGTAATGTTAACCCTGTATACATCTTCa                                >  2:244165/1‑149 (MQ=255)
                                                                                                                            aTTAAAATAGATCGGATTGATATAAGCACACAAAGGGGGAAGTGATTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTATATACATCTTCAAGCTAGGCTTTCCGTATtga            >  1:108435/1‑149 (MQ=255)
                                                                                                                                   tAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTGTATACATCTTCAAGCTGGGCTTTCAGTATTGATGATGTc     <  2:148681/149‑1 (MQ=255)
                                                                                                                                      atcggatcgATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTGTATACATCTTCAAGCTGGGCTTTCAGTATTGATGATGTCGCa  >  1:300870/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GAAGTCGATGCTTCTGTCTTTAGGAGAAGCACGGAAAGTGAAAACGGTTGCAATCAGGTGCTTAATCCATGAGCCAGCGTGCTGAACGATACCGGGATTCTGTTGTCGGAATGGCTGGTTATCCATTAAAATAGATCGGATCGATATAAGCACACAAAGGGGGAAGTGCTTACTAATTATGAAACATAAACTACAAATGATGAAAATGCGTTGGTTGAGTGCTGCAGTAATGTTAACCCTGTATACATCTTCAAGCTGGGCTTTCAGTATTGATGATGTCGCA  >  CP009273/1104595‑1104877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: