Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,363,356 G→A S232N (AGC→AAC)  yhcD → predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,363,3560GA88.3% 265.9 / 22.6 94S232N (AGC→AAC) yhcDpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/0);  new base A (26/57);  total (37/57)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGC  >  W3110S.gb/3363297‑3363420
                                                           |                                                                
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTTGAAg                                                       <  1:489141/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:3777406/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:2253428/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:2058155/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:2053307/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:1519362/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:1482206/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:1399397/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:3148614/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:1131152/71‑1 (MQ=255)
aCGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       <  1:3768342/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:3174290/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:2583056/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:1994145/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:2583306/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:2837631/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:3757946/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:2968644/71‑1 (MQ=255)
   tACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAg                                                    <  1:4121644/71‑1 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:1575041/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:2218931/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:1101242/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:2977759/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:3455660/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:1858460/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:183501/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:468266/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:594818/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:275712/1‑70 (MQ=255)
     ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGc                                                   >  1:2855335/1‑70 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:2705174/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:3697611/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:46691/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:937656/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:3152605/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:8549/1‑56 (MQ=255)
               tCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAg                                                       >  1:1737123/1‑56 (MQ=255)
                           cATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTgg                                                          >  1:3001594/1‑41 (MQ=255)
                           cATGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTgg                                                          >  1:3612340/1‑41 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACg                                        >  1:686869/1‑58 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACg                                        >  1:1668290/1‑58 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACg                                        >  1:718918/1‑58 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGc                           <  1:1426275/71‑1 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGc                           <  1:731692/71‑1 (MQ=255)
                            aTGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGc                           <  1:1228249/71‑1 (MQ=255)
                             ttcGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCa                          <  1:3273233/69‑1 (MQ=255)
                             tGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCa                          <  1:2740758/71‑1 (MQ=255)
                             tGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCa                          <  1:2735967/71‑1 (MQ=255)
                             tGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCa                          <  1:1805445/71‑1 (MQ=255)
                             tGCGAGCAGGTGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCa                          <  1:3777819/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:4023841/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:3086695/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:2342708/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:3824159/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:2531903/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:2793470/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:2788871/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:2783106/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:1781120/71‑1 (MQ=255)
                                       tGAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                <  1:591946/71‑1 (MQ=255)
                                        gAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                >  1:1287224/1‑70 (MQ=255)
                                        gAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                >  1:3154164/1‑70 (MQ=255)
                                        gAAGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACg                >  1:1793773/1‑70 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1426611/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1202398/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1948419/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1203283/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:2691861/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:2740113/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1358991/61‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAt             <  1:3865157/71‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAt             <  1:3422231/71‑1 (MQ=255)
                                          aGTTAACCTCAATAATAACGTGTTGAAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1089527/61‑1 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:3801585/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:1030814/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:858552/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:1193626/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:121609/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:2677011/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:1318219/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:2667976/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:496622/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCAtt            >  1:1969054/1‑71 (MQ=255)
                                           gTTAACCTCAATAATAGCGTGTAGGAAGGAGCCAGTATTTACGGTATCGAAATCGCCCCgg                      >  1:3549018/1‑61 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:3734301/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1491269/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:606182/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:1783927/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:474289/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:2884986/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:2909502/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:321940/51‑1 (MQ=255)
                                                    aataatAACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCg                       <  1:3985385/51‑1 (MQ=255)
                                                     ataataACGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGc  <  1:970987/71‑1 (MQ=255)
                                                           |                                                                
ACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGC  >  W3110S.gb/3363297‑3363420

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