Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,892,325 T→C 18.2% intergenic (+23/‑117) ytfQ → / → ytfR predicted sugar transporter subunit/predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,892,3250TC18.2% 184.8 / 35.1 88intergenic (+23/‑117)ytfQ/ytfRpredicted sugar transporter subunit/predicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base C (8/8);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCG  >  minE/2892257‑2892395
                                                                    |                                                                      
tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTc                                                                      >  1:217898/1‑71 (MQ=255)
 aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCg                                                                     >  1:2844173/1‑71 (MQ=255)
  ccGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCg                                                                     <  1:1008660/70‑1 (MQ=255)
     cAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAgg                                                                 >  1:1318385/1‑71 (MQ=255)
     cAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAgg                                                                 <  1:3408682/71‑1 (MQ=255)
      aaaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGc                                                                <  1:1289023/71‑1 (MQ=255)
       aaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGcc                                                               <  1:3840403/71‑1 (MQ=255)
        aagaagAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCg                                                              <  1:950534/71‑1 (MQ=255)
               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGAt                                                           <  1:2492100/67‑1 (MQ=255)
               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAgg                                                       <  1:2562856/71‑1 (MQ=255)
               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAg                                                        <  1:1189055/70‑1 (MQ=255)
                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                      <  1:2093138/71‑1 (MQ=255)
                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                      >  1:3213144/1‑71 (MQ=255)
                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                      >  1:1839881/1‑71 (MQ=255)
                gtggcaaatgcctgatgcgctacgcttatcaggcctacatattcttgcgccATTCGTAGGCCGGATAAGgc                                                      >  1:1537681‑M2/52‑71 (MQ=255)
                gtggcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGAt                                                           <  1:2557604‑M2/15‑1 (MQ=255)
                 tgacaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGgc                                                      >  1:2898407‑M2/51‑70 (MQ=255)
                  agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcg                                                     >  1:2436020/1‑70 (MQ=255)
                  agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcg                                                     <  1:3308685/70‑1 (MQ=255)
                   gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAaggcgtt                                                   >  1:250992‑M2/49‑69 (MQ=255)
                    caaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTa                                                                   >  1:3572100‑M2/48‑54 (MQ=255)
                    caaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTc                                                  <  1:1854170‑M2/24‑1 (MQ=255)
                    caaatgcctgatgcgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATaaggcgtcc                                                  <  1:3576853‑M2/24‑4 (MQ=255)
                    agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATa                                                          <  1:3497309/63‑1 (MQ=255)
                     aaatgcctgatgcgctacgcttatcaggccttcagattcttgcgccATTCGTAGGCCGGATAAg                                                        >  1:1007612‑M2/47‑64 (MQ=255)
                     aaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAaggcgtt                                                   <  1:2376689‑M2/23‑3 (MQ=255)
                     aaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTtca                                                 >  1:115877‑M2/47‑70 (MQ=255)
                     aaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTtca                                                 >  1:2739671‑M2/47‑70 (MQ=255)
                      aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTcg                                                 <  1:297379/70‑1 (MQ=255)
                        tgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCg                                                                     <  1:2982146‑M2/5‑1 (MQ=255)
                        tgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCg                                                                     <  1:220638‑M2/5‑1 (MQ=255)
                        tgcctgatgcgctacgcttatcaggcctacagattcttgcgccATTCg                                                                     <  1:1858988‑M2/5‑1 (MQ=255)
                         tgcctgatgcgctacgcttatcaggcctacaattctcgcgccACTCg                                                                     >  1:3347542‑M2/43‑47 (MQ=255)
                         tgcctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                             >  1:212568‑M2/43‑71 (MQ=255)
                         tgcctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                             >  1:3638845‑M2/43‑71 (MQ=255)
                         gcctgatgcgctacgcttatcaggcctacagattcttgagccATTCGTAGGCCGGATAAGGCGTTCacgc                                              >  1:434558‑M2/43‑70 (MQ=255)
                         aagaaAAATATGTGGTATTGAGGGTTGCTATGCCTGATGCCGAtt                                                                       <  1:2897881/45‑1 (MQ=255)
                          agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                              >  1:2152367/1‑69 (MQ=255)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTa                                                                   >  1:1338455/1‑47 (MQ=255)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTa                                                                   >  1:2711116/1‑47 (MQ=255)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAg                                                        <  1:3722136/58‑1 (MQ=255)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAg                                                        <  1:423714/58‑1 (MQ=255)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAg                                                        <  1:428465/58‑1 (MQ=255)
                           ctgatgcgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCAcgccgc                                           <  1:2768907‑M2/31‑1 (MQ=255)
                           cctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                             <  1:1742409‑M2/29‑1 (MQ=255)
                             gatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                         >  1:3424724‑M2/39‑71 (MQ=255)
                              gatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATa                                                          <  1:3200186‑M2/16‑1 (MQ=255)
                              atgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                        >  1:1698281‑M2/38‑71 (MQ=255)
                              atgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                        >  1:1417802‑M2/38‑71 (MQ=255)
                              atgcgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATc                                        <  1:3367228‑M2/34‑1 (MQ=255)
                                tgcgctacgcttatcaggcctacaattctcgcgccACTCGTa                                                                   >  1:1545194‑M2/36‑42 (MQ=255)
                                tgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                      >  1:1593977‑M2/36‑71 (MQ=255)
                                 cgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                      <  1:2959088‑M2/36‑1 (MQ=255)
                                 cgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCgg                                     <  1:2360872‑M2/37‑1 (MQ=255)
                                  ttGGGTTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATa                                                          <  1:3541669/48‑1 (MQ=255)
                                   tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                    <  1:104092/70‑1 (MQ=255)
                                   ctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCGgca                                   <  1:1879157‑M2/39‑2 (MQ=255)
                                    tacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCggcat                                  >  1:2696707‑M2/32‑69 (MQ=255)
                                    tacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCggcag                                  >  1:1085164‑M2/32‑69 (MQ=255)
                                     tacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                 >  1:1709259‑M2/31‑71 (MQ=255)
                                     gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                  <  1:47971/70‑1 (MQ=255)
                                      ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                  <  1:712316/69‑1 (MQ=255)
                                      cgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCggcag                                  <  1:2136661‑M2/40‑3 (MQ=255)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCg                                                              >  1:1292980/1‑39 (MQ=255)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCg                                                              >  1:463623/1‑39 (MQ=255)
                                         ttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCa                                          >  1:3110272‑M2/27‑58 (MQ=255)
                                          ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                             >  1:2705674/1‑70 (MQ=255)
                                           tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgccgc                                           >  1:3850029/1‑55 (MQ=255)
                                           tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                  >  1:1298471/1‑64 (MQ=255)
                                           atcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTg                             <  1:585171‑M2/45‑1 (MQ=255)
                                           atcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTg                             <  1:339089‑M2/45‑1 (MQ=255)
                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCg                                      >  1:3652576/1‑59 (MQ=18)
                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATcc                                       <  1:1574127/58‑1 (MQ=255)
                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                              >  1:995263/1‑67 (MQ=255)
                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                           <  1:3526392/70‑1 (MQ=255)
                                            gAGGGGTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTcgcg                                               <  1:1045342/50‑1 (MQ=255)
                                             aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                         <  1:1563226/71‑1 (MQ=255)
                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                      <  1:3731858/71‑1 (MQ=255)
                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                             >  1:2165178/1‑47 (MQ=255)
                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                             >  1:3394226/1‑47 (MQ=255)
                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                <  1:933582/60‑1 (MQ=255)
                                                    acagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCggcaa                                  >  1:41356‑M2/16‑53 (MQ=255)
                                                     tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAGGcat                  >  1:1511977/1‑68 (MQ=255)
                                                      aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                >  1:3206027/1‑71 (MQ=255)
                                                       tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                 <  1:3604240/69‑1 (MQ=255)
                                                        gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                             <  1:3466736/56‑1 (MQ=255)
                                                             aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg         >  1:3843917/1‑71 (MQ=255)
                                                                  gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc    >  1:365581/1‑71 (MQ=255)
                                                                   aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc   >  1:72494/1‑71 (MQ=255)
                                                                    ttCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg  >  1:3310600/1‑71 (MQ=255)
                                                                    |                                                                      
TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCG  >  minE/2892257‑2892395

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: