Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 813,910 A→G 44.6% C60C (TGT→TGC prsA ← phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE813,9100AG44.6% 29.1 / 131.9 96C60C (TGT→TGCprsAphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base G (22/20);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCG  >  minE/813845‑813976
                                                                 |                                                                  
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGt                                                               >  1:86423/1‑71 (MQ=255)
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGt                                                               >  1:1544091/1‑71 (MQ=255)
 cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTg                                                              <  1:1780505/71‑1 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAg                                                                >  1:417040/1‑67 (MQ=255)
         cgacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGgat                                                      >  1:406555/1‑71 (MQ=255)
         cgacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGgat                                                      >  1:2527988/1‑71 (MQ=255)
                   cATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAata                                            <  1:2004841/71‑1 (MQ=255)
                   cATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAata                                            <  1:119781/71‑1 (MQ=255)
                   cATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAata                                            <  1:2272481/71‑1 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                           >  1:2936586/1‑71 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                           >  1:3777274/1‑71 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGAACAAGTGGACTGGATGATGAAAata                                            >  1:2616961/1‑70 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGcaca                                                                  <  1:4220634/47‑1 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAat                                             <  1:3156896/68‑1 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAat                                             <  1:2595125/68‑1 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAat                                             <  1:468838/68‑1 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAat                                             <  1:3416829/68‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAatat                                           >  1:1564033/1‑68 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAatat                                           >  1:393256/1‑68 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAatat                                           >  1:3976732/1‑68 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        >  1:402864/1‑71 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        <  1:399061/71‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcac                                        <  1:1124986/71‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcac                                        <  1:3968658/71‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcac                                        <  1:1757355/71‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcac                                        <  1:4089036/71‑1 (MQ=255)
                        gccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGa                                                            >  1:3868085/2‑50 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTATGGGCACAAGTGGa                                                            >  1:2104473/1‑50 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        <  1:2941083/70‑1 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGa                                                            >  1:2354567/1‑50 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGa                                                            >  1:410215/1‑50 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGa                                                            >  1:1327469/1‑50 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcacc                                       <  1:1076530/71‑1 (MQ=255)
                         ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAg                                                                <  1:733618/45‑1 (MQ=255)
                         ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAg                                                                <  1:1367325/45‑1 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGgatgat                                                   <  1:3974164/56‑1 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGgatgat                                                   <  1:2749691/56‑1 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGgatgat                                                   <  1:1312225/56‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcacc                                       <  1:3717022/67‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcacc                                       <  1:4118794/67‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcacc                                       <  1:3779264/67‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcacc                                       <  1:3814907/67‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTg                                                              <  1:3320635/44‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTg                                                              <  1:3039227/44‑1 (MQ=255)
                             aaCTACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACaact                                                               >  1:2589001/1‑40 (MQ=38)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACcgc                                  >  1:4173636/1‑71 (MQ=255)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACcgc                                  >  1:4053965/1‑71 (MQ=255)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACcgc                                  >  1:1317147/1‑71 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACcgcg                                 >  1:2045954/1‑71 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACcgcg                                 >  1:1466184/1‑71 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACcgcg                                 >  1:2998924/1‑71 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                             <  1:495557/58‑1 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                             <  1:2849373/58‑1 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                             <  1:2509283/58‑1 (MQ=255)
                                gACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACTAGTGGACTGGATGATGAAAATATCACCACCGCGTa                               <  1:511134/71‑1 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                              <  1:3475351/71‑1 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                              <  1:1909907/71‑1 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                              <  1:814788/71‑1 (MQ=255)
                                  cTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                              <  1:3480359/70‑1 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCTCAAGTGGACTGGATGATGAAAATATcac                                        >  1:3582207/1‑57 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        >  1:1841352/1‑57 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        >  1:3519183/1‑57 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATcac                                        >  1:4170487/1‑57 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                           >  1:3585028/1‑54 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                           >  1:3358509/1‑54 (MQ=255)
                                      ttCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACAtttt                         <  1:1190188/71‑1 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTc                        >  1:619970/1‑70 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                       >  1:1695248/1‑71 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                       >  1:933634/1‑71 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                       >  1:2382057/1‑71 (MQ=255)
                                          aTCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                     <  1:2190274/71‑1 (MQ=255)
                                          aTCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                     <  1:763115/71‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTa                    <  1:1254562/71‑1 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTaa                   >  1:2257584/1‑71 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTaa                   >  1:1475816/1‑71 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTaa                   >  1:3783108/1‑71 (MQ=255)
                                               gTTGTCGTTAGTAGGTGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                     <  1:3345220/66‑1 (MQ=255)
                                               gTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                     <  1:580893/66‑1 (MQ=255)
                                               gTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAtt                 <  1:3048996/70‑1 (MQ=255)
                                                ttGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaa                                                <  1:3020795/38‑1 (MQ=255)
                                                 tGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATca                                         >  1:2338917/1‑44 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTa                               >  1:1434028/1‑54 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTa                               >  1:3198212/1‑54 (MQ=255)
                                                  gTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaaaa                                              >  1:650343/1‑38 (MQ=255)
                                                  gTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaaaa                                              >  1:91978/1‑38 (MQ=255)
                                                     gTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGc          >  1:3429919/1‑71 (MQ=255)
                                                     gTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                     <  1:2925083/60‑1 (MQ=255)
                                                      ttAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                  <  1:1708762/62‑1 (MQ=255)
                                                      ttAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                  <  1:4218416/62‑1 (MQ=255)
                                                       tagtatGGGCGCAAGTTGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                  <  1:1555202/61‑1 (MQ=255)
                                                            ggggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACtt    <  1:476097/70‑1 (MQ=255)
                                                            ggggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACtt    <  1:579497/70‑1 (MQ=255)
                                                             gggCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                  <  1:2349778/55‑1 (MQ=255)
                                                             gggCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                  <  1:932136/55‑1 (MQ=255)
                                                              ggCACAAGTGGCCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTc   >  1:2516398/1‑69 (MQ=255)
                                                              ggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTc   >  1:3340508/1‑69 (MQ=255)
                                                              ggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCg  >  1:3779225/1‑70 (MQ=255)
                                                                 |                                                                  
GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCG  >  minE/813845‑813976

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: