Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 991,514 C→T 8.0% S157S (TCC→TCT pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb991,5140CT 8.0% 153.3 / 3.6 50S157S (TCC→TCTpepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CGGACGTGCTGGCGCGTTTTACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCC  >  W3110S.gb/991447‑991579
                                                                   |                                                                 
cGGACGTGCTGGCGCGTTTTACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAAt                                                                <  1:439502/71‑1 (MQ=255)
cGGACGTGCTGGCGCGTTTTACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAAt                                                                >  1:3569411/1‑71 (MQ=255)
                ttttACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTcc                                                                   <  1:4678307/52‑1 (MQ=255)
                  ttACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCTAATGGTAACCGCGTTgcgc                                                >  1:4183983/1‑69 (MQ=255)
                  ttACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTgcgc                                                >  1:781199/1‑69 (MQ=255)
                    accaccAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTcc                                                                   <  1:3187238/48‑1 (MQ=255)
                      caccaAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGc                                           >  1:1526205/1‑70 (MQ=255)
                      caccaAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCg                                          <  1:4747976/71‑1 (MQ=255)
                       accaAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGa                                         <  1:2304106/71‑1 (MQ=255)
                        ccaAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        <  1:827403/71‑1 (MQ=255)
                         caAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTcc                                                                   <  1:3092357/43‑1 (MQ=255)
                          aaaaTTATTGCCGATAAAATCATATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAAc                                       >  1:184098/1‑70 (MQ=255)
                          aaaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCaa                                                                 >  1:4549941/1‑44 (MQ=255)
                          aaaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        >  1:5106292/1‑69 (MQ=255)
                           aaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCaa                                                                 >  1:4045645/1‑43 (MQ=255)
                           aaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCaa                                                                 >  1:2689634/1‑43 (MQ=255)
                           aaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCaa                                                                 >  1:3694833/1‑43 (MQ=255)
                           aaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCaa                                                                 >  1:4359387/1‑43 (MQ=255)
                           aaaTTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTg                                     >  1:510403/1‑71 (MQ=255)
                              ttattGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGaa                                  <  1:2415926/71‑1 (MQ=255)
                              ttattGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGaa                                  <  1:3167700/71‑1 (MQ=255)
                               tattGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGaa                                  >  1:2978479/1‑70 (MQ=255)
                                attGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGaaa                                 <  1:1574643/70‑1 (MQ=255)
                                   gCCGATAAAATCAAATATCCCTTCCTGCTTTTCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:2986317/70‑1 (MQ=255)
                                   gCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              >  1:5093829/1‑70 (MQ=255)
                                   gCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:3175461/70‑1 (MQ=255)
                                   gCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:1823610/70‑1 (MQ=255)
                                   gCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:4504864/70‑1 (MQ=255)
                                    ccGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGa                            >  1:1297437/1‑71 (MQ=255)
                                     cGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTgc                                                  <  1:527793/48‑1 (MQ=255)
                                     cGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTgc                                                  <  1:2459516/48‑1 (MQ=255)
                                     cGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTgc                                                  <  1:4898561/48‑1 (MQ=255)
                                        tAAAATCAAATATCCCTTCCTGCTTTCTAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGcc                        <  1:2199841/71‑1 (MQ=255)
                                         aaaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        >  1:2688927/1‑54 (MQ=255)
                                         aaaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        >  1:1991034/1‑54 (MQ=255)
                                         aaaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        >  1:4515121/1‑54 (MQ=255)
                                         aaaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGaa                                        >  1:4616007/1‑54 (MQ=255)
                                         aaaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCa                       >  1:174195/1‑71 (MQ=255)
                                           aaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:2281240/62‑1 (MQ=255)
                                            aTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:1989301/61‑1 (MQ=255)
                                            aTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACg                              <  1:3391328/61‑1 (MQ=255)
                                             tCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTgg                   >  1:942775/1‑71 (MQ=255)
                                             tCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTgg                   >  1:463102/1‑71 (MQ=255)
                                              cAAATATCCCTTCCTGCTTTCTAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTggg                  >  1:5103193/1‑71 (MQ=255)
                                                  tatCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACa              <  1:3414881/71‑1 (MQ=255)
                                                  tatCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACa              <  1:1554961/71‑1 (MQ=255)
                                                    tCCCTTCCTGCTTTCTAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGt            <  1:4113888/71‑1 (MQ=255)
                                                      ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTg           <  1:2738084/70‑1 (MQ=255)
                                                        ttCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGaa                                  <  1:131007/45‑1 (MQ=255)
                                                          ccTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAg       >  1:1659682/1‑70 (MQ=255)
                                                           cTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAgg      <  1:598088/70‑1 (MQ=255)
                                                              cTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGAccc  >  1:1361756/1‑71 (MQ=255)
                                                                 tCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACTCCa                       <  1:796941/47‑1 (MQ=255)
                                                                   |                                                                 
CGGACGTGCTGGCGCGTTTTACCACCAAAATTATTGCCGATAAAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCC  >  W3110S.gb/991447‑991579

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: