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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE813,9100AG26.0% 82.6 / 31.7 50C60C (TGT→TGCprsAphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (15/22);  new base G (9/4);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGC  >  minE/813845‑813977
                                                                 |                                                                   
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGt                                                                >  1:591546/1‑71 (MQ=255)
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGt                                                                >  1:1168802/1‑71 (MQ=255)
 cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAg                                                                 <  1:1026309/69‑1 (MQ=255)
 cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAg                                                                 <  1:1934981/69‑1 (MQ=255)
       cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTgg                                                         <  1:1342203/71‑1 (MQ=255)
             gCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGatgatg                                                   <  1:1822108/71‑1 (MQ=255)
               aGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaa                                                 <  1:364358/71‑1 (MQ=255)
                   cATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAata                                             <  1:1920082/71‑1 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAatat                                            >  1:320377/1‑71 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAata                                             >  1:1768174/1‑70 (MQ=255)
                     tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                              <  1:1219995/68‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaaaa                                               >  1:386017/1‑65 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaaaa                                               >  1:1956265/1‑65 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGaaaa                                               >  1:863940/1‑65 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACaa                                                                  >  1:1017742/1‑46 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACaa                                                                  >  1:1857123/1‑46 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcacc                                        <  1:1594597/71‑1 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATcac                                         >  1:1172662/1‑70 (MQ=255)
                         ccATAACGACTAGTTCCATCAGGTTGTCGTTAGTAGGGGCACAAg                                                                 >  1:19220/1‑45 (MQ=255)
                         ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                              <  1:378597/64‑1 (MQ=255)
                         ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAat                                              <  1:1903803/64‑1 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAata                                             >  1:925714/1‑63 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCaccacc                                     >  1:1754859/1‑71 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaaa                                               <  1:913475/60‑1 (MQ=255)
                            tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaaa                                               <  1:812520/60‑1 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGAc                                                            >  1:781935/1‑45 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACcgc                                   >  1:1430952/1‑70 (MQ=255)
                                gACTAATTCCATCAGGTTGTCGTTAGTAGGGGcaca                                                                   <  1:81161/36‑1 (MQ=255)
                                gACTAATTCCATCAGGTTGTCGTTAGTAGGGGcaca                                                                   <  1:47631/36‑1 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGAtgatga                                                  >  1:1035020/1‑52 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTGGACTGGATGATGAAAATATCACCACCGCGTa                                >  1:280883/1‑70 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                               <  1:1661340/71‑1 (MQ=255)
                                  cTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                               <  1:283703/70‑1 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaa                                                >  1:2105689/1‑52 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGtc                                                               <  1:1404165/35‑2 (MQ=38)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCGCAAGTg                                                               <  1:2157173/35‑1 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                            >  1:1688681/1‑54 (MQ=255)
                                         cATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAatat                                            <  1:2002589/50‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaa                                                 <  1:1259996/43‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTa                     <  1:595607/71‑1 (MQ=255)
                                           tCAGGGTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaa                                                 <  1:2184557/43‑1 (MQ=255)
                                             aGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                   >  1:932185/1‑71 (MQ=255)
                                                ttGTCGTTAGTAGGGGCACAAGTGTACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTg                <  1:751137/71‑1 (MQ=255)
                                                  gTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaaa                                               >  1:534470/1‑38 (MQ=255)
                                                  gTCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaa                                                >  1:1223523/1‑37 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGaaa                                                <  1:1930550/36‑1 (MQ=255)
                                                     gTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGc           >  1:146797/1‑71 (MQ=255)
                                                        agtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTg         <  1:350989/70‑1 (MQ=255)
                                                        agtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTg         <  1:1568319/70‑1 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGt               <  1:1975443/62‑1 (MQ=255)
                                                              ggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGc  <  1:92431/71‑1 (MQ=255)
                                                               gCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                        >  1:651889/1‑48 (MQ=255)
                                                                 |                                                                   
GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGC  >  minE/813845‑813977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: