Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,865,292 G→A 19.6% Q267* (CAG→TAG)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,865,2920GA19.6% 123.6 / 23.7 56Q267* (CAG→TAG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (4/7);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC  >  W3110S.gb/2865229‑2865362
                                                               |                                                                      
gcgcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                                                  >  1:4375915/1‑70 (MQ=255)
  gcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                                               <  1:1527983/71‑1 (MQ=255)
       gATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAacac                                                          <  1:2524628/71‑1 (MQ=255)
        atatTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                                         <  1:3113269/71‑1 (MQ=255)
              aGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCAc                                                    <  1:1767045/70‑1 (MQ=255)
               gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                                                  >  1:1083917/1‑71 (MQ=255)
                   cTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                                              <  1:377133/71‑1 (MQ=255)
                    tGCGTTTGCAGGATTTCGCGCAAACGGCGAAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                                             >  1:1777521/1‑71 (MQ=255)
                     gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                                                     >  1:2924433/1‑46 (MQ=255)
                     gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                                            >  1:2217424/1‑71 (MQ=255)
                     gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                                             >  1:304162/1‑70 (MQ=255)
                     gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                                             >  1:1786782/1‑70 (MQ=255)
                      cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                                                     >  1:3918867/1‑61 (MQ=255)
                      cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                                                     >  1:1539130/1‑61 (MQ=255)
                      cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                                                     >  1:3596099/1‑61 (MQ=255)
                        tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                                         <  1:450880/71‑1 (MQ=255)
                        tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGcc                                          >  1:2666497/1‑70 (MQ=255)
                        tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAAAACGTTCACGGGTGAGGcc                                          >  1:392943/1‑70 (MQ=255)
                           gCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                                        >  1:1043569/1‑69 (MQ=255)
                            cAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt                                     >  1:114975/1‑71 (MQ=255)
                            cAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt                                     <  1:4043481/71‑1 (MQ=255)
                              tgATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                                                 <  1:1878271/56‑1 (MQ=255)
                               gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:2141719/53‑1 (MQ=255)
                               gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                                                <  1:1097314/57‑1 (MQ=255)
                               gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                                               <  1:234287/58‑1 (MQ=255)
                                 tttCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                                                   >  1:595327/1‑52 (MQ=255)
                                   tCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                                                     >  1:4237505/1‑48 (MQ=255)
                                   tCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                              >  1:3873890/1‑71 (MQ=255)
                                    cgcgCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                              >  1:3452134/1‑70 (MQ=255)
                                      cgcAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc                           <  1:2378390/71‑1 (MQ=255)
                                       gcAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa                          <  1:3875656/71‑1 (MQ=255)
                                       gcAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa                          <  1:3383390/71‑1 (MQ=255)
                                          aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                                          >  1:4363206/1‑36 (MQ=255)
                                          aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt                                     <  1:55978/57‑1 (MQ=255)
                                          aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCt                       >  1:3369568/1‑71 (MQ=255)
                                          aaCGGCGCAGGCCTTCAACCTGAATCTG‑‑GAACACGTTCACGGGTGAGGcc                                          <  1:3198942/50‑1 (MQ=255)
                                             ggCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc                           <  1:4275449/64‑1 (MQ=255)
                                             ggCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:3960717/71‑1 (MQ=255)
                                              gcgcAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCa                   <  1:136746/71‑1 (MQ=255)
                                                gcAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAgt                 <  1:3269529/71‑1 (MQ=255)
                                                    gCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa                                >  1:1674942/1‑52 (MQ=255)
                                                     ccTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc            >  1:1395355/1‑71 (MQ=255)
                                                      cTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc           >  1:1320624/1‑71 (MQ=255)
                                                      cTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    >  1:231901/1‑62 (MQ=255)
                                                       ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                  <  1:3599599/47‑1 (MQ=255)
                                                       ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                  <  1:3131810/47‑1 (MQ=255)
                                                       ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:4005885/61‑1 (MQ=255)
                                                        tCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc                     >  1:3927816/1‑59 (MQ=255)
                                                        tCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGtt              <  1:1227117/66‑1 (MQ=255)
                                                         cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa                                   >  1:546218/1‑44 (MQ=255)
                                                         cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:2002739/1‑71 (MQ=255)
                                                         cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:1618751/1‑71 (MQ=255)
                                                         cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:4534652/1‑71 (MQ=255)
                                                          aaCCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCtt       <  1:1753398/71‑1 (MQ=255)
                                                          aaCCTAAATCTGGCGAACACCTTCACGGGTGAGGCCAATTTCa                                   >  1:532393/1‑43 (MQ=38)
                                                             cTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAgt                 <  1:1535188/58‑1 (MQ=255)
                                                               gAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCgatc  >  1:2922923/1‑68 (MQ=255)
                                                               gAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCgatc  >  1:882879/1‑68 (MQ=255)
                                                               |                                                                      
GCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC  >  W3110S.gb/2865229‑2865362

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: