Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,269,809 T→C 9.4% V174V (GTA→GTG murB ← UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,269,8090TC 9.4% 129.2 / 7.3 53V174V (GTA→GTGmurBUDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (18/30);  new base C (3/2);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCCAGACGAGTTAAGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCACTCTTT  >  minE/2269741‑2269877
                                                                    |                                                                    
tCCAGACGAGTTAAGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTAc                                                                    <  1:1162405/71‑1 (MQ=255)
  cAGACGAGTTAAGTCACCATACGTTAGTTCAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACaa                                                                  >  1:2120001/1‑71 (MQ=255)
  cAGACGAGTTAAGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACaa                                                                  >  1:316679/1‑71 (MQ=255)
          ttAAGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTAc                                                                    <  1:1913185/61‑1 (MQ=255)
           tAAGTCACCATACGTTAGTACAGGTTGCCCCTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGc                                                         >  1:543580/1‑71 (MQ=255)
            aaGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACaa                                                                  >  1:1127223/1‑61 (MQ=255)
            aaGTCACCATACGTTAGTACAGGTTGACACTCTTTTGGCAGACGCAGACCTACGGCTACaa                                                                  >  1:788799/1‑61 (MQ=255)
             aGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCgg                                                       <  1:1660962/71‑1 (MQ=255)
                  ccATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTcc                                                    <  1:1831906/69‑1 (MQ=255)
                   cATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTgg                                                 >  1:2152999/1‑71 (MQ=255)
                     tACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTAc                                                                    <  1:933434/50‑1 (MQ=255)
                     tACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTAc                                                                    <  1:768063/50‑1 (MQ=255)
                     tACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTAc                                                                    <  1:2172883/50‑1 (MQ=255)
                       cGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTAtt                                             <  1:901486/71‑1 (MQ=255)
                       cGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTAtt                                             <  1:35342/71‑1 (MQ=255)
                        gTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTg                                                  <  1:210099/65‑1 (MQ=255)
                          tAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGt                                                      >  1:1780342/1‑59 (MQ=255)
                          tAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCAt                                          >  1:602060/1‑71 (MQ=255)
                             tACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTgg                                                 <  1:1774395/61‑1 (MQ=255)
                                 ggTTGCCACTCTTTTGGCAGACTCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaa                                   <  1:57952/71‑1 (MQ=255)
                                  gTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaa                                   <  1:294916/70‑1 (MQ=255)
                                    tGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaaa                                 >  1:1615758/1‑70 (MQ=255)
                                     gCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAAt                                <  1:2152537/70‑1 (MQ=255)
                                     gCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAAt                                <  1:1515897/70‑1 (MQ=255)
                                        aCTCTTTTGGCAGACGCAGACCTACGGCCACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTg                            >  1:1797702/1‑71 (MQ=255)
                                        aCTCATTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTc                                                     >  1:1419284/1‑46 (MQ=255)
                                          tctTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGc                                                         <  1:158214/40‑1 (MQ=255)
                                          tctTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGt                                                      <  1:1692195/43‑1 (MQ=255)
                                           ctTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTcg                         >  1:1354044/1‑71 (MQ=255)
                                           ctTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTcg                         >  1:432959/1‑71 (MQ=255)
                                            ttttGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTgg                                                 <  1:1650736/46‑1 (MQ=255)
                                                ggCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGAt                     <  1:404931/70‑1 (MQ=255)
                                                 gcagacgcagacCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATa                    <  1:982463/70‑1 (MQ=255)
                                                 gcagacgcagacCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATa                    <  1:1664188/70‑1 (MQ=255)
                                                     acgcagacCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaaa                                 >  1:772030/1‑53 (MQ=255)
                                                     acgcagacCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               <  1:825173/71‑1 (MQ=255)
                                                      cgcagacCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               >  1:1365107/1‑70 (MQ=255)
                                                      cgcagacCTACGGCCACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               >  1:1563924/1‑70 (MQ=255)
                                                      cgcagacCTACGGCCACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               >  1:1163250/1‑70 (MQ=255)
                                                           acCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTcgc                        <  1:1230333/56‑1 (MQ=255)
                                                           acCTACGGCCACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               <  1:2031096/65‑1 (MQ=255)
                                                           acCTACGGCCACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCaa               <  1:923517/65‑1 (MQ=255)
                                                            ccTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACggg         >  1:1721842/1‑69 (MQ=255)
                                                             cTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGtt                                       <  1:1598509/39‑1 (MQ=255)
                                                             cTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaa                                   <  1:520919/43‑1 (MQ=255)
                                                             cTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCAc       >  1:1214168/1‑71 (MQ=255)
                                                             cTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCAc       >  1:1027255/1‑71 (MQ=255)
                                                               aCGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCActc     >  1:164471/1‑71 (MQ=255)
                                                                cGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaaa                                 <  1:1510667/42‑1 (MQ=255)
                                                                  gCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaa                                  <  1:482810/39‑1 (MQ=255)
                                                                  gCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCACTCttt  <  1:1469654/71‑1 (MQ=255)
                                                                   cTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaaa                                 <  1:854728/39‑1 (MQ=255)
                                                                   cTACAATAGCGAAGCGGTCCTGGTATTCATGTTTaaaaa                                 <  1:1208624/39‑1 (MQ=255)
                                                                    |                                                                    
TCCAGACGAGTTAAGTCACCATACGTTAGTACAGGTTGCCACTCTTTTGGCAGACGCAGACCTACGGCTACAATAGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGCCAAAACGGCACTCTTT  >  minE/2269741‑2269877

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: