Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,798,160 T→C 14.9% M49T (ATG→ACG)  yicN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,798,1600TC14.9% 118.6 / 13.6 47M49T (ATG→ACG) yicNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/24);  new base C (3/4);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCAGCGATCGGCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAC  >  W3110S.gb/3798106‑3798227
                                                      |                                                                   
tcAGCGATCGGCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGt                                                     >  1:783065/1‑71 (MQ=255)
tcAGCGATCGGCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGt                                                     >  1:965910/1‑71 (MQ=255)
tcAGCGATCGGCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGt                                                     >  1:262441/1‑71 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:956304/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:1174954/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:795190/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:1211285/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:191117/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:1884629/71‑1 (MQ=255)
          gCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAAtttt                                           <  1:1813364/71‑1 (MQ=255)
                    gATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTAtt                                  <  1:954744/70‑1 (MQ=255)
                    gATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTAtt                                  <  1:895662/70‑1 (MQ=255)
                    gATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTAtt                                  <  1:439105/70‑1 (MQ=255)
                    gATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTAtt                                  <  1:1656168/70‑1 (MQ=255)
                          gTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCg                           <  1:1814094/71‑1 (MQ=255)
                          gTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCg                           <  1:1839002/71‑1 (MQ=255)
                          gTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCg                           <  1:1053358/71‑1 (MQ=255)
                          gTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCg                           <  1:1495975/71‑1 (MQ=255)
                             cccGTGGAGTCGATGGTCGATCAAACGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTcc                        >  1:1041905/1‑71 (MQ=255)
                             cccGTGGAGTCGATGGTCGATCAAACGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTcc                        >  1:335549/1‑71 (MQ=255)
                             cccGTGGAGTCGATGGTCGATCAAACGGGAAAGACAGCTgg                                                      >  1:1743415/1‑41 (MQ=255)
                              ccGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCg                       <  1:1735290/71‑1 (MQ=255)
                                   gAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGAt                    >  1:571456/1‑69 (MQ=255)
                                   gAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATg                   <  1:1748865/70‑1 (MQ=255)
                                   gAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATg                   <  1:151055/70‑1 (MQ=255)
                                     gTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGt                <  1:53179/71‑1 (MQ=255)
                                     gTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGt                <  1:1862896/71‑1 (MQ=255)
                                      tCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCg                       <  1:137647/63‑1 (MQ=255)
                                      tCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCg                       <  1:286657/63‑1 (MQ=255)
                                      tCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCg                       <  1:1632000/63‑1 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:1936224/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:340392/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:1335913/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:471183/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:1331018/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:1533654/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:1240855/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTg                                >  1:162718/1‑52 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGa                     >  1:1837782/1‑63 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGa                     >  1:1093026/1‑63 (MQ=255)
                                        gATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGa                     >  1:1062249/1‑63 (MQ=255)
                                            gTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCt         >  1:1195447/1‑71 (MQ=255)
                                                  cAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTAtt                                  <  1:1399176/40‑1 (MQ=255)
                                                   aaaCGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAc  <  1:969593/71‑1 (MQ=255)
                                                   aaaCGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAc  <  1:808564/71‑1 (MQ=255)
                                                   aaaCGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAc  <  1:526785/71‑1 (MQ=255)
                                                   aaaCGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAc  <  1:1837086/71‑1 (MQ=255)
                                                      |                                                                   
TCAGCGATCGGCTGAACATCGATGTCGTACCCGTGGAGTCGATGGTCGATCAAATGGGAAAGTCAGCTGGTGACGAATTTTTACGTTATTTGCATCGTCCGGATGAGTCGCACCTGCAAAAC  >  W3110S.gb/3798106‑3798227

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