Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,062,130 A→C 32.7% E163A (GAA→GCA)  argK → membrane ATPase/protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,062,1300AC32.7% 23.4 / 52.6 49E163A (GAA→GCA) argKmembrane ATPase/protein kinase
Reads supporting (aligned to +/- strand):  ref base A (17/16);  new base C (10/6);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAGGGCAT  >  W3110S.gb/3062063‑3062199
                                                                   |                                                                     
ttATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACATaa                                                                      <  1:377833/69‑1 (MQ=255)
ttATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGaa                                                                      <  1:2460665/69‑1 (MQ=255)
ttATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAg                                                                     <  1:912405/70‑1 (MQ=255)
ttATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAg                                                                     <  1:1177928/70‑1 (MQ=255)
ttATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAg                                                                     <  1:37691/70‑1 (MQ=255)
      gAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTAGCCCg                                                              <  1:2048450/71‑1 (MQ=255)
      gAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCg                                                              <  1:4269628/71‑1 (MQ=255)
      gAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCg                                                              <  1:1556613/71‑1 (MQ=255)
        agcagcGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGaa                                                                      <  1:3480494/61‑1 (MQ=255)
         gcatcGGGTTATGACGTAGTGACTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAt                                                           <  1:2819667/71‑1 (MQ=255)
         gcagcGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAt                                                           >  1:4227289/1‑71 (MQ=255)
              gggTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCAtggtg                                                       >  1:704144/1‑70 (MQ=255)
                gTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGc                                                                 <  1:3009733/58‑1 (MQ=255)
                gTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGAc                                                    <  1:3008465/71‑1 (MQ=255)
                    tGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGtt                                                >  1:4188920/1‑71 (MQ=255)
                     gACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGttt                                               >  1:2739332/1‑71 (MQ=255)
                       cGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGt                                                                    >  1:3545655/1‑48 (MQ=255)
                        gTAGTGATTGTCGAAACGGTTGGCGTCGGGCTGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTAt                                            <  1:2488859/71‑1 (MQ=255)
                         tAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAg                                                                     <  1:2993573/45‑1 (MQ=255)
                         tAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCg                                                              <  1:3876731/52‑1 (MQ=255)
                          aGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGt                                                                    >  1:130919/1‑45 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGaa                                                                      >  1:1988052/1‑42 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTc                                                                   >  1:1262335/1‑45 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTc                                                                   >  1:36492/1‑45 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCa                                                            >  1:3338153/1‑52 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCa                                                            >  1:3636095/1‑52 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTAtctc                                         >  1:1475108/1‑71 (MQ=255)
                           gTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTAt                                            >  1:2668877/1‑68 (MQ=255)
                            tGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAAcagcag                                                                     <  1:3915160/42‑1 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCAtgg                                                         >  1:4297305/1‑52 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCAtgg                                                         >  1:60674/1‑52 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCAtgg                                                         >  1:1994006/1‑52 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGtt                                      >  1:1270788/1‑71 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGtt                                      >  1:2611717/1‑71 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGt                                       <  1:1091782/70‑1 (MQ=255)
                              aTTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGt                                       <  1:1442840/70‑1 (MQ=255)
                                 gTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGcc                                                                >  1:3831647/1‑42 (MQ=255)
                                  tCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCa                                   <  1:2768926/70‑1 (MQ=255)
                                     aaaCGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAAtt                               <  1:540354/71‑1 (MQ=255)
                                        cGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGcc                            <  1:1057290/71‑1 (MQ=255)
                                          gTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCgg                          >  1:4588203/1‑71 (MQ=255)
                                                 tCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCaaa                                 >  1:2190195/1‑57 (MQ=255)
                                                 tCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCaa                                  >  1:2195476/1‑56 (MQ=255)
                                                      cAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCgatga              <  1:3330586/71‑1 (MQ=255)
                                                        gTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATc            <  1:2463058/71‑1 (MQ=255)
                                                         tCGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTggcggc                   <  1:1182223/63‑1 (MQ=255)
                                                          cGGAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTg          >  1:642564/1‑71 (MQ=255)
                                                           ggAAACAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCga                 >  1:1112508/1‑63 (MQ=255)
                                                            gAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGc         <  1:2570818/70‑1 (MQ=255)
                                                              aaCAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAgg      >  1:2623441/1‑71 (MQ=255)
                                                               aCAGCAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAggg     <  1:4296134/71‑1 (MQ=255)
                                                               aCAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTgg                       >  1:1787433/1‑53 (MQ=255)
                                                               aCAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAggg     >  1:2036848/1‑71 (MQ=255)
                                                                  gAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAGGGCAt  >  1:2884660/1‑71 (MQ=255)
                                                                   |                                                                     
TTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGACTGTTTTATCTCGTTGCAAATTGCCGGTGGCGGCGATGATCTGCAGGGCAT  >  W3110S.gb/3062063‑3062199

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: