Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,044,904 T→C 38.7% 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,9040TC38.7% 53.6 / 85.0 75E65E (GAA→GAGetpphosphotyrosine‑protein phosphatase
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base C (7/22);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 4.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAC  >  W3110S.gb/1044841‑1044971
                                                               |                                                                   
cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               <  1:2531333/70‑1 (MQ=255)
cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCaa                                                                 >  1:2957315/1‑68 (MQ=255)
cATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCaa                                                                 >  1:3656910/1‑68 (MQ=255)
 aTTGGCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAgaga                                                             <  1:1605462/71‑1 (MQ=255)
 aTTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAgaga                                                             <  1:501920/71‑1 (MQ=255)
    gCCATAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATAc                                                          >  1:3191395/1‑71 (MQ=255)
    gCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATAc                                                          >  1:3518402/1‑71 (MQ=255)
         aaTCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCAt                                                     <  1:3610116/71‑1 (MQ=255)
         aaTCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCAt                                                     <  1:1047265/71‑1 (MQ=255)
          aTCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCa                                                                  <  1:188785/57‑1 (MQ=255)
          aTCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATg                                                    <  1:2571319/71‑1 (MQ=255)
           tCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATAc                                                          <  1:430585/64‑1 (MQ=255)
           tCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATAc                                                          <  1:2825882/64‑1 (MQ=255)
           tCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGa                                                   >  1:3974140/1‑71 (MQ=255)
           tCAGATCGTAATTTCGTGCCATATCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGa                                                   >  1:3707044/1‑71 (MQ=255)
               aTCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAaga                                                               <  1:3887621/55‑1 (MQ=255)
               aTCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAaga                                                               <  1:3418905/55‑1 (MQ=255)
               aTCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAaga                                                               <  1:3599979/55‑1 (MQ=255)
               aTCGAAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAaga                                                               <  1:3944976/55‑1 (MQ=39)
                tCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACTCCa                                                      >  1:3607693/1‑63 (MQ=255)
                 cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTa                                                >  1:2686835/1‑68 (MQ=255)
                 cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTa                                                >  1:314671/1‑68 (MQ=255)
                 cGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTa                                                >  1:1805996/1‑68 (MQ=255)
                     aTTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAaga                                                               >  1:798748/1‑49 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATa                                                           >  1:3902506/1‑50 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATa                                                           >  1:727266/1‑50 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATa                                                           >  1:3897092/1‑50 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACa                                       >  1:3226361/1‑70 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      >  1:2061257/1‑71 (MQ=255)
                        tCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      >  1:741335/1‑71 (MQ=255)
                         cGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACAt                                      <  1:178597/70‑1 (MQ=255)
                               aTCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAAc                                        <  1:1543089/62‑1 (MQ=255)
                               aTCTCAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAAc                                        <  1:379815/62‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGc                                                        <  1:1719948/45‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGc                                                        <  1:3766479/45‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCt                                    <  1:1042359/65‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCt                                    <  1:1143785/65‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGAtt                                                 <  1:891604/52‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGAtt                                                 <  1:1408417/52‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGAtt                                                 <  1:2077440/52‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGAtt                                                 <  1:162150/52‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAAc                                        <  1:3304665/61‑1 (MQ=255)
                                tctcAGCAGTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAAc                                        <  1:826524/61‑1 (MQ=255)
                                  tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGc                                  >  1:1016754/1‑65 (MQ=255)
                                  tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGc                                  >  1:273944/1‑65 (MQ=255)
                                  tcAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGt                            >  1:3949883/1‑71 (MQ=255)
                                      cagTGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCa                        <  1:3489807/71‑1 (MQ=255)
                                         tGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACa                                       <  1:4021385/53‑1 (MQ=255)
                                         tGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGc                                           <  1:4002799/49‑1 (MQ=255)
                                         tGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGc                                           <  1:2512493/49‑1 (MQ=255)
                                         tGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACa                                       <  1:35801/53‑1 (MQ=255)
                                         tGAGCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACa                                       <  1:866152/53‑1 (MQ=255)
                                          gAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGa                                                   >  1:933061/1‑40 (MQ=255)
                                           aGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTg                    <  1:3419626/70‑1 (MQ=255)
                                           aGCTTAAGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTg                    >  1:4161/1‑70 (MQ=255)
                                            gCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTg                    >  1:2445883/1‑69 (MQ=255)
                                            gCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTg                    >  1:909420/1‑69 (MQ=255)
                                            gCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCa                  <  1:1180059/71‑1 (MQ=255)
                                            gCTTACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCa                  <  1:1142206/71‑1 (MQ=255)
                                               tACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTa                                                >  1:2943701/1‑38 (MQ=255)
                                               tACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAggg               >  1:349281/1‑71 (MQ=255)
                                               tACGTCCGGCATGTCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAggg               >  1:531426/1‑71 (MQ=255)
                                                aCGTCCGGCATGTCCTTCAAGAGATACGCCATGATTA‑‑GGCAACATCTGc                                  >  1:982028/1‑49 (MQ=38)
                                                 cGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGgtg             <  1:1364099/71‑1 (MQ=255)
                                                     cGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGttt          >  1:2431066/1‑70 (MQ=255)
                                                     cGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGttt          >  1:2204149/1‑70 (MQ=255)
                                                       gCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTAc       <  1:334852/71‑1 (MQ=255)
                                                       gCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTAc       <  1:3273600/71‑1 (MQ=255)
                                                       gCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTAc       <  1:1880531/71‑1 (MQ=255)
                                                         aTGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCa     <  1:3412099/71‑1 (MQ=255)
                                                         aTGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCa     <  1:1280107/71‑1 (MQ=255)
                                                            tCCTTCAAGAGATACGCCATGATTAGCGGCAACAt                                      >  1:818709/1‑35 (MQ=255)
                                                            tCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGttt          >  1:437830/1‑63 (MQ=255)
                                                            tCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGttt          >  1:3399888/1‑63 (MQ=255)
                                                            tCCCTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAc  >  1:2983706/1‑71 (MQ=255)
                                                               |                                                                   
CATTGCCAGAATCAGATCGTAATTTCGTGCCATCTCAGCAGTGAGCTTACGTCCGGCATGTCCTTCAAGAGATACGCCATGATTAGCGGCAACATCTGCCGCTGTCGCATCTGCAGGGTGTTTTACCAGAC  >  W3110S.gb/1044841‑1044971

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: