Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,387,067 G→A 16.0% G264S (GGC→AGC)  elaC → binuclear zinc phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,387,0670GA16.0% 127.6 / 15.5 50G264S (GGC→AGC) elaCbinuclear zinc phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base A (4/4);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

GCCAAAGCCAATAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATG  >  W3110S.gb/2386998‑2387137
                                                                     |                                                                      
gccaaagccaaTAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCgg                                                                       >  1:352642/1‑71 (MQ=255)
     agccaaTAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGc                                                                  <  1:1702792/71‑1 (MQ=255)
     agccaaTAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGc                                                                  <  1:4761027/71‑1 (MQ=255)
       ccaaTAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCt                                                                 <  1:2229575/70‑1 (MQ=255)
        caaTAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTa                                                                <  1:4182293/70‑1 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCaa                                                                    >  1:741069/1‑61 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCac                                                           >  1:1228334/1‑69 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCAtt                                                          <  1:1727264/71‑1 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCAt                                                           >  1:728509/1‑70 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCAt                                                           <  1:866713/70‑1 (MQ=255)
             gTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCAt                                                           >  1:3075594/1‑70 (MQ=255)
              tCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCAtt                                                          <  1:1209126/70‑1 (MQ=255)
                   gCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTa                                                                <  1:688787/59‑1 (MQ=255)
                   gCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCATTACCCAc                                                    <  1:924732/71‑1 (MQ=255)
                      aTAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTc                                                 <  1:4704540/71‑1 (MQ=255)
                        aGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAg                                               >  1:4589815/1‑71 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAg                                               >  1:3697349/1‑70 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAg                                               <  1:1775825/70‑1 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAg                                               <  1:1650144/70‑1 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGc                                              <  1:2936122/71‑1 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGc                                              <  1:982472/71‑1 (MQ=255)
                         gCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGc                                              <  1:4164943/71‑1 (MQ=255)
                          ctctACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCATTACCCACGt                                                  <  1:1999191/66‑1 (MQ=255)
                          ctctACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCATTACCAACGt                                                  <  1:631919/66‑1 (MQ=255)
                             tACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGt                                                  <  1:445600/63‑1 (MQ=255)
                             tACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTcgc                                          >  1:1492057/1‑71 (MQ=255)
                             tACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTcg                                           >  1:2843379/1‑70 (MQ=255)
                               cacGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCa                                                            >  1:3211080/1‑51 (MQ=255)
                               cacGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCa                                                            >  1:1675801/1‑51 (MQ=255)
                                        cTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGAc                               <  1:87354/71‑1 (MQ=255)
                                         tGCGACACTAGCCCGTGAGGCTGGAGTCAGCAAGCTAATCATTACCCACGTCAGCTCGCGCTAtgatga                                >  1:615180/1‑69 (MQ=255)
                                          gCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCa                                                >  1:3646457/1‑52 (MQ=255)
                                          gCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCa                                                >  1:3320310/1‑52 (MQ=255)
                                          gCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACa                              >  1:3269496/1‑70 (MQ=255)
                                                 tAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCa                                                >  1:4692463/1‑45 (MQ=255)
                                                 tAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCATCTCGCGCTATGATGACAAAGGTTGt                      >  1:1999803/1‑71 (MQ=255)
                                                 tAGCCCGTGAGGCTGCAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGt                      >  1:4984069/1‑71 (MQ=255)
                                                  aGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGt                      >  1:985789/1‑70 (MQ=255)
                                                   gCCCGTGAGGCTGGAGTCAGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGt                      <  1:3594714/69‑1 (MQ=255)
                                                   gCCCCTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTTTGATGACAAAGGTTGt                      <  1:905805/69‑1 (MQ=255)
                                                      cGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCa                 >  1:3050123/1‑71 (MQ=255)
                                                      cGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCa                 >  1:533878/1‑71 (MQ=255)
                                                          aGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTg             <  1:3944056/71‑1 (MQ=255)
                                                          aGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTg             <  1:793739/71‑1 (MQ=255)
                                                           ggCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGt            <  1:4038960/71‑1 (MQ=255)
                                                              tGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTAGCGCTATGATGACAAAGGTTGTCAGCACCTGTTAc         >  1:1515058/1‑71 (MQ=255)
                                                               ggAGTCAGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACg        >  1:768656/1‑71 (MQ=255)
                                                                 aGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTg      <  1:4845124/71‑1 (MQ=255)
                                                                   tCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTc                     <  1:1409496/54‑1 (MQ=255)
                                                                   tCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTc                     <  1:2399709/54‑1 (MQ=255)
                                                                   tCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGAAAAAGGTTGTc                     <  1:2853190/54‑1 (MQ=255)
                                                                     ggCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAg                   >  1:3496321/1‑54 (MQ=255)
                                                                     ggCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATg  <  1:4620155/71‑1 (MQ=255)
                                                                     ggCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATg  >  1:4165802/1‑71 (MQ=255)
                                                                     ggCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATg  <  1:140042/71‑1 (MQ=255)
                                                                     ggCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATg  <  1:3148969/71‑1 (MQ=255)
                                                                     |                                                                      
GCCAAAGCCAATAGTCGCGGCCATAGCTCTACACGCCAGGCTGCGACACTAGCCCGTGAGGCTGGAGTCGGCAAGCTAATCATTACCCACGTCAGCTCGCGCTATGATGACAAAGGTTGTCAGCACCTGTTACGTGAATG  >  W3110S.gb/2386998‑2387137

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: