Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,093,658 T→C T907A (ACA→GCA)  ydiJ ← predicted FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,093,6580TC96.5% 299.3 / 0.5 87T907A (ACA→GCA) ydiJpredicted FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (26/58);  total (29/58)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGGCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGTAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTAAACTCGCCA  >  minE/1093590‑1093725
                                                                    |                                                                   
tGGCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAg                                                                   >  1:182367/1‑71 (MQ=255)
tGGCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAg                                                                   >  1:1783812/1‑71 (MQ=255)
tGGCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTAACCTCTGACTGCAg                                                                   >  1:136198/1‑71 (MQ=255)
 ggCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGc                                                                  >  1:1591525/1‑71 (MQ=255)
 ggCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGc                                                                  >  1:1466630/1‑71 (MQ=255)
  gcgcACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGcc                                                                 <  1:1150731/71‑1 (MQ=255)
  gcgcACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGcc                                                                 <  1:2703159/71‑1 (MQ=255)
  gcgcACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGcc                                                                 <  1:2105639/71‑1 (MQ=255)
    gcACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCAc                                                               <  1:2698976/71‑1 (MQ=255)
      aCCCGGCAAGGCTGTAACTTCGGTACAGTGACCAAAGAAATACCATGTTTCACCGCTGACTGCAGCCACTg                                                             <  1:991656/71‑1 (MQ=255)
       cccGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTgg                                                            <  1:3103831/71‑1 (MQ=255)
        ccGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGc                                                           >  1:1568408/1‑71 (MQ=255)
        ccGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGc                                                           >  1:1628573/1‑71 (MQ=255)
         cGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCt                                                          <  1:2066898/71‑1 (MQ=255)
         cGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCt                                                          <  1:488000/71‑1 (MQ=255)
          ggCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCtg                                                         >  1:988358/1‑71 (MQ=255)
          ggCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCtg                                                         >  1:1481117/1‑71 (MQ=255)
          ggCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCtg                                                         >  1:2564117/1‑71 (MQ=255)
           gCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCtgt                                                        <  1:286102/71‑1 (MQ=255)
             aaGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGa                                                      <  1:3044913/71‑1 (MQ=255)
                  ggTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCaa                                                 <  1:939279/71‑1 (MQ=255)
                   gTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAg                                                <  1:1456616/71‑1 (MQ=255)
                   gTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAg                                                <  1:2498688/71‑1 (MQ=255)
                   gTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAg                                                <  1:2937534/71‑1 (MQ=255)
                    tAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCt                                                          <  1:1853493/60‑1 (MQ=255)
                          cGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGc                                                                  <  1:2245474/46‑1 (MQ=255)
                          cGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGc                                                                  <  1:244599/46‑1 (MQ=255)
                          cGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGc                                                                  <  1:144112/46‑1 (MQ=255)
                          cGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCtt                                         >  1:97096/1‑71 (MQ=255)
                          cGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCt                                          <  1:1270677/70‑1 (MQ=255)
                           ggTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCtt                                         <  1:970128/70‑1 (MQ=255)
                            gTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCtt                                         <  1:217078/69‑1 (MQ=255)
                              aCAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgc                                       >  1:2029921/1‑69 (MQ=255)
                               cAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCt                                                          <  1:2165286/49‑1 (MQ=255)
                               cAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccag                                    <  1:806010/71‑1 (MQ=255)
                               cAGTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccag                                    <  1:771433/71‑1 (MQ=255)
                                 gTGACCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   >  1:2292895/1‑70 (MQ=255)
                                    aCCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcatt                               >  1:2619076/1‑71 (MQ=255)
                                    aCCAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcatt                               >  1:1522010/1‑71 (MQ=255)
                                      cAAAGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTgcg                                            <  1:334708/56‑1 (MQ=255)
                                       aaaGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCTAGCcattcat                            <  1:2333060/71‑1 (MQ=255)
                                       aaaGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattcat                            <  1:422657/71‑1 (MQ=255)
                                       aaaGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattcat                            <  1:1377632/71‑1 (MQ=255)
                                        aaGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattcatt                           >  1:1564606/1‑71 (MQ=255)
                                        aaGAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattcat                            <  1:1823719/70‑1 (MQ=255)
                                          gAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcatt                               <  1:1456727/65‑1 (MQ=255)
                                          gAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcatt                               <  1:2095769/65‑1 (MQ=255)
                                          gAAATACCATGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcatt                               <  1:497029/65‑1 (MQ=255)
                                                  aTGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCattcattc                          <  1:352443/62‑1 (MQ=255)
                                                  aTGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCattcattc                          <  1:1905653/62‑1 (MQ=255)
                                                  aTGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCattcattc                          <  1:2507543/62‑1 (MQ=255)
                                                   tGATTCACCGCTGACTGTAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   >  1:1517794/1‑52 (MQ=255)
                                                   tGATTCACCGCTGACTGTAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   >  1:854721/1‑52 (MQ=255)
                                                   tGATTCACCGCTGACTGTAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   >  1:185974/1‑52 (MQ=255)
                                                   tGATTCACCGCTGACTGCAGCCACTGGCTGTGACt                                                    >  1:258734/1‑35 (MQ=255)
                                                   tGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGa                <  1:1516698/71‑1 (MQ=255)
                                                   tGATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGa                <  1:2117801/71‑1 (MQ=255)
                                                    gATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattca                             <  1:1878989/57‑1 (MQ=255)
                                                    gATTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattca                             <  1:2772101/57‑1 (MQ=255)
                                                     aTTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattc                              <  1:3003694/55‑1 (MQ=255)
                                                     aTTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattc                              <  1:3192884/55‑1 (MQ=255)
                                                     aTTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattc                              <  1:2662503/55‑1 (MQ=255)
                                                     aTTCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCcattc                              <  1:1991843/55‑1 (MQ=255)
                                                      ttCACCGCTGACTGCAGCCTCTGGCTGTGACTCAAg                                                >  1:2351158/1‑36 (MQ=37)
                                                      ttCACCGCTGACTGCAGCCACTGGCTGTGACTCAAg                                                >  1:886005/1‑36 (MQ=255)
                                                      ttCACCGCTGACTGCAGCCACTGGCTGTGACTCAAg                                                >  1:1892006/1‑36 (MQ=255)
                                                      ttCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccggc                                   <  1:848671/49‑1 (MQ=255)
                                                      ttCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   <  1:973571/49‑1 (MQ=255)
                                                      ttCACCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTgccagc                                   <  1:1055923/49‑1 (MQ=255)
                                                         aCCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:1818549/70‑1 (MQ=255)
                                                         aCCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:2898612/70‑1 (MQ=255)
                                                         aCCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           >  1:249723/1‑70 (MQ=255)
                                                         aCCGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:1446516/70‑1 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCtt                                         <  1:1405713/39‑1 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCtt                                         <  1:1456601/39‑1 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGc                        >  1:646386/1‑56 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGc                        >  1:1304889/1‑56 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGc                        >  1:1234388/1‑56 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGc                        >  1:2524926/1‑56 (MQ=255)
                                                          ccGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTaaa         <  1:2340064/71‑1 (MQ=255)
                                                           cGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:1933032/68‑1 (MQ=255)
                                                           cGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:2059478/68‑1 (MQ=255)
                                                           cGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTa           <  1:1890404/68‑1 (MQ=255)
                                                           cGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTAAAc        >  1:705980/1‑71 (MQ=255)
                                                           cGCTGACTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTAAAc        <  1:322215/71‑1 (MQ=255)
                                                                 cTGCAGCCTCTGGCTGTGACTCAAGTGCGCTTGCCAGCcattca                             <  1:559358/44‑1 (MQ=38)
                                                                 cTGCAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTAAACTCGCCa  <  1:1851526/71‑1 (MQ=255)
                                                                    |                                                                   
TGGCGCACCCGGCAAGGCGGTAACTTCGGTACAGTGACCAAAGAAATACCATGATTCACCGCTGACTGTAGCCACTGGCTGTGACTCAAGTGCGCTTGCCAGCCATTCATTCGCCAGTAAGACGTTAAACTCGCCA  >  minE/1093590‑1093725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: