Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,110,226 G→T 15.1% A98D (GCC→GAC)  mdoC ← membrane protein required for modification of periplasmic glucan

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,110,2260GT15.1% 93.8 / 4.1 33A98D (GCC→GAC) mdoCmembrane protein required for modification of periplasmic glucan
Reads supporting (aligned to +/- strand):  ref base G (20/8);  new base T (3/2);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACTTTCCGCTTTTCCTTTGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCAAGGGATAGCGTAAAA  >  W3110S.gb/1110160‑1110295
                                                                  |                                                                     
cAACTTTCCGCTTTTCCTTTGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGTCTGt                                                                   >  1:2213428/1‑71 (MQ=255)
 aaCTTTCCGCTTTTCCTTTGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGt                                                                   >  1:2292635/1‑70 (MQ=255)
      tCCGCTTTTCCTTTGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGTCTGATAACAt                                                             >  1:2234932/1‑71 (MQ=255)
                   tGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTaca                                                >  1:1770178/1‑71 (MQ=255)
                   tGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTa                                                  <  1:1646966/69‑1 (MQ=255)
                        tatTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGTCTGTTAACATCgg                                                          <  1:575587/56‑1 (MQ=255)
                        tatTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGTCTGTTAACATCgg                                                          <  1:2309352/56‑1 (MQ=255)
                        tatTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTc                                           <  1:1329214/71‑1 (MQ=255)
                        tatTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTc                                           >  1:1741327/1‑71 (MQ=255)
                          ttGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGtt                                                                  >  1:1726864/1‑46 (MQ=255)
                          ttGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGa                                         >  1:1207692/1‑71 (MQ=255)
                          ttGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGa                                         >  1:2573694/1‑71 (MQ=255)
                          ttGCAGCATAATAAAATGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGa                                         >  1:462056/1‑71 (MQ=255)
                            gcagcaTAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTaa                                                                >  1:941603/1‑46 (MQ=255)
                            gcagcaTAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTaa                                                                >  1:1421136/1‑46 (MQ=255)
                              agcaTAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAAc                                                               <  1:683154/45‑1 (MQ=255)
                              agcaTAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACg                                     <  1:1957946/71‑1 (MQ=255)
                                  taataaATTGCGGTAATGTCAGTAGGGGGATGTCTGTTAACATCGGGATACCTACACGTTCGACACg                                     >  1:1083440/1‑67 (MQ=255)
                                  taataaATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACt                                 >  1:2300848/1‑71 (MQ=255)
                                   aataaATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACtt                                >  1:2299411/1‑71 (MQ=255)
                                    ataaATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGa                                         >  1:150450/1‑61 (MQ=255)
                                    ataaATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACttt                               >  1:1346043/1‑71 (MQ=255)
                                          tGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTccacca                         >  1:306784/1‑71 (MQ=255)
                                           gCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTa                                   >  1:1221059/1‑60 (MQ=255)
                                           gCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCAt                        >  1:1409448/1‑71 (MQ=255)
                                             ggTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCAtt                       <  1:1383338/70‑1 (MQ=255)
                                                  tGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTa                                                  <  1:1741031/38‑1 (MQ=255)
                                                    tCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGact                                       >  1:1773740/1‑46 (MQ=255)
                                                    tCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCAAg               >  1:1473694/1‑71 (MQ=255)
                                                        tAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCAAg               <  1:1527606/67‑1 (MQ=255)
                                                         aGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACtt                                >  1:1018786/1‑49 (MQ=255)
                                                          gggggATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCaa                >  1:522831/1‑64 (MQ=255)
                                                              gATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTc                  <  1:143035/58‑1 (MQ=255)
                                                                  gCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCAAGGGATAGCGTaaaa  >  1:80615/1‑70 (MQ=255)
                                                                  |                                                                     
CAACTTTCCGCTTTTCCTTTGACATATTGCAGCATAATAAATTGCGGTAATGTCAGTAGGGGGATGGCTGTTAACATCGGGATACCTACACGTTCGACACGTACTTTCCACCATTTTTTCAAGGGATAGCGTAAAA  >  W3110S.gb/1110160‑1110295

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