Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,951,622 A→G 7.4% G205G (GGT→GGC aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,951,6220AG 7.4% 142.3 / 2.8 54G205G (GGT→GGCaspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

TCAGGCTGACGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTTGTG  >  W3110S.gb/1951554‑1951690
                                                                    |                                                                    
tCAGGCTGCCGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAAcc                                                                    >  1:2959886/1‑71 (MQ=255)
tCAGGCTGACGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAAAGATCTGATAGTAACGGTCAAAAcc                                                                    >  1:3588744/1‑71 (MQ=255)
  aGGCTGACGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCgg                                                                  <  1:2144134/71‑1 (MQ=255)
  aGGCTGACGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCgg                                                                  <  1:1696917/71‑1 (MQ=255)
        aCGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAGCGGTCAAAGCCGGAcatca                                                            <  1:3282497/71‑1 (MQ=255)
        aCGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAGCGGTCAAAGCCGGAcatca                                                            <  1:1500101/71‑1 (MQ=255)
             cAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATcag                                                        >  1:1184182/1‑70 (MQ=255)
                 aCGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATcagca                                                      >  1:1868267/1‑68 (MQ=255)
                 aCGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATcagca                                                      >  1:327824/1‑68 (MQ=255)
                  cGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTg                                                  >  1:1885410/1‑71 (MQ=255)
                       gTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATcagca                                                      <  1:1289803/62‑1 (MQ=255)
                       gTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTg                                              <  1:3350590/70‑1 (MQ=255)
                       gTCCTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTg                                              <  1:17902/70‑1 (MQ=255)
                        tCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGa                                             <  1:2553641/70‑1 (MQ=255)
                        tCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGa                                             <  1:703512/70‑1 (MQ=255)
                         cTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCagcagc                                                    >  1:1087324/1‑62 (MQ=255)
                            cGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGCTTGaac                                           <  1:3655000/68‑1 (MQ=255)
                              tCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTg                                      <  1:134974/71‑1 (MQ=255)
                              tCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTg                                      <  1:1566291/71‑1 (MQ=255)
                                aCGGAAGCATTTAACGATCTGATAGTAGCGGTCAAAGCCGGACATCATCAGTAGCTGTTTGAACAACTGCg                                    >  1:3449378/1‑71 (MQ=255)
                                aCGGAAGCATTTAACGATCTGATAGTAGCGGTCAAAGCCGGACATCATCAGTAGCTGTTTGAACAACTGCg                                    >  1:1543137/1‑71 (MQ=255)
                                  ggAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCggg                                  <  1:3163686/71‑1 (MQ=255)
                                  ggAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCggg                                  >  1:3045631/1‑71 (MQ=255)
                                  ggAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGCAcatc                                                             >  1:1540613/1‑44 (MQ=255)
                                   gAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTg                                              <  1:380388/58‑1 (MQ=255)
                                       cATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTg                             <  1:1713782/71‑1 (MQ=255)
                                       cATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTg                             <  1:2690133/71‑1 (MQ=255)
                                           tAACGATCTGATAGTAACGGTCAAAACCGGAcatcat                                                           >  1:1718474/1‑37 (MQ=255)
                                           tAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGaac                                           >  1:869795/1‑53 (MQ=255)
                                           tAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCgg                          >  1:859441/1‑70 (MQ=255)
                                            aaCGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGaac                                           <  1:1809587/52‑1 (MQ=255)
                                            aaCGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGc                         <  1:632415/70‑1 (MQ=255)
                                            aaCGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGc                         <  1:953202/70‑1 (MQ=255)
                                                 tCTGATAGTAACGGTCAGAACCGGACATCATCAGCAGCTGTTTGaacaa                                         >  1:3224213/1‑49 (MQ=255)
                                                 tCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGaacaa                                         >  1:3376178/1‑49 (MQ=255)
                                                 tCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGaacaa                                         >  1:1103062/1‑49 (MQ=255)
                                                 tCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCg                   >  1:3344347/1‑71 (MQ=255)
                                                   tGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTg                             <  1:3404542/59‑1 (MQ=255)
                                                   tGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTg                             <  1:3083209/59‑1 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGaa                                            >  1:2451484/1‑42 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGc                    <  1:3464892/66‑1 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAg                >  1:754011/1‑70 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGa               <  1:327779/71‑1 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGa               <  1:1038855/71‑1 (MQ=255)
                                                     aTAGTAACGGTCAAAACCGGACATCATAAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAg                <  1:530245/70‑1 (MQ=255)
                                                         tAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAAttt           <  1:985050/71‑1 (MQ=255)
                                                          aaCGGTCAAAACCGGACATCATCAGCAGCTGTTTGaacaa                                         >  1:3276041/1‑40 (MQ=255)
                                                           aCGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAg                       >  1:2780313/1‑57 (MQ=255)
                                                           aCGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTAc         >  1:2044245/1‑71 (MQ=255)
                                                             ggTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCt       >  1:2643383/1‑71 (MQ=255)
                                                              gTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCtt      >  1:2257913/1‑71 (MQ=255)
                                                              gTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCtt      <  1:122414/71‑1 (MQ=255)
                                                                cAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtg    >  1:2734154/1‑71 (MQ=255)
                                                                cAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtg    >  1:1247719/1‑71 (MQ=255)
                                                                 aaaaCCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtg    <  1:2736635/70‑1 (MQ=255)
                                                                 aaaaCCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtg    <  1:114938/70‑1 (MQ=255)
                                                                  aaaCCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtgtg  <  1:3791205/71‑1 (MQ=255)
                                                                  aaaCCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTtgtg  <  1:3211633/71‑1 (MQ=255)
                                                                  aaaCCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACATTtgtg  >  1:1081704/1‑71 (MQ=255)
                                                                    |                                                                    
TCAGGCTGACGGTCAGCACGCAGGTCTTCGTCACGGAAGCATTTAACGATCTGATAGTAACGGTCAAAACCGGACATCATCAGCAGCTGTTTGAACAACTGCGGGGATTGCGGCAGTGCGTAGAATTTACCTTTGTG  >  W3110S.gb/1951554‑1951690

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: