Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,052,400 T→C 12.8% Y33H (TAC→CAC)  ugpQ → glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,052,4000TC12.8% 185.0 / 17.3 70Y33H (TAC→CAC) ugpQglycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base T (39/22);  new base C (3/6);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GTCGCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATG  >  W3110S.gb/4052331‑4052466
                                                                     |                                                                  
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAAt                                                                    >  1:2192028/1‑70 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:2629598/1‑71 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:1823451/1‑71 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:2738527/1‑71 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:3344767/1‑71 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:621225/1‑71 (MQ=255)
gtcgCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATa                                                                   >  1:916267/1‑71 (MQ=255)
      cATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCa                                                             >  1:2727845/1‑71 (MQ=255)
          gTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCAt                                                            >  1:1244162/1‑68 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATc                                                   >  1:1578746/1‑69 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCg                                                  >  1:490154/1‑70 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCg                                                  >  1:649185/1‑70 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCg                                                  >  1:2693033/1‑70 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCg                                                  >  1:1987240/1‑70 (MQ=255)
                  ggTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCg                                                  >  1:3073358/1‑70 (MQ=255)
                        cTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAAttt                                             >  1:2034387/1‑69 (MQ=255)
                                        cccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGaaa                           <  1:2251991/71‑1 (MQ=255)
                                        cccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGaaa                           <  1:217098/71‑1 (MQ=255)
                                        cccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGaaa                           <  1:945281/71‑1 (MQ=255)
                                        cccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGaaa                           <  1:1809022/71‑1 (MQ=255)
                                        cccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGaaa                           <  1:55290/71‑1 (MQ=255)
                                         ccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAg                          <  1:3109061/71‑1 (MQ=255)
                                         ccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAg                          <  1:3341492/71‑1 (MQ=255)
                                         ccTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAg                          <  1:1332539/71‑1 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:2991618/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:3199877/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:3300084/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:571190/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:2657385/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:1167939/1‑71 (MQ=255)
                                            ggCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATg                       >  1:1524508/1‑71 (MQ=255)
                                                tCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:1057383/55‑1 (MQ=255)
                                                tCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:1392096/55‑1 (MQ=255)
                                                tCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:1801553/55‑1 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:2556222/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:2584746/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:519610/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:440160/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:1841069/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:2808964/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:1509299/1‑55 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:1108223/70‑1 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:2391665/70‑1 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:1736038/70‑1 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:2974640/70‑1 (MQ=255)
                                                  aaTCGACGTCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:2608981/70‑1 (MQ=255)
                                                      gACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCtt             >  1:2847682/1‑71 (MQ=255)
                                                      gACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCtt             >  1:777462/1‑71 (MQ=255)
                                                      gACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCtt             >  1:3205715/1‑71 (MQ=255)
                                                      gACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCtt             >  1:1352110/1‑71 (MQ=255)
                                                      gACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCtt             >  1:1916710/1‑71 (MQ=255)
                                                          tCGGGGCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCgag                  <  1:1148161/62‑1 (MQ=255)
                                                            ggggCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCt       >  1:1526429/1‑71 (MQ=255)
                                                            ggggCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCt       >  1:2834530/1‑71 (MQ=255)
                                                            ggggCAAAACACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCt       >  1:143808/1‑71 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACTCGAAGt                                   <  1:52831/38‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:2482691/71‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:739957/71‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:1357352/71‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:405306/71‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:3161581/71‑1 (MQ=255)
                                                                 aaaaTACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATg  <  1:2867531/71‑1 (MQ=255)
                                                                  aaaTACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:2390176/37‑1 (MQ=255)
                                                                  aaaTACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:1460102/37‑1 (MQ=255)
                                                                  aaaTACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:2819697/37‑1 (MQ=255)
                                                                  aaaTACGGTCATAAGATGATCGAATTTGACGCGAAGt                                   <  1:897823/37‑1 (MQ=255)
                                                                    aTACGGTCATAAGATGATCGTATTTGACGCGAAGTTa                                 >  1:2932463/1‑37 (MQ=255)
                                                                    aTACGGTCATAAGATGATCGATTTTGACGCGAAGTTa                                 >  1:2030483/1‑37 (MQ=255)
                                                                    aTACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:387905/1‑37 (MQ=255)
                                                                    aTACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:1253857/1‑37 (MQ=255)
                                                                    aTACGGTCATAAGATGATCGAATTTGACGCGAAGTTa                                 >  1:3298618/1‑37 (MQ=255)
                                                                     |                                                                  
GTCGCTCATCGTGGCGGCGGTAAGCTGGCCCCGGAAAACACCCTGGCGTCAATCGACGTCGGGGCAAAATACGGTCATAAGATGATCGAATTTGACGCGAAGTTATCGAAAGATGGCGAGATCTTCCTGCTCCATG  >  W3110S.gb/4052331‑4052466

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: