Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,044,904 T→C 20.4% E65E (GAA→GAG etp ← phosphotyrosine‑protein phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,044,9040TC20.4% 111.3 / 25.3 54E65E (GAA→GAGetpphosphotyrosine‑protein phosphatase
Reads supporting (aligned to +/- strand):  ref base T (25/18);  new base C (5/6);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ACTCCATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAC  >  W3110S.gb/1044837‑1044971
                                                                   |                                                                   
aCTCCATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCa                                                                  <  1:2127862/71‑1 (MQ=255)
  tCCATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTc                                                                   >  1:3095725/1‑68 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTTCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               <  1:3126022/70‑1 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCaa                                                                 >  1:1668616/1‑68 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCaa                                                                 >  1:897910/1‑68 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               <  1:2878714/70‑1 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               <  1:4042449/70‑1 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               <  1:313781/70‑1 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCaa                                                                 >  1:2753231/1‑68 (MQ=255)
    cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCaa                                                                 >  1:3637037/1‑68 (MQ=255)
        gCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAgaga                                                             >  1:112620/1‑68 (MQ=255)
         ccAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCaa                                                                 >  1:1842467/1‑63 (MQ=255)
          cAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGc                                                        <  1:1948170/71‑1 (MQ=255)
               tCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGa                                                   >  1:1587109/1‑71 (MQ=255)
                 aGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGAtt                                                 <  1:248725/71‑1 (MQ=255)
                     cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGAt                                                  >  1:3462035/1‑66 (MQ=255)
                     cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGAt                                                  >  1:586812/1‑66 (MQ=255)
                     cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCg                                             >  1:3222429/1‑71 (MQ=255)
                       tAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCgg                                            >  1:4115956/1‑70 (MQ=255)
                       tAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCgg                                            >  1:925600/1‑70 (MQ=255)
                       tAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGc                                           <  1:1419209/71‑1 (MQ=255)
                        aaTTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATg                                                    <  1:1491643/61‑1 (MQ=255)
                           ttCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAagag                                                              >  1:925808/1‑48 (MQ=255)
                            tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAag                                                                <  1:4436032/45‑1 (MQ=255)
                            tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      >  1:2633998/1‑71 (MQ=255)
                            tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      >  1:3878398/1‑71 (MQ=255)
                             cGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      <  1:748440/70‑1 (MQ=255)
                               tGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAgaga                                                             <  1:3559461/45‑1 (MQ=255)
                               tGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGcc                                                       <  1:1219294/51‑1 (MQ=255)
                                 ccATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               >  1:4510693/1‑41 (MQ=255)
                                    tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCaa                                                                 <  1:4036013/36‑1 (MQ=255)
                                    tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAgaga                                                             <  1:3588784/40‑1 (MQ=255)
                                    tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAag                                                                <  1:402274/37‑1 (MQ=255)
                                      tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGAt                                                            >  1:342822/1‑39 (MQ=255)
                                      tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGcc                                                       >  1:4425868/1‑44 (MQ=255)
                                      tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGc                                                        <  1:4391517/43‑1 (MQ=255)
                                      tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCaa                                         >  1:2617282/1‑58 (MQ=255)
                                       cagcagTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTc                           >  1:1117113/1‑71 (MQ=255)
                                       cagcagTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTc                           >  1:4165139/1‑71 (MQ=255)
                                        agcagTGAGCTTACGTCCGGCATGTCCTTCAAGTGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCg                          <  1:1475849/71‑1 (MQ=255)
                                             tGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGc                                                        <  1:3335373/36‑1 (MQ=255)
                                             tGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGa                                                   <  1:3255343/41‑1 (MQ=255)
                                             tGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGc                         <  1:503967/67‑1 (MQ=255)
                                               aGCTTACGTCCGGCATGTCCCTCCAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGc                   <  1:989561/71‑1 (MQ=255)
                                               aGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTc                           >  1:766088/1‑63 (MQ=255)
                                                 cTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGc                         <  1:3579049/63‑1 (MQ=255)
                                                   tACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAggg               >  1:1559114/1‑71 (MQ=255)
                                                      gTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGgtgt            <  1:3226314/71‑1 (MQ=255)
                                                           gCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATATGCCGCTGTCGCATCTGCAGGGTGTTTTAc       >  1:4054492/1‑71 (MQ=255)
                                                                tCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTg                                   >  1:1207933/1‑38 (MQ=255)
                                                                tCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAc  >  1:3857216/1‑71 (MQ=255)
                                                                  cTTCAAGAGATACGCCATGATTAGCGGCAACATCTGc                                  >  1:4485464/1‑37 (MQ=255)
                                                                  cTTCAAGAGATACGCCATGATTAGCGGCAACATCTGc                                  >  1:2780408/1‑37 (MQ=255)
                                                                  cTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTa        >  1:3021472/1‑63 (MQ=255)
                                                                   ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCt                              <  1:3170293/40‑1 (MQ=255)
                                                                   ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAc  <  1:1371369/68‑1 (MQ=255)
                                                                   |                                                                   
ACTCCATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAC  >  W3110S.gb/1044837‑1044971

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: