Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,430,365 T→C 45.2% I111M (ATA→ATG rfe ← UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,430,3650TC45.2% 11.9 / 79.3 62I111M (ATA→ATGrfeUDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base T (20/14);  new base C (13/15);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CCGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAGC  >  minE/2430297‑2430434
                                                                    |                                                                     
ccGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATta                                                                      >  1:530950/1‑70 (MQ=255)
 cGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     >  1:3246988/1‑70 (MQ=255)
 cGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     <  1:2649503/70‑1 (MQ=255)
 cGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     <  1:39854/70‑1 (MQ=255)
 cGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     <  1:449073/70‑1 (MQ=255)
  gaacACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     <  1:2056754/66‑1 (MQ=255)
    gCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACa                                                                 <  1:2931554/71‑1 (MQ=255)
      aCCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAAt                                                               <  1:616394/71‑1 (MQ=255)
        cATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtat                                                                     <  1:2497100/63‑1 (MQ=255)
          tctcCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGcc                                                            >  1:411343/1‑70 (MQ=255)
              ccAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATtata                                                                    <  1:2997940/58‑1 (MQ=255)
              ccAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCaa                                                          <  1:2055069/68‑1 (MQ=255)
                aGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCg                                                     >  1:3022482/1‑71 (MQ=255)
                aGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCg                                                     >  1:3419524/1‑71 (MQ=255)
                aGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCg                                                     <  1:3702487/71‑1 (MQ=255)
                aGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCg                                                     <  1:2772344/71‑1 (MQ=255)
                 ggAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCg                                                     <  1:1805364/70‑1 (MQ=255)
                     ccAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGcc                                                            >  1:3229912/1‑59 (MQ=255)
                       aaaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTa                                              <  1:1004658/71‑1 (MQ=255)
                       aaaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTa                                              <  1:1051562/71‑1 (MQ=255)
                       aaaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTa                                              >  1:3702551/1‑71 (MQ=255)
                        aaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATg                                                              <  1:1520545/54‑1 (MQ=255)
                         aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAtg                                            >  1:1628340/1‑71 (MQ=255)
                         aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAtg                                            >  1:1111434/1‑71 (MQ=255)
                         aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATg                                                              <  1:1076174/53‑1 (MQ=255)
                          gATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACa                                                                 <  1:1159066/49‑1 (MQ=255)
                          gATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGcc                                                            >  1:2222422/1‑54 (MQ=255)
                          gATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAtgg                                           >  1:328943/1‑71 (MQ=255)
                             ataACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTa                                              <  1:1177442/65‑1 (MQ=255)
                              taACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtggtgg                                        <  1:2065153/70‑1 (MQ=255)
                                 cccAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACg                                    >  1:740747/1‑71 (MQ=255)
                                  ccAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACa                                                        <  1:2460143/50‑1 (MQ=255)
                                  ccAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACa                                                        <  1:441269/50‑1 (MQ=255)
                                     ggCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtt                                >  1:3351822/1‑71 (MQ=255)
                                      gCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAttt                               <  1:666280/71‑1 (MQ=255)
                                      gCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtggt                                          <  1:2526485/60‑1 (MQ=255)
                                       cTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTa                                              >  1:2263930/1‑55 (MQ=255)
                                       cTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                              >  1:3887394/1‑71 (MQ=255)
                                        tACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGc                                                      >  1:847405/1‑46 (MQ=255)
                                        tACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                              >  1:3194303/1‑70 (MQ=255)
                                        tACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                              >  1:2301753/1‑70 (MQ=255)
                                        tACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                              >  1:1456874/1‑70 (MQ=255)
                                          cTGAGATAAAGCTTGCCGAACACCATCATAACAATGcc                                                            <  1:698754/38‑1 (MQ=255)
                                           tGAGATAAAGCTTGCCGAACACCATCATAACAATGcc                                                            <  1:247316/37‑1 (MQ=255)
                                            gagaTAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGt                                               >  1:1031309/1‑49 (MQ=255)
                                            gagaTAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGt                                               >  1:3202306/1‑49 (MQ=255)
                                            gagaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAg                                                       >  1:1998880/1‑41 (MQ=255)
                                             agaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGTATTTTTACGCt                        <  1:3865893/71‑1 (MQ=255)
                                              gaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCg                                                     <  1:2891928/41‑1 (MQ=255)
                                              gaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCg                                                     <  1:2075274/41‑1 (MQ=255)
                                                 aaaGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGata                    >  1:809726/1‑71 (MQ=255)
                                                 aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtggtgg                                        >  1:877963/1‑51 (MQ=255)
                                                 aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtgg                                           >  1:2581099/1‑48 (MQ=255)
                                                 aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtt                                <  1:1721548/59‑1 (MQ=255)
                                                   aGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGatat                   >  1:3182296/1‑70 (MQ=255)
                                                      ttGCCGAACACCATTATAACAATGCCAACAGCGGCCt                                                 >  1:2094376/1‑37 (MQ=255)
                                                      ttGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaa                <  1:872748/70‑1 (MQ=255)
                                                        gCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAtggtg                                         >  1:1782349/1‑43 (MQ=255)
                                                         ccGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACg            <  1:2469688/71‑1 (MQ=255)
                                                          cGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGata                    >  1:3401858/1‑62 (MQ=255)
                                                            aaCACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAttt                               >  1:3351420/1‑49 (MQ=255)
                                                            aaCACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAttt                               >  1:2387357/1‑49 (MQ=255)
                                                             acacCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCa        <  1:202076/71‑1 (MQ=255)
                                                             acacCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCa        <  1:1690225/71‑1 (MQ=255)
                                                             acacCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCa        <  1:700376/71‑1 (MQ=255)
                                                              cacCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCAt       >  1:2826386/1‑71 (MQ=255)
                                                                 cATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAt                                 >  1:1298250/1‑42 (MQ=255)
                                                                 cATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCa    <  1:1036553/71‑1 (MQ=255)
                                                                  aTTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTc         >  1:3839558/1‑65 (MQ=255)
                                                                  aTTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTc         >  1:981462/1‑65 (MQ=255)
                                                                   ttATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAgc  >  1:608175/1‑71 (MQ=255)
                                                                   ttATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAgc  >  1:559657/1‑71 (MQ=255)
                                                                   ttATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAgc  >  1:975040/1‑71 (MQ=255)
                                                                    tataACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTc         >  1:665075/1‑63 (MQ=255)
                                                                    |                                                                     
CCGAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAGC  >  minE/2430297‑2430434

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: