Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,970 G→A 5.9% G108S (GGC→AGC)  pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,9700GA 5.9% 278.1 / 5.0 84G108S (GGC→AGC) pyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (34/45);  new base A (3/2);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGG  >  W3110S.gb/3824901‑3825037
                                                                     |                                                                   
cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                                                                    <  1:4897471/71‑1 (MQ=255)
   gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                                                  <  1:1372096/70‑1 (MQ=255)
   gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAg                                                                     >  1:3027003/1‑67 (MQ=255)
         catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATc                                                            >  1:809945/1‑70 (MQ=255)
         catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATc                                                            >  1:3382597/1‑70 (MQ=255)
         catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCt                                                           <  1:4412349/71‑1 (MQ=255)
          atcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtg                                                          <  1:535606/71‑1 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc                                                                >  1:2629154/1‑62 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc                                                                >  1:1490569/1‑62 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc                                                                >  1:3800021/1‑62 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc                                                                >  1:3487678/1‑62 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                                                  >  1:4223994/1‑60 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                                                  >  1:708664/1‑60 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                                                  >  1:915576/1‑60 (MQ=255)
             atGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt                                                       >  1:336096/1‑71 (MQ=255)
                accttgacctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttggt                                                    <  1:1363076/71‑1 (MQ=255)
                accttgacctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttggt                                                    <  1:643700/71‑1 (MQ=255)
                accttgacctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg                                                     >  1:2566667/1‑70 (MQ=255)
                   ttgacctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                                        >  1:2768917/1‑64 (MQ=255)
                    tgacctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgcg                                                >  1:2585346/1‑71 (MQ=255)
                      acctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAGTCTGGTTGGTAGcgcg                                              >  1:3745522/1‑71 (MQ=255)
                      acctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                                              >  1:2486817/1‑71 (MQ=255)
                      acctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                                              >  1:1797508/1‑71 (MQ=255)
                       cctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGt                                             <  1:2749202/71‑1 (MQ=255)
                        ctGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGtt                                            <  1:2688517/71‑1 (MQ=255)
                           ccGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACa                                         <  1:2271168/71‑1 (MQ=255)
                              tACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                                       <  1:561261/70‑1 (MQ=255)
                                cTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttggt                                                    <  1:4694505/55‑1 (MQ=255)
                                cTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAc                                    >  1:1993466/1‑71 (MQ=255)
                                cTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAc                                    >  1:1765588/1‑71 (MQ=255)
                                 tGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                                              <  1:4546431/60‑1 (MQ=255)
                                 tGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                                              <  1:2675185/60‑1 (MQ=255)
                                   cTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcgcg                                 >  1:2409512/1‑71 (MQ=255)
                                   cTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcgcg                                 >  1:4343231/1‑71 (MQ=255)
                                      tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:2289671/69‑1 (MQ=255)
                                      tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:3690780/69‑1 (MQ=255)
                                      tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:4412154/69‑1 (MQ=255)
                                      tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTaa                              >  1:774416/1‑71 (MQ=255)
                                        aCCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATg                            <  1:943075/71‑1 (MQ=255)
                                          cGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcgcg                                 <  1:4016340/64‑1 (MQ=255)
                                          cGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCt                          <  1:3618391/71‑1 (MQ=255)
                                          cGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCt                          <  1:3284066/71‑1 (MQ=255)
                                           gCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTg                         >  1:2637173/1‑71 (MQ=255)
                                              aaGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTa                      <  1:3698224/71‑1 (MQ=255)
                                                gagGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTaa                              >  1:2655652/1‑61 (MQ=255)
                                                gagGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTaa                              >  1:1996479/1‑61 (MQ=255)
                                                gagGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTaa                              >  1:2126021/1‑61 (MQ=255)
                                                 agGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:4385545/1‑70 (MQ=255)
                                                 agGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:1731366/1‑70 (MQ=255)
                                                 agGCAAAAGACCACGGTGAAAGCGGCAATCTGGTTGGTAGcgcg                                              <  1:5142899/44‑1 (MQ=255)
                                                      aaaGACCACGGTGAAGGCGGCAATCTGGTTGGTAgcgc                                               <  1:265350/38‑1 (MQ=255)
                                                      aaaGACCACGGTGAAGGCGGCAATCTGGTTGGTAgcgc                                               <  1:859257/38‑1 (MQ=255)
                                                      aaaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                                              <  1:4233199/39‑1 (MQ=255)
                                                      aaaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgt              <  1:1552793/71‑1 (MQ=255)
                                                      aaaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgt              <  1:2177622/71‑1 (MQ=255)
                                                       aaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACa                                         <  1:2681195/43‑1 (MQ=255)
                                                       aaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAgatg                  >  1:1769497/1‑66 (MQ=255)
                                                       aaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgtg             <  1:1539715/71‑1 (MQ=255)
                                                       aaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgtg             <  1:1398871/71‑1 (MQ=255)
                                                       aaGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgt              >  1:3557480/1‑70 (MQ=255)
                                                        aGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAgatga                 >  1:2305477/1‑66 (MQ=255)
                                                        aGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAgatga                 >  1:5034266/1‑66 (MQ=255)
                                                         gACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGAt           <  1:2772328/71‑1 (MQ=255)
                                                          aCCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAgatga                 <  1:1309175/64‑1 (MQ=255)
                                                          aCCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAgatga                 <  1:1834799/64‑1 (MQ=255)
                                                          aCCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATc          <  1:4352945/71‑1 (MQ=255)
                                                          aCCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATc          <  1:3069458/71‑1 (MQ=255)
                                                          aCCACGGTGAAAGCGGCTATCTGGTTGGTAGcgcg                                              <  1:1874566/35‑1 (MQ=37)
                                                           ccACGGTGAAGGCGGCAATCTGGTTGTTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGAt           <  1:4765314/69‑1 (MQ=255)
                                                           ccACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGAt           <  1:3199891/69‑1 (MQ=255)
                                                           ccACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGAt           <  1:401400/69‑1 (MQ=255)
                                                           ccACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATc          <  1:1648333/70‑1 (MQ=255)
                                                           ccACGGTGAAGACGGCAATCTGGTTGTTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGAtgt              <  1:112293/66‑1 (MQ=255)
                                                            cACGGTGAAGGCGGCAATCTTGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:2309250/47‑1 (MQ=255)
                                                            cACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:2598608/47‑1 (MQ=255)
                                                            cACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:3333692/47‑1 (MQ=255)
                                                            cACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTa                      <  1:151528/57‑1 (MQ=255)
                                                            cACGGTGAAAGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:3334072/1‑59 (MQ=255)
                                                            cACGGTGAAAGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:277384/1‑59 (MQ=255)
                                                            cACGGTGAAAGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:3337400/1‑59 (MQ=255)
                                                              cGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:755383/1‑57 (MQ=255)
                                                              cGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:763610/1‑57 (MQ=255)
                                                              cGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAga                    >  1:1346927/1‑57 (MQ=255)
                                                                gTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                                       <  1:931702/36‑1 (MQ=255)
                                                                gTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                                       <  1:922670/36‑1 (MQ=255)
                                                                gTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                                       <  1:4445191/36‑1 (MQ=255)
                                                                gTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCAccgcc    >  1:5024864/1‑71 (MQ=255)
                                                                gTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCAccgcc    >  1:2832139/1‑71 (MQ=255)
                                                                  gAAGGCGGCAATCTGGTTGGTAGCGCGTTGCAAGGACGCGTAATGCTGGTAGATGATGTGATCAccgccg   >  1:3808536/1‑70 (MQ=255)
                                                                  gAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTa                               >  1:1180460/1‑42 (MQ=255)
                                                                  gAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCAccgccg   >  1:1877665/1‑70 (MQ=255)
                                                                  gAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCAccgccg   >  1:4436046/1‑70 (MQ=255)
                                                                   aaGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt                                <  1:28413/40‑1 (MQ=255)
                                                                   aaGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCgg  >  1:83288/1‑70 (MQ=255)
                                                                     |                                                                   
CTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGG  >  W3110S.gb/3824901‑3825037

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: