Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,042,843 Δ1 bp 27.1% coding (1790/2181 nt) yccC ← cryptic autophosphorylating protein tyrosine kinase Etk

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,042,8430A.27.1% 138.5 / 68.0 70coding (1790/2181 nt)yccCcryptic autophosphorylating protein tyrosine kinase Etk
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base . (9/10);  total (33/37)
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.40e-01

AGGTGGCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGTTATGC  >  W3110S.gb/1042776‑1042908
                                                                   |                                                                 
aGGTCGCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATTCTGGATCACTTTGTTGAGCt                                                                >  1:3439818/1‑71 (MQ=255)
    ggCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCa                                                              >  1:966075/1‑69 (MQ=255)
    ggCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCa                                                              >  1:2716353/1‑69 (MQ=255)
    ggCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTAGCTCa                                                              >  1:286413/1‑68 (MQ=255)
       aCCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCttt                                                         <  1:1053059/71‑1 (MQ=255)
       aCCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCttt                                                         <  1:3920588/71‑1 (MQ=255)
       aCCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGAGCACTTTGTTGAGCTCATCttt                                                         <  1:1286081/71‑1 (MQ=255)
        ccTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTa                                                        <  1:981947/71‑1 (MQ=255)
        ccTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTa                                                        <  1:3096458/71‑1 (MQ=255)
         cTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAc                                                       <  1:3399350/71‑1 (MQ=255)
           gACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCt                                                     >  1:914491/1‑71 (MQ=255)
           gACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCt                                                     >  1:499502/1‑71 (MQ=255)
                 gAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGa                                               <  1:1947091/71‑1 (MQ=255)
                 gAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGa                                               <  1:2722688/71‑1 (MQ=255)
                 gAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGa                                               <  1:82973/71‑1 (MQ=255)
                   gTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGatat                                            <  1:3125190/71‑1 (MQ=255)
                   gTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGatat                                            <  1:3660945/71‑1 (MQ=255)
                   gTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGatat                                            <  1:1991312/71‑1 (MQ=255)
                    tAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGatat                                            <  1:2795972/70‑1 (MQ=255)
                      aTCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAcc                                                      >  1:1342106/1‑59 (MQ=255)
                      aTCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAcc                                                      >  1:4020292/1‑59 (MQ=255)
                        cacaTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCttt                                                         <  1:3412537/54‑1 (MQ=255)
                         acaTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCAt                                                             >  1:2114898/1‑48 (MQ=255)
                          caTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGAc                                      <  1:2380354/71‑1 (MQ=255)
                          caTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGAc                                      <  1:92927/71‑1 (MQ=255)
                           aTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGAc                                      >  1:2215518/1‑69 (MQ=255)
                           aTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGAc                                      >  1:2254242/1‑69 (MQ=255)
                           aTCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGAc                                      >  1:2642230/1‑69 (MQ=255)
                            tCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCttt                                                         <  1:513459/50‑1 (MQ=255)
                             cAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGata                                             >  1:2162331/1‑61 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAc                                                       >  1:2745596/1‑50 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAc                                                       >  1:1337732/1‑50 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAc                                                       >  1:369448/1‑50 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTAc                                                       >  1:3030298/1‑50 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATAtt                                           >  1:2365273/1‑62 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCATATATTCCGACAAGc                                  >  1:795978/1‑70 (MQ=255)
                              aaaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGc                                  >  1:1846449/1‑70 (MQ=255)
                               aaGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGcc                                 <  1:102329/70‑1 (MQ=255)
                                 gCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGc                                                                 >  1:3823651/1‑37 (MQ=255)
                                  cctcctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGc                                  <  1:1168230/67‑1 (MQ=255)
                                  cctcctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGc                                  <  1:3198013/67‑1 (MQ=255)
                                  cctcctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGc                                  <  1:3259338/67‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTTAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGtt                            <  1:1000520/69‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTTAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGtt                            <  1:1965601/69‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTc                           <  1:896657/70‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTc                           <  1:1215307/70‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCa                          >  1:2442445/1‑71 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCAt                         <  1:35050/71‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCAt                         <  1:3517142/71‑1 (MQ=255)
                                      ctTTGCCAAAATGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCAt                         <  1:1315371/71‑1 (MQ=255)
                                       tttGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGtt                            <  1:3371763/68‑1 (MQ=255)
                                         tGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATg                              >  1:2610306/1‑64 (MQ=255)
                                         tGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCa                          >  1:3619538/1‑68 (MQ=255)
                                         tGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCa                          >  1:2237462/1‑68 (MQ=255)
                                                 tGCTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTa              <  1:743574/71‑1 (MQ=255)
                                                   cTGGATCACTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTaaa            <  1:661472/71‑1 (MQ=255)
                                                         cACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACTAGCCATGTTc                           <  1:3991696/51‑1 (MQ=255)
                                                         cACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGt       >  1:1397478/1‑71 (MQ=255)
                                                         cACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGt       >  1:2667291/1‑71 (MQ=255)
                                                         cACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGt       >  1:1141130/1‑71 (MQ=255)
                                                         cACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGt       >  1:3312377/1‑71 (MQ=255)
                                                          aCTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGt       <  1:4033739/70‑1 (MQ=255)
                                                          aCTTTGTTGGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGt                             >  1:195171/1‑47 (MQ=255)
                                                            tttGTTGAGCTCATCTTTACCTTCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACa          <  1:1526664/65‑1 (MQ=255)
                                                            tttGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACa          <  1:409449/65‑1 (MQ=255)
                                                            tttGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACa          <  1:2811125/65‑1 (MQ=255)
                                                            tttGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACa          <  1:2190442/65‑1 (MQ=255)
                                                             ttgttgGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGTTATGc  >  1:1335311/1‑71 (MQ=255)
                                                             ttgttgGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGTTATGc  >  1:3355944/1‑71 (MQ=255)
                                                              tgttgAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGTTATGc  >  1:3167571/1‑71 (MQ=255)
                                                                   |                                                                 
AGGTGGCACCTGACCGCGAGTAATCACATCAAAGCCTCCTTTGCCAAAATGCTGGATCACTTTGTTGAGCTCATCTTTACCTGCCAGATATTCCGACAAGCCATGTTCATTACTCACGGTAAACAGGTTATGC  >  W3110S.gb/1042776‑1042908

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: