Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,117,343 C→A 10.3% A142A (GCG→GCT yghJ ← predicted inner membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,117,3430CA10.3% 244.5 / 8.7 77A142A (GCG→GCTyghJpredicted inner membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base A (5/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACGTCCTCGCCAG  >  W3110S.gb/3117276‑3117413
                                                                   |                                                                      
cGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTAcgcg                                                                       >  1:3132032/1‑69 (MQ=255)
cGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTAcgcg                                                                       >  1:3159561/1‑69 (MQ=255)
cGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTAcgcg                                                                       >  1:1642687/1‑69 (MQ=255)
cGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACgcgc                                                                      <  1:2179527/70‑1 (MQ=255)
 ggAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCa                                                                     <  1:2103845/70‑1 (MQ=255)
 ggAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCa                                                                     <  1:1929224/70‑1 (MQ=255)
   aGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCt                                                                  <  1:1037500/71‑1 (MQ=255)
      ccGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTct                                                               <  1:184178/71‑1 (MQ=255)
          cAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctgact                                                           >  1:1164630/1‑71 (MQ=255)
          cAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTctgact                                                           >  1:2681216/1‑71 (MQ=255)
          cAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTctgact                                                           >  1:2687247/1‑71 (MQ=255)
           aGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCtt                                                                 <  1:2425027/64‑1 (MQ=255)
           aGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTct                                                               <  1:474549/66‑1 (MQ=255)
           aGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTct                                                               <  1:3878672/66‑1 (MQ=255)
                 tGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCTCACGCAAGCTACGCGCAGCTTCTGACTGAGTGtt                                                    >  1:129481/1‑71 (MQ=255)
                  gcgcGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGtt                                                    <  1:4109102/70‑1 (MQ=255)
                    gcgTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCt                                                                  <  1:2929296/54‑1 (MQ=255)
                    gcgTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCt                                                                  <  1:5051886/54‑1 (MQ=255)
                      gTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctga                                                             >  1:4286500/1‑57 (MQ=255)
                      gTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAgg                                               <  1:4314689/71‑1 (MQ=255)
                       tCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCt                                                                  >  1:1838712/1‑51 (MQ=255)
                         cTCAAGGCTAAACGACACTTTTTCAACCGCAGGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTgg                                            >  1:2522042/1‑71 (MQ=255)
                         cTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTgg                                            >  1:1793256/1‑71 (MQ=255)
                         cTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTgg                                            >  1:54626/1‑71 (MQ=255)
                         cTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTgg                                            >  1:3264237/1‑71 (MQ=255)
                         cTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTgg                                            >  1:1310914/1‑71 (MQ=255)
                            aaGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctga                                                             >  1:3264057/1‑51 (MQ=255)
                            aaGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCaa                                         >  1:893197/1‑71 (MQ=255)
                            aaGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCa                                          <  1:4002120/70‑1 (MQ=255)
                            aaGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCa                                          <  1:2587953/70‑1 (MQ=255)
                             aGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAAt                                        <  1:2803844/71‑1 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCtt                                                                 >  1:3330396/1‑45 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCtt                                                                 >  1:962875/1‑45 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCtt                                                                 >  1:2117993/1‑45 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAg                                                >  1:753135/1‑62 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAg                                                >  1:2175316/1‑62 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAAtt                                       <  1:4249457/71‑1 (MQ=255)
                              ggCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAAtt                                       <  1:1256413/71‑1 (MQ=255)
                               gCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctgact                                                           <  1:4959157/50‑1 (MQ=255)
                               gCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctgact                                                           <  1:2336983/50‑1 (MQ=255)
                                cTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTctgact                                                           <  1:1307889/49‑1 (MQ=255)
                                  aaaCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCt                                                                  <  1:4555142/40‑1 (MQ=255)
                                  aaaCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCg                                   >  1:1506207/1‑71 (MQ=255)
                                  aaaCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCg                                   >  1:3113525/1‑71 (MQ=255)
                                   aaCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTc                                    <  1:802082/69‑1 (MQ=255)
                                   aaCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgt                                  >  1:61702/1‑71 (MQ=255)
                                   aaCGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCg                                   <  1:2203978/70‑1 (MQ=255)
                                    aCGACACTTTTTCAACCGCACGCAAGCTACGCGCa                                                                     <  1:2214393/35‑1 (MQ=255)
                                    aCGACACTTTTTCAACCGCACGCAAGCTACGCGCa                                                                     <  1:1840904/35‑1 (MQ=255)
                                    aCGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgtg                                 <  1:2281087/71‑1 (MQ=255)
                                    aCGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTctga                                                             >  1:750112/1‑43 (MQ=255)
                                     cGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgtgt                                <  1:2852812/71‑1 (MQ=255)
                                     cGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgtgt                                <  1:2846967/71‑1 (MQ=255)
                                     cGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTCTGACTGAGTGTTGAAGGTg                                             >  1:2247250/1‑58 (MQ=255)
                                     cGACACTTTTTCAACCGCACGCAAGCTACGAGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgtgt                                >  1:3466739/1‑71 (MQ=255)
                                      gACACTTTTTAAACCGCACGCAAGCTACGCGCAGCTTctgactga                                                         >  1:3780330/1‑45 (MQ=255)
                                        cacTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTAc                             >  1:4766375/1‑71 (MQ=255)
                                        cacTTTTTCAACCGCACGAAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTAc                             >  1:747647/1‑71 (MQ=255)
                                          cTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGa                                                  <  1:3741215/48‑1 (MQ=255)
                                          cTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTAccc                           >  1:507689/1‑71 (MQ=255)
                                          cTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTAccc                           >  1:4109190/1‑71 (MQ=255)
                                           tttttCAACCTCAAGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGa                                                  <  1:2388885/47‑1 (MQ=255)
                                             tttCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAgg                                               >  1:1452440/1‑48 (MQ=255)
                                             tttCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGc                         >  1:2331076/1‑70 (MQ=255)
                                             tttCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGc                         >  1:3235646/1‑70 (MQ=255)
                                              ttCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCa                       >  1:4730354/1‑71 (MQ=255)
                                              ttCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCa                       >  1:1602396/1‑71 (MQ=255)
                                               tCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGcc                        <  1:4861526/69‑1 (MQ=255)
                                                 aaCCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGc                    <  1:2615660/71‑1 (MQ=255)
                                                 aaCCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGc                    <  1:2989644/71‑1 (MQ=255)
                                                     gcacgcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGt                <  1:468085/71‑1 (MQ=255)
                                                     gcacgcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGt                <  1:713555/71‑1 (MQ=255)
                                                     gcacgcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGt                <  1:1186060/71‑1 (MQ=255)
                                                      cacgcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCaa                  >  1:1877180/1‑68 (MQ=255)
                                                      cacgcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCAACGCAAGt                >  1:1034187/1‑70 (MQ=255)
                                                      cacgcaAGCTACGAGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCgtgt                                <  1:838291/54‑1 (MQ=255)
                                                         gcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACg            <  1:4043431/71‑1 (MQ=255)
                                                         gcaAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACg            <  1:3540344/71‑1 (MQ=255)
                                                          caAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCa                                          <  1:2955072/40‑1 (MQ=255)
                                                           aaGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACGTc          <  1:955001/71‑1 (MQ=255)
                                                               tACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACGTCCTCg      >  1:361286/1‑71 (MQ=255)
                                                                  gcgcAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTaa              >  1:671381/1‑60 (MQ=255)
                                                                   cgcAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACGTCCTCGCCAg  >  1:2640220/1‑71 (MQ=255)
                                                                   |                                                                      
CGGAGCCCGCCAGTTCTTGCGCGTCCTCAAGGCTAAACGACACTTTTTCAACCGCACGCAAGCTACGCGCAGCTTCTGACTGAGTGTTGAAGGTGGCAATTGTCGTGTTACCCGCCACGCAAGTAACGTCCTCGCCAG  >  W3110S.gb/3117276‑3117413

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: