Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,813,784 C→T 8.8% A438T (GCC→ACC)  gshA ← gamma‑glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,813,7840CT 8.8% 208.4 / 8.4 68A438T (GCC→ACC) gshAgamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (4/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGC  >  W3110S.gb/2813717‑2813847
                                                                   |                                                               
aTACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAAc                                                              >  1:1232688/1‑71 (MQ=255)
  aCCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAAc                                                              >  1:3480821/1‑69 (MQ=255)
  aCCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAAc                                                              >  1:3700910/1‑69 (MQ=255)
     aGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAg                                                           >  1:499535/1‑69 (MQ=255)
      gTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAAc                                                              <  1:4113313/65‑1 (MQ=255)
       tATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtca                                                       <  1:1393388/71‑1 (MQ=255)
         tcaatcaTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGTAACCAGTtcat                                                      >  1:2410408/1‑70 (MQ=255)
         tcaatcaTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtcat                                                      <  1:2936062/70‑1 (MQ=255)
         tcaatcaTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtcat                                                      <  1:4534969/70‑1 (MQ=255)
             tcaTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtcat                                                      <  1:3491466/66‑1 (MQ=255)
             tcaTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtcat                                                      <  1:1350680/66‑1 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGTAACCAGTTCATcacaca                                                <  1:4118083/70‑1 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGtt                                                         >  1:2253505/1‑61 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtcatca                                                    >  1:898482/1‑66 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCacacac                                               >  1:1571799/1‑71 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCacacac                                               >  1:400307/1‑71 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCacacac                                               >  1:909754/1‑71 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCacacac                                               >  1:915008/1‑71 (MQ=255)
               aTAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCacacac                                               >  1:3083066/1‑71 (MQ=255)
                tAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACt                                              <  1:3265008/71‑1 (MQ=255)
                     cTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtca                                                       <  1:4134527/57‑1 (MQ=255)
                        aaGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtca                                                       >  1:976642/1‑54 (MQ=255)
                        aaGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGa                                        <  1:2840441/69‑1 (MQ=255)
                        aaGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGa                                        <  1:843049/69‑1 (MQ=255)
                        aaGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATa                                      >  1:1335520/1‑71 (MQ=255)
                         aGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATa                                      <  1:354498/70‑1 (MQ=255)
                            tACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAg                                                           >  1:1259528/1‑46 (MQ=255)
                            tACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATAACACACTTTCTGATACGCt                                  <  1:3655464/71‑1 (MQ=255)
                             aCGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCtt                                 >  1:4462133/1‑71 (MQ=255)
                             aCGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCtt                                 >  1:1320935/1‑71 (MQ=255)
                             aCGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCtt                                 >  1:3242440/1‑71 (MQ=255)
                              cGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTc                                <  1:4398829/71‑1 (MQ=255)
                                  cAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTcgccg                            >  1:951558/1‑71 (MQ=255)
                                  cAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCCGATACGCTTcgccg                            >  1:3348277/1‑71 (MQ=255)
                                   agagAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATcac                                                   <  1:465419/47‑1 (MQ=255)
                                    gagaAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTcgcc                             <  1:1189033/68‑1 (MQ=255)
                                     agaAAGTCAGATCGGGATTATCGAAGCAGGCAACCa                                                            >  1:1760953/1‑36 (MQ=255)
                                     agaAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGtt                                                         <  1:1164829/39‑1 (MQ=255)
                                     agaAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGtt                                                         <  1:547020/39‑1 (MQ=255)
                                     agaAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGtt                                                         <  1:2843124/39‑1 (MQ=255)
                                      gaAAGTCAGATCGGGATTATCGAAGCAGGTAACCAGTTCATCACACACTTTCTGATACTCTTCGCCGCCGt                        <  1:1751594/71‑1 (MQ=255)
                                       aaaGTCAGATCGGGATTATCGAAGCAGGCAACCAGtt                                                         <  1:1350372/37‑1 (MQ=255)
                                       aaaGTCAGATCGGGATTATCGAAGCAGGCAACCAGTtc                                                        <  1:3764365/38‑1 (MQ=255)
                                       aaaGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTTTGATACGCTTCGCCGCCGt                        >  1:3771038/1‑70 (MQ=255)
                                       aaaGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATa                                      >  1:4076003/1‑56 (MQ=255)
                                       aaaGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGtt                       <  1:2515395/71‑1 (MQ=255)
                                         aGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTc                                           <  1:2138887/49‑1 (MQ=255)
                                          gTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAAt                    <  1:3579454/71‑1 (MQ=255)
                                           tCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATa                   >  1:869917/1‑71 (MQ=255)
                                             aGATCGGGATTATCGAAGCAGGTAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATAc                  >  1:380225/1‑70 (MQ=255)
                                             aGATCGGGATTATCGAAGCAGGTAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATAc                  >  1:3059866/1‑70 (MQ=255)
                                              gATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCt                                          <  1:1305139/45‑1 (MQ=255)
                                                tCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGa                                        <  1:218267/45‑1 (MQ=255)
                                                tCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTaa                     <  1:1523383/64‑1 (MQ=255)
                                                 cGGGATTATCGAAGCAGGTAACCAGTTCATCACACACTTTCTGATACGc                                   >  1:2842030/1‑49 (MQ=255)
                                                 cGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATAc                                     <  1:454538/47‑1 (MQ=255)
                                                      ttATCGAAGCAGGCAACCAGTTCATCACACACttt                                            <  1:2430396/35‑1 (MQ=255)
                                                      ttATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGt        >  1:2209902/1‑71 (MQ=255)
                                                       tATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGtt       >  1:1518235/1‑71 (MQ=255)
                                                       tATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGtt       >  1:1512807/1‑71 (MQ=255)
                                                        aTCGAAGCAGGCAACCAGTTCATCACACACTTTCTGa                                        >  1:4230142/1‑37 (MQ=255)
                                                        aTCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGttt      >  1:3847420/1‑71 (MQ=255)
                                                         tCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTg     <  1:3417979/71‑1 (MQ=255)
                                                           gAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTa                      >  1:3822886/1‑52 (MQ=255)
                                                            aaGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTgcgc  <  1:817653/71‑1 (MQ=255)
                                                            aaGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTgcgc  <  1:2716216/71‑1 (MQ=255)
                                                            aaGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTgcgc  <  1:3216633/71‑1 (MQ=255)
                                                                  gCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGtt                       <  1:593844/44‑1 (MQ=255)
                                                                   |                                                               
ATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCAGGCAACCAGTTCATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGC  >  W3110S.gb/2813717‑2813847

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: