Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 515,934 C→T 49.4% A43V (GCG→GTG)  ybbM → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb515,9340CT49.4% 0.9 / 126.9 87A43V (GCG→GTG) ybbMpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (16/28);  new base T (22/21);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAG  >  W3110S.gb/515881‑515997
                                                     |                                                               
ccATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatta                                                <  1:2532352/71‑1 (MQ=255)
ccATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatta                                                <  1:147247/71‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:2426561/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:2430486/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:2534481/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:1847815/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:1686344/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:339819/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCt                                                     <  1:549246/63‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGat                                                  <  1:1063470/66‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGat                                                  <  1:2356829/66‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGat                                                  <  1:3245296/66‑1 (MQ=255)
   tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGat                                                  <  1:1393253/66‑1 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAATGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:3005508/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:2606576/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:2129658/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:3212561/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:1277896/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:1437993/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:250967/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:2602386/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:3070702/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:324925/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:2962809/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                           >  1:2543102/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                           >  1:603037/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                           >  1:2864569/1‑71 (MQ=255)
     aaGAAAAACTGGCGCTGGAGAAAGATATTCGCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                           >  1:2687163/1‑71 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:1910415/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:3199827/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:693927/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:6275/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:1819960/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:1797193/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:2841717/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:1562064/1‑70 (MQ=255)
         aaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGc                                        >  1:56557/1‑70 (MQ=255)
          aaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAAt                                                           >  1:2627854/1‑50 (MQ=255)
             cTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtgt                                   <  1:2354852/71‑1 (MQ=255)
             cTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtgt                                   <  1:293696/71‑1 (MQ=255)
             cTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtgt                                   <  1:2350827/71‑1 (MQ=255)
             cTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtgt                                   <  1:1653742/71‑1 (MQ=255)
             cTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCAGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtgt                                   <  1:1914565/71‑1 (MQ=255)
                 cgcTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTg                               <  1:796508/71‑1 (MQ=255)
                  gcTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTg                               <  1:345000/70‑1 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCg                                          >  1:2040903/1‑50 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:3166299/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:346584/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:2691079/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:611396/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:739277/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:1068754/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:2006955/1‑71 (MQ=255)
                           aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                     >  1:149805/1‑71 (MQ=255)
                           aGATATTATCTGGAGCGTCGGGCGAGTGATAATTCAGCTGa                                                   >  1:2297059/1‑41 (MQ=38)
                                     tGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg             <  1:2417738/69‑1 (MQ=255)
                                     tGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg             <  1:3295412/69‑1 (MQ=255)
                                     tGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg             <  1:420609/69‑1 (MQ=255)
                                     tGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg             <  1:1382516/69‑1 (MQ=255)
                                      ggAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg             <  1:2882460/68‑1 (MQ=255)
                                       gAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGTCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1646993/70‑1 (MQ=255)
                                       gAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1290657/70‑1 (MQ=255)
                                       gAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:631978/70‑1 (MQ=255)
                                       gAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1075769/70‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:992220/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1427499/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1263200/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1121979/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:96159/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1776016/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2627745/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2574628/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2271139/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:30546/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2738646/67‑1 (MQ=255)
                                          cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2565075/67‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:18212/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:265387/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:782149/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2403325/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:2171543/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1579686/65‑1 (MQ=255)
                                            tCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat          <  1:1620668/65‑1 (MQ=255)
                                              gggCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgatga        >  1:2564530/1‑65 (MQ=255)
                                              gggCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgatga        >  1:142453/1‑65 (MQ=255)
                                              gggCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgatga        >  1:897628/1‑65 (MQ=255)
                                               ggCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAg  <  1:2428693/70‑1 (MQ=255)
                                                     |                                                               
CCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAG  >  W3110S.gb/515881‑515997

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: