Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 199,072 T→C 16.0% V382A (GTC→GCC)  yaeT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb199,0720TC16.0% 108.0 / 16.3 51V382A (GTC→GCC) yaeTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/27);  new base C (6/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 5.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

CTCGAAAGATGCCGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTCGATACCGATACCCAACG  >  W3110S.gb/199004‑199138
                                                                    |                                                                  
cTCGAAAGATGCCGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTc                                                                   <  1:341943/70‑1 (MQ=255)
    aaaGATGCCGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCa                                                              <  1:630902/71‑1 (MQ=255)
        aTGCCGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGt                                                          <  1:2746843/71‑1 (MQ=255)
           ccGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGAt                                                                >  1:692942/1‑62 (MQ=255)
           ccGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAg                                                       <  1:1549797/71‑1 (MQ=255)
            cGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAgg                                                      <  1:1236707/71‑1 (MQ=255)
              tCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGt                                                          <  1:615861/65‑1 (MQ=255)
                cTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGc                                                   >  1:351111/1‑70 (MQ=255)
                cTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAGGAGc                                                   <  1:1014915/70‑1 (MQ=255)
                  gCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTc                                                >  1:142894/1‑71 (MQ=255)
                   cGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCt                                               >  1:997992/1‑71 (MQ=255)
                   cGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCt                                               >  1:775837/1‑71 (MQ=255)
                   cGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCt                                               >  1:635483/1‑71 (MQ=255)
                     tCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCg                                                                  <  1:658288/50‑1 (MQ=255)
                     tCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCg                                                                  <  1:2182903/50‑1 (MQ=255)
                     tCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCg                                                                  <  1:1791551/50‑1 (MQ=255)
                      cgcgAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATAAGGGTAAGGAGCGTCTGa                                             <  1:987110/70‑1 (MQ=255)
                       gcgAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAg                                                       >  1:1486653/1‑59 (MQ=255)
                       gcgAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAg                                                       >  1:1389963/1‑59 (MQ=255)
                       gcgAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAg                                                       >  1:1846735/1‑59 (MQ=255)
                       gcgAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAGGAGCGTCTGAAt                                           >  1:1832274/1‑71 (MQ=255)
                          aaaTGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGt                                        <  1:633692/71‑1 (MQ=255)
                           aaTGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGt                                        <  1:1185944/70‑1 (MQ=255)
                           aaTGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTc                                       >  1:936092/1‑71 (MQ=255)
                           aaTGCGCCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTc                                       >  1:2235428/1‑71 (MQ=255)
                                gTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGCCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGc                                  <  1:1073176/71‑1 (MQ=255)
                                 tCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGc                                  <  1:874746/70‑1 (MQ=255)
                                 tCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGGCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGct                                 <  1:1536932/71‑1 (MQ=255)
                                   aGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAAt                                           >  1:2729967/1‑59 (MQ=255)
                                        gAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGa                          <  1:2782880/71‑1 (MQ=255)
                                          aGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGcttc                               <  1:1053583/64‑1 (MQ=255)
                                           ggTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTc                                       <  1:892829/55‑1 (MQ=255)
                                           ggTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTc                                       <  1:520520/55‑1 (MQ=255)
                                           ggTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTc                                       <  1:1949260/55‑1 (MQ=255)
                                             tGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCCTCTGAATCGTc                                       >  1:22229/1‑53 (MQ=255)
                                                 tGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATc                                          <  1:807355/46‑1 (MQ=255)
                                                 tGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGaa                         <  1:1353376/63‑1 (MQ=255)
                                                       gggAGCGATCTTGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgatacc             <  1:1691069/69‑1 (MQ=255)
                                                       gggAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgatacc             <  1:1624295/69‑1 (MQ=255)
                                                        ggAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGa                          <  1:183500/55‑1 (MQ=255)
                                                        ggAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgataccgat          <  1:1224226/71‑1 (MQ=255)
                                                          aGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgataccgatac        >  1:1507513/1‑71 (MQ=255)
                                                          aGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgataccgatac        >  1:1208380/1‑71 (MQ=255)
                                                          aGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgataccgatac        >  1:2560803/1‑71 (MQ=255)
                                                          aGCGATCTGGCCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgataccgatac        >  1:1308678/1‑71 (MQ=255)
                                                            cGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGa                          <  1:593387/51‑1 (MQ=255)
                                                              aTCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTcgatacc             <  1:1957699/62‑1 (MQ=255)
                                                                cTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTCGATACCGATACCCAACg  >  1:2516839/1‑71 (MQ=255)
                                                                cTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTCGATACCGATACCCAACg  >  1:1257749/1‑71 (MQ=255)
                                                                cTGGCCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTCGATACCGATACCCAACg  >  1:1149489/1‑71 (MQ=255)
                                                                   gTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTc                   <  1:1991721/51‑1 (MQ=255)
                                                                    |                                                                  
CTCGAAAGATGCCGTCCTGCGTCGCGAAATGCGTCAGATGGAAGGTGCATGGCTGGGGAGCGATCTGGTCGATCAGGGTAAGGAGCGTCTGAATCGTCTGGGCTTCTTTGAAACTGTCGATACCGATACCCAACG  >  W3110S.gb/199004‑199138

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: