Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,708,514 G→T 11.3% intergenic (+52/‑238) 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,5140GT11.3% 153.1 / 5.9 53intergenic (+52/‑238)trpS/yhfZtryptophanyl‑tRNA synthetase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base T (3/3);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

TTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCATCAGG  >  minE/2708446‑2708584
                                                                    |                                                                      
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:1674458/70‑1 (MQ=255)
                        ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGc                                               <  1:56510/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:3127233/3‑45 (MQ=11)
                                 acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                                               >  1:2113157/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:1146221/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                    <  1:1383325/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                    <  1:66273/42‑1 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCaa                                    <  1:1975071/42‑1 (MQ=255)
                                                                 aaCGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCAt      >  1:3370133/1‑70 (MQ=255)
                                                                   cGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                             <  1:355489/45‑1 (MQ=255)
                                                                   cGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                             >  1:299035/1‑45 (MQ=255)
                                                                    gCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg                             <  1:16924/44‑1 (MQ=255)
                                                                    gCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCATCAgg  >  1:813129/1‑71 (MQ=255)
                                                                    gCCTTATACGGCCTACAAAGTAATGCACATATTGCAAGGaaa                               >  1:1495741/1‑42 (MQ=255)
                                                                    |                                                                      
TTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATGTAGGCCTGATAAGCGTAGCGCATCAGG  >  minE/2708446‑2708584

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: