Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,865,292 G→A 11.8% 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,2920GA11.8% 144.1 / 9.1 51Q267* (CAG→TAG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (2/4);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC  >  W3110S.gb/2865224‑2865362
                                                                    |                                                                      
gAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                                                      >  1:1604678/1‑71 (MQ=255)
   cAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt                                                                   >  1:3211674/1‑71 (MQ=255)
   cAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                                                  >  1:3798942/1‑71 (MQ=255)
     gcgcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                                                 >  1:348095/1‑71 (MQ=255)
     gcgcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                                                 >  1:1414014/1‑71 (MQ=255)
           cGATATTCAGCCCCTGCTTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                                            <  1:3244086/70‑1 (MQ=255)
           cGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                                            <  1:2857687/70‑1 (MQ=255)
            gATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                                           <  1:2826238/70‑1 (MQ=255)
               atTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                                                       >  1:1409745/1‑71 (MQ=255)
               atTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                                                       >  1:1794262/1‑71 (MQ=255)
                ttCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                                                       <  1:1089986/70‑1 (MQ=255)
                ttCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGtt                                                       <  1:879003/70‑1 (MQ=255)
                ttCAGCACCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGtt                                                       <  1:2217518/70‑1 (MQ=255)
                   aGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:3236030/70‑1 (MQ=255)
                   aGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                                                   <  1:993231/71‑1 (MQ=255)
                         tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                                                       <  1:3865701/61‑1 (MQ=255)
                          gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCCACCTGAATCTGGCGAACACGTTCACGGGTGAgg                                            >  1:1171428/1‑71 (MQ=255)
                          gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                                          >  1:873318/1‑57 (MQ=255)
                          gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                                          >  1:1437550/1‑57 (MQ=255)
                          gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                                            >  1:163491/1‑71 (MQ=255)
                          gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                                            >  1:23529/1‑71 (MQ=255)
                           cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                                           >  1:2123464/1‑71 (MQ=255)
                             tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                                         <  1:3869930/71‑1 (MQ=255)
                                 cAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt                                     <  1:439155/71‑1 (MQ=255)
                                     aTTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg                                 <  1:2058220/71‑1 (MQ=255)
                                         cgcgCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                              <  1:2499627/70‑1 (MQ=255)
                                          gcgcAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa                            <  1:2845170/71‑1 (MQ=255)
                                                aCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:2059464/41‑1 (MQ=255)
                                                aCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:1782436/41‑1 (MQ=255)
                                                aCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:2611580/41‑1 (MQ=255)
                                                aCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:2692801/41‑1 (MQ=255)
                                                aCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                                                    <  1:1466133/41‑1 (MQ=255)
                                                     gcAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAgt                 <  1:1828234/71‑1 (MQ=255)
                                                           cTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc           >  1:1871013/1‑71 (MQ=255)
                                                           cTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc           >  1:3611680/1‑71 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGTGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:2728185/61‑1 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:3015424/61‑1 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:1219743/61‑1 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:805007/61‑1 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:675351/61‑1 (MQ=255)
                                                            ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                    <  1:411966/61‑1 (MQ=255)
                                                              cAAGCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:2423343/1‑71 (MQ=255)
                                                              cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:1509556/1‑71 (MQ=255)
                                                              cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:1137296/1‑71 (MQ=255)
                                                              cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:546129/1‑71 (MQ=255)
                                                              cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:2134454/1‑71 (MQ=255)
                                                              cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt        >  1:224017/1‑71 (MQ=255)
                                                                  cTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt                                     <  1:24339/38‑1 (MQ=255)
                                                                  cTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGtt              <  1:1202644/61‑1 (MQ=255)
                                                                  cTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGtt              <  1:154985/61‑1 (MQ=255)
                                                                   tGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc                     >  1:1965851/1‑53 (MQ=255)
                                                                    gAATCTGGCGAACACGTTGACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAc  >  1:1042275/1‑71 (MQ=255)
                                                                    |                                                                      
GAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC  >  W3110S.gb/2865224‑2865362

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: