Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,820,171 G→A 10.0% A166V (GCA→GTA)  alaS ← alanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,820,1710GA10.0% 151.4 / 5.4 50A166V (GCA→GTA) alaSalanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/16);  new base A (2/3);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCC  >  W3110S.gb/2820104‑2820239
                                                                   |                                                                    
aTCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCg                                                                    <  1:4150943/70‑1 (MQ=255)
aTCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCg                                                                    >  1:481108/1‑70 (MQ=255)
  cGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTAt                                                                 >  1:2606611/1‑71 (MQ=255)
  cGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTAt                                                                 >  1:123875/1‑71 (MQ=255)
              aggTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCa                                                           >  1:334842/2‑65 (MQ=255)
                gTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGtt                                                   <  1:3398840/71‑1 (MQ=255)
                gTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                                    >  1:3363297/1‑70 (MQ=255)
                     cGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGc                                              >  1:1048712/1‑71 (MQ=255)
                     cGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGc                                              >  1:2053922/1‑71 (MQ=255)
                     cGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGc                                              >  1:4597038/1‑71 (MQ=255)
                     cGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGc                                              >  1:1357424/1‑71 (MQ=255)
                      gggCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGcc                                             >  1:2345711/1‑71 (MQ=255)
                       ggCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCg                                            <  1:4287274/71‑1 (MQ=255)
                       ggCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCg                                            <  1:2727563/71‑1 (MQ=255)
                         ccGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGAt                                          >  1:1757564/1‑71 (MQ=255)
                            cACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGc                                              <  1:3740950/64‑1 (MQ=255)
                              cGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCg                                                                    <  1:2825138/40‑1 (MQ=255)
                              cGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                                    >  1:1606833/1‑56 (MQ=255)
                              cGGAACAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                                    >  1:1287588/1‑56 (MQ=255)
                                     gtgtCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACg                              >  1:2768045/1‑71 (MQ=255)
                                     gtgtCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACg                              >  1:3213355/1‑71 (MQ=255)
                                      tgtCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGcc                                             <  1:3400683/55‑1 (MQ=255)
                                        tCACCCATCTGCCAGAAGTTGTCAGATACGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGtt                            >  1:5136916/1‑70 (MQ=255)
                                        tCACCCATCTGCCAGAAGTTGTCAGATACGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTc                           >  1:192371/1‑71 (MQ=255)
                                          aCCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcg                          <  1:3677315/70‑1 (MQ=255)
                                          aCCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcg                          <  1:1159228/70‑1 (MQ=255)
                                           cccATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGtt                                                   <  1:1792688/44‑1 (MQ=255)
                                            ccATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTg                                                     <  1:4261450/41‑1 (MQ=255)
                                             cATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcgcgcg                      >  1:910712/1‑71 (MQ=255)
                                              aTCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                                    <  1:4958365/40‑1 (MQ=255)
                                               tCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaataa                                  <  1:5141420/57‑1 (MQ=255)
                                               tCTGCCAGAAGTTGTCAGATACGTATGGCGCACCTTTGtt                                                   <  1:1398721/40‑1 (MQ=255)
                                                 tGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaata                                   >  1:3642126/1‑54 (MQ=255)
                                                 tGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaata                                   >  1:3770356/1‑54 (MQ=255)
                                                 tGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaata                                   >  1:751935/1‑54 (MQ=255)
                                                   ccAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg                    <  1:4020680/67‑1 (MQ=255)
                                                   ccAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg                    <  1:227241/67‑1 (MQ=255)
                                                      gAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaa                                     >  1:1453256/1‑47 (MQ=255)
                                                      gAAGTTGTCAGATGCGTATGGCGCACATTTGTTATCGCCGATGCGaa                                     >  1:5135396/1‑47 (MQ=255)
                                                       aaGTTGTCAGATACGTATGGCGCACCTTTGTTATCGCCGATGCGaataa                                  <  1:2038085/49‑1 (MQ=255)
                                                       aaGTTGTCAGATACGTATGGCGCACCTTTGTTATCGCCGATGCGaataa                                  <  1:2031353/49‑1 (MQ=255)
                                                        aGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcgcgcg                      >  1:1630215/1‑60 (MQ=255)
                                                         gTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTagc           >  1:1066526/1‑68 (MQ=255)
                                                            gTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcg                          >  1:5014403/1‑52 (MQ=255)
                                                            gTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcg                          >  1:4021881/1‑52 (MQ=255)
                                                            gTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcg                          >  1:2609853/1‑52 (MQ=255)
                                                              cAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTActtctt       <  1:482178/69‑1 (MQ=255)
                                                                 aTGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTcc  >  1:4070480/1‑71 (MQ=255)
                                                                 aTGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTcc  >  1:586973/1‑71 (MQ=255)
                                                                 aTGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTcc  >  1:590301/1‑71 (MQ=255)
                                                                   |                                                                    
ATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCC  >  W3110S.gb/2820104‑2820239

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: