Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,941,030 T→G 5.3% T150T (ACT→ACG ygfZ → predicted folate‑dependent regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,941,0300TG 5.3% 188.2 / 5.5 75T150T (ACT→ACGygfZpredicted folate‑dependent regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (43/28);  new base G (2/2);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGC  >  minE/1940961‑1941097
                                                                     |                                                                   
cgcgcTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACtt                                                                    <  1:3500647/71‑1 (MQ=255)
 gcgcTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACttt                                                                   >  1:1047728/1‑71 (MQ=255)
  cgcTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg                                                                  >  1:1285810/1‑71 (MQ=255)
  cgcTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg                                                                  >  1:2609530/1‑71 (MQ=255)
   gcTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGc                                                                 <  1:703601/71‑1 (MQ=255)
    cTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGc                                                                 <  1:2045905/70‑1 (MQ=255)
       gCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg                                                             >  1:800316/1‑71 (MQ=255)
        cAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg                                                            >  1:3013991/1‑71 (MQ=255)
             ctctTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTg                                                        >  1:2094857/1‑70 (MQ=255)
              tctTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:1112120/1‑71 (MQ=255)
              tctTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:1400866/1‑71 (MQ=255)
                tttAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:4214258/1‑69 (MQ=255)
                  tAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACt                                                                     >  1:520176/1‑52 (MQ=255)
                  tAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACt                                                                     >  1:1139657/1‑52 (MQ=255)
                  tAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:1821679/1‑67 (MQ=255)
                  tAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAacac                                                   <  1:1519541/70‑1 (MQ=255)
                  tAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACAcc                                                  >  1:1538093/1‑71 (MQ=255)
                   aGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:2608461/1‑66 (MQ=255)
                   aGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACAcc                                                  <  1:542491/70‑1 (MQ=255)
                     cGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCgg                                               <  1:1699639/71‑1 (MQ=255)
                      gAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAG‑‑GCGACCACTTTGCTATGGTTTGAACACCCGGCAg                                            <  1:119068/71‑1 (MQ=255)
                         cTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTa                                                               >  1:2482274/1‑51 (MQ=255)
                         cTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:1403981/1‑60 (MQ=255)
                          tGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg                                                                  <  1:699981/47‑1 (MQ=255)
                          tGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg                                                            <  1:4358726/53‑1 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGGGACCACTTTGCTATGGtt                                                          >  1:1332520/1‑53 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACtt                                                                    >  1:2736844/1‑43 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACt                                                                     >  1:818846/1‑42 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGaa                                                      >  1:3144951/1‑57 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACg                                        >  1:1543446/1‑71 (MQ=255)
                            ccTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACg                                        >  1:438279/1‑71 (MQ=255)
                             cTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGt                                       <  1:1768656/71‑1 (MQ=255)
                               tCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAAc                                         >  1:742811/1‑67 (MQ=255)
                                cGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGtt                                      <  1:2596298/69‑1 (MQ=255)
                                cGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTc                                    >  1:4278848/1‑71 (MQ=255)
                                cGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTc                                    >  1:1736268/1‑71 (MQ=255)
                                 gaaagaaaAACAGGTAGTCAAAGAAGGCGCGACCACtt                                                                    >  1:2957587/1‑38 (MQ=255)
                                 gaaagaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTc                                    <  1:4091661/70‑1 (MQ=255)
                                  agagaaaAACAGGTAGTCAAAGAAGGCGCGACCACGTTGCTATGGTTTGAACACCCGGCAGAACGt                                       <  1:4559330/64‑1 (MQ=255)
                                  agagaaaAACAGGTAGTCAAAGAAGGCGCGACCACGTTGCTATGGTTTGAACACCCGGCAGAACGt                                       <  1:4220025/64‑1 (MQ=255)
                                  aaagaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAg                                            <  1:4110454/61‑1 (MQ=255)
                                  aaagaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCt                                  >  1:3379651/1‑71 (MQ=255)
                                   aagaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCt                                  >  1:4061097/1‑70 (MQ=255)
                                     gaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAAc                                         <  1:3641200/61‑1 (MQ=255)
                                     gaaaAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAAc                                         <  1:4345917/61‑1 (MQ=255)
                                            aGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGt                                                           <  1:2218723/36‑1 (MQ=255)
                                            aGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAAcc                        <  1:2366322/71‑1 (MQ=255)
                                                 gTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATg                    <  1:84072/70‑1 (MQ=255)
                                                  tCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAg                                            <  1:2585863/45‑1 (MQ=255)
                                                  tCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGtt                                      <  1:2398680/51‑1 (MQ=255)
                                                   cAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTg                                 <  1:635469/55‑1 (MQ=255)
                                                   cAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAg                 >  1:4444878/1‑71 (MQ=255)
                                                   cAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAg                 >  1:1300192/1‑71 (MQ=255)
                                                    aaaGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGa                                           >  1:3335001/1‑44 (MQ=255)
                                                    aaaGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAAc                         >  1:3234423/1‑62 (MQ=255)
                                                    aaaGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGc                <  1:1750528/71‑1 (MQ=255)
                                                    aaaGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGc                >  1:4205584/1‑71 (MQ=255)
                                                     aagaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGa                   <  1:2716001/67‑1 (MQ=255)
                                                      agaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTa              >  1:2615260/1‑71 (MQ=255)
                                                      agaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTa              >  1:2025364/1‑71 (MQ=255)
                                                       gaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCa                                             >  1:1167703/1‑39 (MQ=255)
                                                       gaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCa                                             >  1:4336908/1‑39 (MQ=255)
                                                       gaagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCa                                             >  1:4068146/1‑39 (MQ=255)
                                                        aagGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTaa                          >  1:2177649/1‑57 (MQ=255)
                                                         agGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTg           >  1:1389178/1‑71 (MQ=255)
                                                         agGCGCGACCACGTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTaa                          >  1:666586/1‑56 (MQ=255)
                                                         agGCGCGACCACGTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTg           >  1:4352275/1‑71 (MQ=255)
                                                          ggCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGc          >  1:1046923/1‑71 (MQ=255)
                                                           gcgcGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGt                            <  1:2712953/52‑1 (MQ=255)
                                                            cgcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCt         >  1:3458702/1‑70 (MQ=255)
                                                            cgcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCt         >  1:4695368/1‑70 (MQ=255)
                                                            cgcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTa        <  1:602725/71‑1 (MQ=255)
                                                            cgcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTa        <  1:984388/71‑1 (MQ=255)
                                                            cgcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTa        <  1:1556483/71‑1 (MQ=255)
                                                             gcgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAg                 >  1:484724/1‑61 (MQ=255)
                                                              cgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCt               <  1:3472217/62‑1 (MQ=255)
                                                              cgACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCt               <  1:1520268/62‑1 (MQ=255)
                                                                aCCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTAc             >  1:3822081/1‑62 (MQ=255)
                                                                 ccACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATg   >  1:1836600/1‑71 (MQ=255)
                                                                  cACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGc  >  1:4325760/1‑71 (MQ=255)
                                                                  cACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGc  >  1:3786701/1‑71 (MQ=255)
                                                                  cACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGc  >  1:1662558/1‑71 (MQ=255)
                                                                     tttGCTATGGTTTGAACACCCGGCAGAACGTTTCCTg                                 <  1:3537573/37‑1 (MQ=255)
                                                                     tttGCTATGGTTTGAACACCCGGCAGAACGTTTCCTg                                 <  1:3137477/37‑1 (MQ=255)
                                                                     tttGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCt         <  1:4677120/61‑1 (MQ=255)
                                                                     |                                                                   
CGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGC  >  minE/1940961‑1941097

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: