Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,635,048 T→C 12.7% K322R (AAG→AGG)  der ← predicted GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,635,0480TC12.7% 137.6 / 12.5 55K322R (AAG→AGG) derpredicted GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base C (3/4);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

AAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCC  >  W3110S.gb/2634982‑2635110
                                                                  |                                                              
aaGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTAt                                                             <  1:2289395/70‑1 (MQ=255)
     cAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcacca                                                       >  1:759702/1‑71 (MQ=255)
     cAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcacca                                                       >  1:465543/1‑71 (MQ=255)
        aCGGAAGTCCAGCGTTTCTTTCTCCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcac                                                         <  1:2967077/66‑1 (MQ=255)
          ggAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTc                                                           >  1:3164296/1‑62 (MQ=255)
          ggAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAAt                                                   >  1:4590405/1‑70 (MQ=255)
          ggAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCACACTTATTcacca                                                       >  1:3103268/1‑66 (MQ=255)
              gTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                                          >  1:1647615/1‑59 (MQ=255)
               tCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTAtt                                                            <  1:2692969/56‑1 (MQ=255)
               tCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACa                                               <  1:3061988/69‑1 (MQ=255)
                 cAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcac                                                         <  1:4252692/57‑1 (MQ=255)
                 cAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             >  1:3633966/1‑69 (MQ=255)
                  aGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcac                                                         <  1:2339578/56‑1 (MQ=255)
                    cGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgc                                        >  1:3566482/1‑71 (MQ=255)
                    cGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACCTATTCACCACAATGACAAGTGAgc                                        <  1:3585961/71‑1 (MQ=255)
                        tctttcACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGcc                                     >  1:3356089/1‑70 (MQ=255)
                            tcACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACt                                                                >  1:1658046/1‑39 (MQ=255)
                            tcACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             >  1:3924586/1‑58 (MQ=255)
                            tcACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCAc                                  <  1:3949299/69‑1 (MQ=255)
                            tcACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTa                                >  1:2160467/1‑71 (MQ=255)
                              aCCTGCTCTTTCACTTCCTGACTTAGGCCATCCCACCTATTCACCACAATGACAAGTGAGCGCCCACTAt                               <  1:4217957/70‑1 (MQ=255)
                              aCCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTc                                                           >  1:3665232/1‑42 (MQ=255)
                               ccTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg                             <  1:1572500/71‑1 (MQ=255)
                                cTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACa                                               <  1:4599610/52‑1 (MQ=255)
                                   ctctTTCACTTCCTGACTCAGGCCATCCCACTTAtt                                                            >  1:4133758/1‑36 (MQ=255)
                                   ctctTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGaa                         >  1:639925/1‑71 (MQ=255)
                                   ctctTTCACTTCCTGACTCAGGCCATCCCACCTATTCACCACaa                                                    >  1:1380122/1‑44 (MQ=255)
                                   ctctTTCACTTCCTGACTCAGGCCATCCCACCTATTCACCACAATGACAAGTGAGCGCCCACTATTgag                           <  1:4079043/69‑1 (MQ=255)
                                    tctTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCAc                                  >  1:2571072/1‑61 (MQ=255)
                                    tctTTCACTTCCTGACTCAGGCCATCCCACCTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGaa                         >  1:4117581/1‑70 (MQ=255)
                                      tttCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAAt                                                   <  1:4303641/42‑1 (MQ=255)
                                      tttCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAAt                        <  1:282322/69‑1 (MQ=255)
                                        tCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             >  1:3240276/1‑46 (MQ=255)
                                        tCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGt                                            >  1:1712017/1‑47 (MQ=255)
                                        tCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATaa                      <  1:638556/69‑1 (MQ=255)
                                         cACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGc                   >  1:4141213/1‑71 (MQ=255)
                                          aCTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGACCGCCCACTATTGAGAATAAAGc                   <  1:2163872/70‑1 (MQ=255)
                                           cTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGAc                                                <  1:1015679/40‑1 (MQ=255)
                                           cTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGAc                                                <  1:79371/40‑1 (MQ=255)
                                           cTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGccc                 <  1:1652776/71‑1 (MQ=255)
                                             tCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             <  1:4547325/41‑1 (MQ=255)
                                             tCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCcag               >  1:3512894/1‑71 (MQ=255)
                                              ccTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             >  1:941837/1‑40 (MQ=255)
                                              ccTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                                             >  1:4508690/1‑40 (MQ=255)
                                               cTGACTCAGGCCATCCCACTTATTCACCACAATGACa                                               <  1:158725/37‑1 (MQ=255)
                                               cTGACTCAGGCCATCCCACTTATTCACCACAATGACa                                               <  1:1007645/37‑1 (MQ=255)
                                                   cTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCga         <  1:1315898/71‑1 (MQ=255)
                                                     cAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCgaga       <  1:3045325/71‑1 (MQ=255)
                                                      aGGCCATCCCACCTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCgag        <  1:2269928/69‑1 (MQ=255)
                                                         ccATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTgaga                          <  1:3045011/48‑1 (MQ=255)
                                                         ccATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTgaga                          <  1:587272/48‑1 (MQ=255)
                                                          cATCCCACCTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATcc  >  1:3207751/1‑71 (MQ=255)
                                                           aTCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAAt                        >  1:2870837/1‑48 (MQ=255)
                                                              ccACTTATTCACCACAATGACAAGTGAGCGCCCACt                                 <  1:159358/36‑1 (MQ=255)
                                                               cACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGaa                         <  1:4362080/43‑1 (MQ=255)
                                                                aCTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCgaga       <  1:1175275/60‑1 (MQ=255)
                                                                  ttATTCACCACAATGACAAGTGAGCGCCCACTATTGAGaa                         <  1:4291332/40‑1 (MQ=255)
                                                                  ttATTCACCACAATGACAAGTGAGCGCCCACTATAGAGAATAAAg                    >  1:2988420/1‑45 (MQ=255)
                                                                  |                                                              
AAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCC  >  W3110S.gb/2634982‑2635110

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: