Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,489,152 T→C L494L (TTA→TTG pepA ← aminopeptidase A, a cyteinylglycinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,489,1520TC87.4% 244.6 / 21.0 95L494L (TTA→TTGpepAaminopeptidase A, a cyteinylglycinase
Reads supporting (aligned to +/- strand):  ref base T (11/1);  new base C (0/83);  total (11/84)
Fisher's exact test for biased strand distribution p-value = 1.54e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTG  >  W3110S.gb/4489095‑4489191
                                                         |                                       
cGGCAATAACAGCCTTGCCTGACGCATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2431455/62‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCTCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1523657/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAg                                    <  1:1692872/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2449226/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2810406/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2798770/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2794823/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:279076/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2662605/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2637160/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2559324/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2509151/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:245055/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2894867/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:239536/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2372297/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2312408/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2298500/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2291664/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2259208/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2227072/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2187814/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2147665/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:760184/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:95712/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:955389/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:882697/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:866516/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:856125/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:827584/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:810303/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:804861/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:790783/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:776178/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2882901/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:682555/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:570048/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:474496/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:451845/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:421488/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:367860/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:3451/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:335923/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2922632/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2153852/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1353332/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1678003/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1676485/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1671267/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1579208/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1554164/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1510039/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:148679/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:147700/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1454217/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1578785/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1309108/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1288819/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1266166/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1235613/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1132755/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:11194/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1028128/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1026289/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:101266/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2130503/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2128947/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2120708/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2066234/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2063175/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2033434/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2023839/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1686652/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1970070/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1949386/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1795436/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1789401/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1777553/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2133876/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGATTCAACAg                                    <  1:324317/61‑1 (MQ=255)
  gCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCATTAAACCCAGCGCGGTTCAACAg                                    <  1:1890575/61‑1 (MQ=255)
   cAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:1896425/60‑1 (MQ=255)
             cTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:512987/50‑1 (MQ=255)
                         aaTTACTCTTCGCCGTTAAACCCAGCGCGGTTCAACAg                                    <  1:2051195/38‑1 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:137289/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:1043799/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:1107214/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:899700/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:1004983/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:663757/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:497385/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:1424505/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:2144269/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:2854005/1‑62 (MQ=255)
                                   cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTg  >  1:1759457/1‑62 (MQ=255)
                                                         |                                       
CGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTG  >  W3110S.gb/4489095‑4489191

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: