Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,430,365 T→C 74.6% 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,3650TC74.6% 89.7 / 29.5 55I111M (ATA→ATGrfeUDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base C (24/17);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAG  >  minE/2430299‑2430433
                                                                  |                                                                    
gAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACcatcat                                                                    <  1:4688081/69‑1 (MQ=255)
        tctcCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGcg                                                           >  1:1591545/1‑69 (MQ=255)
          tcCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCaa                                                         >  1:3953817/1‑70 (MQ=255)
               ggAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGACAACAGCgg                                                   >  1:1876026/1‑71 (MQ=255)
                gAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGc                                                  <  1:3143576/71‑1 (MQ=255)
                 aGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGcc                                                 <  1:4639075/71‑1 (MQ=255)
                 aGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGcc                                                 <  1:3735624/71‑1 (MQ=255)
                     aaaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTa                                             >  1:3041887/1‑71 (MQ=255)
                     aaaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTa                                             <  1:996059/71‑1 (MQ=255)
                      aaGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAt                                            <  1:1363460/71‑1 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACa                                                       >  1:1682930/1‑59 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACa                                                       >  1:3401074/1‑59 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACa                                                       >  1:2955528/1‑59 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAg                                                      >  1:2353773/1‑60 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAg                                                      >  1:2851912/1‑60 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCg                                                    <  1:950974/62‑1 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGc                                                  >  1:3629612/1‑64 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtg                                           <  1:3881189/71‑1 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAtg                                           >  1:3320812/1‑71 (MQ=255)
                       aGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAt                                            >  1:3636176/1‑70 (MQ=255)
                           ataACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACa                                                                <  1:2839133/46‑1 (MQ=255)
                               cccAGGCTACTGAGATAAAGCTTGCCTAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCAc                                    >  1:4647286/1‑70 (MQ=255)
                               cccAGGCTACTGAGAAAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACg                                   >  1:3595707/1‑71 (MQ=255)
                                ccAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGc                                                     >  1:1234656/1‑52 (MQ=255)
                                  aGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAt                                >  1:3704513/1‑71 (MQ=255)
                                  aGGCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAt                                >  1:3746406/1‑71 (MQ=255)
                                    gCTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtt                               <  1:5627160/70‑1 (MQ=255)
                                     cTACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                             <  1:11425/71‑1 (MQ=255)
                                      tACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTa                                             >  1:3052358/1‑54 (MQ=255)
                                      tACTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                             >  1:3017120/1‑70 (MQ=255)
                                        cTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTAc                          <  1:2165112/71‑1 (MQ=255)
                                        cTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTg                                               <  1:2518382/50‑1 (MQ=255)
                                        cTGAGATAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTg                                               <  1:1269056/50‑1 (MQ=255)
                                         tGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACg                         >  1:4908/1‑71 (MQ=255)
                                          gagaTAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTAtg                                           >  1:726098/1‑52 (MQ=255)
                                           agaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGc                                                  >  1:152776/1‑44 (MQ=255)
                                            gaTAAAGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTg                                               <  1:1980333/46‑1 (MQ=255)
                                               aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCgg                                                   >  1:190055/1‑39 (MQ=255)
                                               aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACg                                   >  1:377346/1‑55 (MQ=255)
                                               aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGata                   >  1:5697217/1‑71 (MQ=255)
                                               aaaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGata                   >  1:4343376/1‑71 (MQ=255)
                                                aaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGAtttt                             <  1:44716/60‑1 (MQ=255)
                                                aaGCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGata                   <  1:5457600/70‑1 (MQ=255)
                                                  gCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAt                                            <  1:5621114/43‑1 (MQ=255)
                                                  gCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAt                                            <  1:1645848/43‑1 (MQ=255)
                                                  gCTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTAt                                            <  1:267178/43‑1 (MQ=255)
                                                   cTTGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaa               <  1:129569/71‑1 (MQ=255)
                                                     tGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTa                                             >  1:1524836/1‑39 (MQ=255)
                                                     tGCCGAACACCATCATAACAATGCCAACAGCGGCCTGTa                                             >  1:3865971/1‑39 (MQ=255)
                                                      gCCGAACACCATTATAACACTGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaaaa             <  1:5215364/70‑1 (MQ=255)
                                                      gCCGAACACCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaaa              <  1:5081049/69‑1 (MQ=255)
                                                           acacCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaa               >  1:2856093/1‑63 (MQ=255)
                                                           acacCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaa               >  1:3799802/1‑63 (MQ=255)
                                                           acacCATCATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCaa               >  1:1183979/1‑63 (MQ=255)
                                                               catcatAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATcc    >  1:243200/1‑70 (MQ=255)
                                                               cATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATc     >  1:1072832/1‑69 (MQ=255)
                                                                aTTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAg  >  1:5588071/1‑71 (MQ=255)
                                                                  |                                                                    
GAGCACCATCTCCCAGGAGCCAAAGATATAACCCAGGCTACTGAGATAAAGCTTGCCGAACACCATTATAACAATGCCAACAGCGGCCTGTATGGTGGCACGGATTTTTACGCTGATATCAAAACGGTCATCCAG  >  minE/2430299‑2430433

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: