Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,708,513 C→T 9.8% intergenic (+51/‑239) trpS → / → yhfZ tryptophanyl‑tRNA synthetase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,708,5130CT 9.8% 140.2 / 4.8 51intergenic (+51/‑239)trpS/yhfZtryptophanyl‑tRNA synthetase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/28);  new base T (3/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCAT  >  minE/2708446‑2708580
                                                                   |                                                                   
ttGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                  >  1:797493/1‑71 (MQ=255)
 tgtgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                  <  1:5037687/70‑1 (MQ=255)
   tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                <  1:610389/70‑1 (MQ=255)
   tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                               <  1:1204911/71‑1 (MQ=255)
   tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                               <  1:736673/71‑1 (MQ=255)
   tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                               <  1:186081/71‑1 (MQ=255)
      cGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                <  1:1095803/67‑1 (MQ=255)
        aaGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                               <  1:2277731/66‑1 (MQ=255)
           ccGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                        >  1:4915527/1‑70 (MQ=255)
             gTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                      >  1:5523793/1‑70 (MQ=255)
              tAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                     >  1:23368/1‑70 (MQ=255)
                aTTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                      <  1:5000217/67‑1 (MQ=255)
                 ttCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAg                                                 <  1:4359624/71‑1 (MQ=255)
                       aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                              >  1:4534983/1‑52 (MQ=11)
                       aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                              >  1:5060043/1‑52 (MQ=11)
                       aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                              >  1:3300924/1‑52 (MQ=11)
                        ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                 >  1:5399815/1‑48 (MQ=255)
                        ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGc                                           <  1:56510/70‑1 (MQ=255)
                        ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGc                                           <  1:1674458/70‑1 (MQ=255)
                          gCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGcaca                                        >  1:4685801/1‑71 (MQ=255)
                            tGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                               <  1:859954/46‑1 (MQ=255)
                             aacAACTTGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCgg                                                          <  1:5467278/48‑1 (MQ=11)
                                 acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                                           >  1:3592912/3‑45 (MQ=11)
                                 acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                                           >  1:2113157/3‑45 (MQ=11)
                                 acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                                           >  1:3127233/3‑45 (MQ=11)
                                 aaTTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCa                                 <  1:1939108/71‑1 (MQ=255)
                                  cttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                                           <  1:2178324/43‑1 (MQ=11)
                                   ttGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAg                               <  1:421498/71‑1 (MQ=255)
                                    tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTaa                                              >  1:580117/1‑55 (MQ=38)
                                     gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAgg                              >  1:4676054/1‑70 (MQ=255)
                                      cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGaa                            <  1:1651142/71‑1 (MQ=255)
                                        cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAAt                          <  1:3270070/71‑1 (MQ=255)
                                         tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                         >  1:1982313/1‑71 (MQ=255)
                                           tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                         <  1:4888478/69‑1 (MQ=255)
                                                aTGCCGGATGCGGCGTGAACTCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGa                             <  1:1860145/60‑1 (MQ=255)
                                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGa                             <  1:3954048/59‑1 (MQ=255)
                                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGa                >  1:2272304/1‑70 (MQ=255)
                                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACat                                       <  1:2739686/46‑1 (MQ=255)
                                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATa              <  1:5031994/71‑1 (MQ=255)
                                                        tGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCg          >  1:3603607/1‑71 (MQ=255)
                                                         gcggcgTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATa              >  1:2940464/1‑66 (MQ=255)
                                                          cggcgTGAACGCCTTATCCGGCCTACAAAGTAATg                                            >  1:4387541/1‑35 (MQ=255)
                                                            gcgTGAACGCGTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAgc      <  1:1384881/71‑1 (MQ=255)
                                                            gcgTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTg                                   <  1:4038131/42‑1 (MQ=255)
                                                             cgTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGt                        >  1:3409375/1‑52 (MQ=255)
                                                             cgTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAgcg     <  1:4653431/71‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                <  1:1975071/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                <  1:66273/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                <  1:1383325/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                <  1:1146221/42‑1 (MQ=255)
                                                                 aaCGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCAt  >  1:3370133/1‑70 (MQ=255)
                                                                   cGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                         >  1:299035/1‑45 (MQ=255)
                                                                   cGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                         <  1:355489/45‑1 (MQ=255)
                                                                   |                                                                   
TTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCAT  >  minE/2708446‑2708580

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: