Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,225,948 G→C 21.6% V129V (GTG→GTC dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,225,9480GC21.6% 177.2 / 46.3 88V129V (GTG→GTCdusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base G (28/41);  new base C (13/6);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 3.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAA  >  minE/2225878‑2226017
                                                                      |                                                                     
aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg                                                                       <  1:3783823/71‑1 (MQ=255)
 tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg                                                                       >  1:371905/1‑70 (MQ=255)
     aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  <  1:1197369/71‑1 (MQ=255)
     aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  <  1:824064/71‑1 (MQ=255)
     aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  <  1:2968828/71‑1 (MQ=255)
       gCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  <  1:3736900/69‑1 (MQ=255)
         tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg                                                              >  1:249402/1‑71 (MQ=255)
           gCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa                                                            <  1:846568/71‑1 (MQ=255)
             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg                                                              <  1:149032/67‑1 (MQ=255)
             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg                                                              <  1:3351605/67‑1 (MQ=255)
              ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGt                                                         <  1:1773993/71‑1 (MQ=255)
                 tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTcc                                                      >  1:474043/1‑71 (MQ=255)
                  cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCt                                                     >  1:1836973/1‑71 (MQ=255)
                   aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCt                                                     <  1:1100646/70‑1 (MQ=255)
                     cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGaa                                                                     >  1:3746082/1‑52 (MQ=255)
                        tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGAAGTTc                                                       <  1:1390836/63‑1 (MQ=255)
                          ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAg                                                                >  1:844470/1‑52 (MQ=255)
                          ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa                                                            >  1:109602/1‑56 (MQ=255)
                          ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCaa                                                                     >  1:737736/1‑47 (MQ=255)
                          ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCaa                                                                     >  1:4067667/1‑47 (MQ=255)
                           tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg                                                                       <  1:3999125/44‑1 (MQ=255)
                           tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGt                                                         <  1:3767477/58‑1 (MQ=255)
                           tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCt                                             >  1:3194082/1‑70 (MQ=255)
                            gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg                                            >  1:2407788/1‑70 (MQ=255)
                            gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg                                            <  1:3061687/70‑1 (MQ=255)
                            gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                            >  1:4089697/1‑70 (MQ=255)
                              aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  <  1:302763/46‑1 (MQ=255)
                              aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg                                            <  1:4049192/68‑1 (MQ=255)
                              aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg                                            <  1:3317275/68‑1 (MQ=255)
                                tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa                                                            <  1:3367987/50‑1 (MQ=255)
                                tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa                                                 <  1:3619858/61‑1 (MQ=255)
                                tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                                   <  1:563923/43‑1 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGATTGCAGTGGACGTTCCTGTTACCCTGAAGAtt                                      <  1:1251302/70‑1 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa                                                 <  1:603562/59‑1 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa                                                 <  1:904980/59‑1 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                     >  1:552476/1‑71 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                     >  1:3070399/1‑71 (MQ=255)
                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGt                                                   >  1:1257279/1‑57 (MQ=255)
                                   ccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCg                                    >  1:961502/1‑71 (MQ=255)
                                    cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg                                                                       <  1:2527596/35‑1 (MQ=255)
                                    cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa                                        <  1:978580/66‑1 (MQ=255)
                                    cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc                                   <  1:480107/71‑1 (MQ=255)
                                     ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc                                   >  1:1019297/1‑70 (MQ=255)
                                      gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc                                                                  >  1:2336690/1‑38 (MQ=255)
                                      gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGaa                                          <  1:862912/62‑1 (MQ=255)
                                      gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGaa                                          <  1:1559226/62‑1 (MQ=255)
                                       aTGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCGTGAAGATTCGCAcc                                >  1:3568028/1‑71 (MQ=255)
                                        tGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAcc                                               >  1:2867657/1‑55 (MQ=255)
                                        tGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                               >  1:872906/1‑71 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGtt                                                        <  1:3296814/44‑1 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGtt                                                        <  1:4118384/44‑1 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                               <  1:2935189/69‑1 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                               <  1:638173/69‑1 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc                             <  1:2194316/71‑1 (MQ=255)
                                          tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc                             >  1:1892846/1‑71 (MQ=255)
                                           cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCt                            >  1:4072572/1‑71 (MQ=255)
                                            gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGa                                           <  1:26493/55‑1 (MQ=255)
                                            gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                           <  1:371425/71‑1 (MQ=255)
                                            gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                           <  1:2442608/71‑1 (MQ=255)
                                            gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGATg                           >  1:1704525/1‑71 (MQ=255)
                                            gTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGtt                                                  >  1:834632/1‑48 (MQ=255)
                                              tAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                           >  1:2925710/1‑69 (MQ=255)
                                               aaaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                         >  1:902497/1‑70 (MQ=255)
                                               aaaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                         <  1:948761/70‑1 (MQ=255)
                                               aaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                         >  1:1040573/1‑70 (MQ=255)
                                                aaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                        <  1:4113760/70‑1 (MQ=255)
                                                   cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAcc                                               >  1:1707126/1‑44 (MQ=255)
                                                   cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa                                        <  1:861811/51‑1 (MQ=255)
                                                   cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                    >  1:2639476/1‑71 (MQ=255)
                                                   cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                    >  1:2734110/1‑71 (MQ=255)
                                                    gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa                                        >  1:3753712/1‑50 (MQ=255)
                                                    gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                   <  1:4043621/71‑1 (MQ=255)
                                                    gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                   <  1:894718/71‑1 (MQ=255)
                                                    gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                   <  1:3449094/71‑1 (MQ=255)
                                                    gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                    >  1:3700855/1‑70 (MQ=255)
                                                     aTCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt                                                   <  1:706093/38‑1 (MQ=255)
                                                      tCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt                                                   >  1:2208856/1‑37 (MQ=255)
                                                      tCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAc                      >  1:2570900/1‑66 (MQ=255)
                                                        cTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt                                                   <  1:1812174/35‑1 (MQ=255)
                                                          tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac              >  1:3741284/1‑70 (MQ=255)
                                                          tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca               >  1:3100097/1‑69 (MQ=255)
                                                          tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc             >  1:759367/1‑71 (MQ=255)
                                                           aCCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCg            >  1:3348407/1‑71 (MQ=255)
                                                            ccGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt           <  1:1447235/71‑1 (MQ=255)
                                                            ccGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt           <  1:3423814/71‑1 (MQ=255)
                                                            ccGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt           <  1:3193126/71‑1 (MQ=255)
                                                             cGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTa          <  1:2980963/71‑1 (MQ=255)
                                                                  tCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                       <  1:3172347/37‑1 (MQ=255)
                                                                   cGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc                             >  1:2270427/1‑46 (MQ=255)
                                                                   cGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGc     >  1:766059/1‑70 (MQ=255)
                                                                    gTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCg    >  1:1238930/1‑70 (MQ=255)
                                                                     tGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa  <  1:3317308/71‑1 (MQ=255)
                                                                      |                                                                     
ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAA  >  minE/2225878‑2226017

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: