Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,789 G→A 6.4% N101N (AAC→AAT rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,7890GA 6.4% 294.6 / 6.6 94N101N (AAC→AATrffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base A (4/2);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATTTTTTATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAT  >  minE/2425719‑2425857
                                                                      |                                                                    
aTTTTTTATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATg                                                                      >  1:1807915/1‑71 (MQ=255)
  tttttATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGtt                                                                    >  1:957597/1‑71 (MQ=255)
     ttATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGt                                                                 >  1:349846/1‑71 (MQ=255)
     ttATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGt                                                                 >  1:1328281/1‑71 (MQ=255)
            ccGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTc                                                              <  1:1628859/67‑1 (MQ=255)
            ccGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAAt                                                           <  1:1159179/70‑1 (MQ=255)
             cGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCa                                                             >  1:248737/1‑67 (MQ=255)
             cGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaa                                                         <  1:3219417/71‑1 (MQ=255)
             cGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaa                                                         <  1:134917/71‑1 (MQ=255)
              gTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaaa                                                        >  1:357337/1‑71 (MQ=255)
               tCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaa                                                         <  1:2823633/69‑1 (MQ=255)
               tCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAAc                                                       <  1:104806/71‑1 (MQ=255)
               tCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAAc                                                       <  1:2874746/71‑1 (MQ=255)
                cAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACg                                                      >  1:2291183/1‑71 (MQ=255)
                 aGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGc                                                     <  1:1669872/71‑1 (MQ=255)
                   cgcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTg                                                   <  1:1253654/71‑1 (MQ=255)
                    gcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTgg                                                                  <  1:1247140/55‑1 (MQ=255)
                    gcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaaa                                                        >  1:2233087/1‑65 (MQ=255)
                    gcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTg                                                   >  1:3280083/1‑70 (MQ=255)
                    gcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGc                                                  >  1:1393256/1‑71 (MQ=255)
                    gcATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGc                                                  >  1:381667/1‑71 (MQ=255)
                     cATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGcc                                                 >  1:261162/1‑71 (MQ=255)
                     cATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGcc                                                 >  1:350627/1‑71 (MQ=255)
                      aTTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCg                                                >  1:2524158/1‑71 (MQ=255)
                      aTTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCg                                                >  1:2704211/1‑71 (MQ=255)
                      aTTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCg                                                <  1:2777507/71‑1 (MQ=255)
                       ttCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCgg                                               >  1:3226780/1‑71 (MQ=255)
                        tCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGt                                                                 >  1:2952461/1‑52 (MQ=255)
                        tCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGt                                                                 >  1:294005/1‑52 (MQ=255)
                        tCCAGTAAGCCCGCGCCGATTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCgg                                               >  1:3884299/1‑70 (MQ=255)
                         ccAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGc                                             >  1:1891702/1‑71 (MQ=255)
                          cAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGc                                                  >  1:1035074/1‑65 (MQ=255)
                          cAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGcc                                            <  1:105697/71‑1 (MQ=255)
                            gTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCg                                           <  1:3314546/70‑1 (MQ=255)
                             tAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaaa                                                        <  1:1546688/56‑1 (MQ=255)
                              aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                         <  1:1631177/70‑1 (MQ=255)
                               aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                       >  1:687279/1‑71 (MQ=255)
                                gCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                       <  1:2441393/70‑1 (MQ=255)
                                 cccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                       >  1:529760/1‑69 (MQ=255)
                                  ccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTg                                                                   >  1:2413329/1‑40 (MQ=255)
                                  ccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                    <  1:1885767/71‑1 (MQ=255)
                                   cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATa                                     >  1:115415/1‑69 (MQ=255)
                                   cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                    >  1:2762674/1‑70 (MQ=255)
                                    gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                  >  1:264388/1‑71 (MQ=255)
                                      gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                <  1:699644/71‑1 (MQ=255)
                                      gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                <  1:1797915/71‑1 (MQ=255)
                                      gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTCCCGGGCCGTCAATAGAACg                                >  1:1662591/1‑71 (MQ=255)
                                       ccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                <  1:1733003/70‑1 (MQ=255)
                                        ctcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCaa                                                            <  1:3052140/39‑1 (MQ=255)
                                         gcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                             <  1:3654842/71‑1 (MQ=255)
                                            tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaa                                                         >  1:2135611/1‑40 (MQ=255)
                                            tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATaa                                                         >  1:1793314/1‑40 (MQ=255)
                                            tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAAt                                      <  1:1894911/59‑1 (MQ=255)
                                            tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCaa                           >  1:834530/1‑70 (MQ=255)
                                             cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                         <  1:1397775/55‑1 (MQ=255)
                                             cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                         >  1:1823733/1‑71 (MQ=255)
                                             cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                         >  1:3807353/1‑71 (MQ=255)
                                             cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                         >  1:1956674/1‑71 (MQ=255)
                                              aaGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                         <  1:2355245/70‑1 (MQ=255)
                                              aaGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                        <  1:2614328/71‑1 (MQ=255)
                                               aGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                        <  1:594339/70‑1 (MQ=255)
                                                  aaTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                        <  1:3639418/67‑1 (MQ=255)
                                                  aaTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                        <  1:1425877/67‑1 (MQ=255)
                                                   aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                          <  1:2374994/48‑1 (MQ=255)
                                                   aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCtt                   >  1:1387052/1‑71 (MQ=255)
                                                    tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                  <  1:2269879/71‑1 (MQ=255)
                                                    tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                  <  1:1103193/71‑1 (MQ=255)
                                                     gTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTc                 >  1:2844453/1‑71 (MQ=255)
                                                      tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTc                 <  1:2995065/70‑1 (MQ=255)
                                                      cataAGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCttt                  >  1:3209261/2‑69 (MQ=255)
                                                      cataAGTACCAACAATATTGGTTTCAATAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTcag               <  1:861050/70‑3 (MQ=255)
                                                        taAGTCCCGACAATGTTGGTTTCAATAAACGCTGcc                                                 >  1:1344396/1‑36 (MQ=255)
                                                        taAGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGcctg                                          >  1:1241989/1‑41 (MQ=16)
                                                         aaGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg               <  1:762616/69‑1 (MQ=255)
                                                         aaGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGc              >  1:1957527/1‑70 (MQ=255)
                                                          aGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTg            <  1:172685/71‑1 (MQ=255)
                                                          aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACg                                >  1:2668448/1‑51 (MQ=255)
                                                           gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                             <  1:3487039/53‑1 (MQ=255)
                                                           gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc           >  1:2278821/1‑71 (MQ=255)
                                                           gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc           >  1:59207/1‑71 (MQ=255)
                                                            tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                        >  1:3646852/1‑57 (MQ=255)
                                                            tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGcc          <  1:2126022/71‑1 (MQ=255)
                                                              ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc           >  1:1535743/1‑68 (MQ=255)
                                                              ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa         <  1:851055/70‑1 (MQ=255)
                                                              ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAg        <  1:3323253/71‑1 (MQ=255)
                                                                gACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                             >  1:765827/1‑48 (MQ=255)
                                                                 aCAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTg     <  1:764629/71‑1 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACTCTGCCGGGCCGTCAATAGAACGGTCAAc                          >  1:3013270/1‑49 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACGCTGCCTGGCCGTCAATAGAACGGTCAACAt                        <  1:145190/51‑1 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATg     >  1:2659163/1‑70 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATg     >  1:1830248/1‑70 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc    <  1:2277127/71‑1 (MQ=255)
                                                                  cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc    <  1:154459/71‑1 (MQ=255)
                                                                   aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa   <  1:3150288/71‑1 (MQ=255)
                                                                    aTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt  >  1:3006506/1‑71 (MQ=255)
                                                                      |                                                                    
ATTTTTTATCTTCCGTCAGCGCATTCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAT  >  minE/2425719‑2425857

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: