Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,331,670 C→T 10.8% Q298* (CAA→TAA)  atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,331,6700CT10.8% 136.4 / 6.1 46Q298* (CAA→TAA) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/3);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGC  >  W3110S.gb/2331604‑2331740
                                                                  |                                                                      
aTTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCaaaa                                                                    <  1:827875/71‑1 (MQ=255)
aTTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCaaa                                                                     >  1:4020536/1‑70 (MQ=255)
  tAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg                                                                  >  1:522254/1‑71 (MQ=255)
   aaaaGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc                                                                 >  1:4042896/1‑71 (MQ=255)
    aaaGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc                                                                 <  1:3739160/70‑1 (MQ=255)
        ttATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTa                                                             >  1:4896537/1‑70 (MQ=255)
        ttATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAc                                                            <  1:207413/71‑1 (MQ=255)
         tATGCCAGCGGTGGCGTGCCCCCCGCATTATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACa                                                           <  1:4464146/70‑1 (MQ=255)
                gCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCg                                                    <  1:4474953/71‑1 (MQ=255)
                        gTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCa                                            <  1:4365173/71‑1 (MQ=255)
                           ccccccGCATTGATGGGTATGGGGCCAGTACCTGCCACGTAAAAAGCGTTACaa                                                          >  1:3055166/1‑54 (MQ=255)
                              cccGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCa                                            <  1:621371/65‑1 (MQ=255)
                               ccGCATTGATGGGTATGGGGCCAGTACCTGCCACGTAAAAAGCGtt                                                              <  1:2765323/46‑1 (MQ=255)
                                cGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCg                                     >  1:1308469/1‑70 (MQ=255)
                                 gCATTGATGGGTATGGGGCCAGTACCTGCCACGTAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGa                                   <  1:3116342/71‑1 (MQ=255)
                                 gCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGa                                   >  1:86422/1‑71 (MQ=255)
                                 gCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGa                                   <  1:1646750/71‑1 (MQ=255)
                                  cATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGc                                             <  1:1681882/60‑1 (MQ=255)
                                    ttGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc                                                                 <  1:3592267/38‑1 (MQ=255)
                                    ttGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc                                                                 <  1:1049759/38‑1 (MQ=255)
                                      gATGGGTATGGGGCCAGTACCTGCCACGTAAAAAGCGTTa                                                             >  1:3639632/1‑40 (MQ=255)
                                      gATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGGGGATAtt                               >  1:1196177/1‑70 (MQ=255)
                                      gATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATAtt                               >  1:4584138/1‑70 (MQ=255)
                                      gATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATAtt                               >  1:1745821/1‑70 (MQ=255)
                                        tGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAt                            >  1:4510426/1‑71 (MQ=255)
                                           gTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAtctc                         <  1:2509293/71‑1 (MQ=255)
                                           gTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAtctc                         <  1:249195/71‑1 (MQ=255)
                                              tGGGGCCAGTACCTGCCACGCAAAAAGCGTTACaa                                                          <  1:3668247/35‑1 (MQ=255)
                                              tGGGGCCAGTACCTGCCACGCAAAAAGCGTTACaa                                                          <  1:5015617/35‑1 (MQ=255)
                                              tGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCAtt                      >  1:2762303/1‑71 (MQ=255)
                                                      gTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAtct                          <  1:4035522/59‑1 (MQ=255)
                                                      gTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTaa              >  1:712674/1‑71 (MQ=255)
                                                      gTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTa               >  1:127476/1‑70 (MQ=255)
                                                       tACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTg                     <  1:836280/63‑1 (MQ=255)
                                                       tACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAAt             >  1:1454443/1‑71 (MQ=255)
                                                        aCCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAAt             >  1:4863067/1‑70 (MQ=255)
                                                         ccTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAt                            >  1:3782140/1‑54 (MQ=255)
                                                         ccTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAt                            >  1:4183769/1‑54 (MQ=255)
                                                         ccTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAt                            >  1:3995058/1‑54 (MQ=255)
                                                          cTGCCACGTAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGaa          <  1:5057635/71‑1 (MQ=255)
                                                             ccACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCa       <  1:4801056/71‑1 (MQ=255)
                                                             ccACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCa       >  1:1142389/1‑71 (MQ=255)
                                                             ccACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAAGGAAGCa       <  1:4926/71‑1 (MQ=255)
                                                              cACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCAt      >  1:3005589/1‑71 (MQ=255)
                                                                cGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCAttt    <  1:1821021/71‑1 (MQ=255)
                                                                 gCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTtg   >  1:3067086/1‑71 (MQ=255)
                                                                  cAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTtgc  >  1:4428151/1‑71 (MQ=255)
                                                                  |                                                                      
ATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGC  >  W3110S.gb/2331604‑2331740

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: