Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,877,709 T→A 25.8% intergenic (+112/‑197) fklB → / → cycA FKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,877,7090TA25.8% 109.8 / 40.5 62intergenic (+112/‑197)fklB/cycAFKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (10/6);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCTA  >  minE/2877641‑2877778
                                                                    |                                                                     
cAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTAttttt                                                                     >  1:1598629/1‑71 (MQ=255)
   cAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtgt                                                                   >  1:517612/1‑70 (MQ=255)
       aaaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCt                                                              >  1:3710726/1‑71 (MQ=255)
        aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCt                                                              <  1:1791081/70‑1 (MQ=255)
        aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCt                                                              <  1:218082/70‑1 (MQ=255)
          cccTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGtt                                                           >  1:1603885/1‑71 (MQ=255)
             ttCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTc                                                         >  1:3050059/1‑70 (MQ=255)
               cAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTcgc                                                       <  1:1212365/70‑1 (MQ=255)
                aGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGt                                                     <  1:2225178/71‑1 (MQ=255)
                aGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                     >  1:3970668/1‑71 (MQ=255)
                  tttGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTa                                                   >  1:2045108/1‑71 (MQ=255)
                    tGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTAt                                                  >  1:1449040/1‑70 (MQ=255)
                      ccAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTc                                                         <  1:1961778/61‑1 (MQ=255)
                       cAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCAcc                                              >  1:172171/1‑71 (MQ=255)
                           cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGa                                          <  1:3637281/71‑1 (MQ=255)
                           cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGa                                          <  1:1682942/71‑1 (MQ=255)
                           cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTCGTGAGCTGTTTCGCGTTATCACCGTGa                                          <  1:1738159/71‑1 (MQ=255)
                            ggCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGa                                          <  1:1022774/70‑1 (MQ=255)
                              cgcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGatat                                       <  1:3080940/71‑1 (MQ=255)
                              cgcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGatat                                       <  1:2717779/71‑1 (MQ=255)
                              cgcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGatat                                       <  1:3736481/71‑1 (MQ=255)
                                cTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                    <  1:3628240/56‑1 (MQ=255)
                                cTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                     <  1:2242853/71‑1 (MQ=255)
                                cTATTTGGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATg                                      >  1:1863/1‑70 (MQ=255)
                                 tATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                     <  1:1401159/54‑1 (MQ=255)
                                  aTTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGata                                        >  1:3215279/1‑66 (MQ=255)
                                  aTTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                          >  1:724179/1‑64 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGtt                                                           <  1:3675735/46‑1 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATc                                                 <  1:1544117/56‑1 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:3231060/1‑70 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:1016952/1‑70 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 <  1:1090014/70‑1 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:239918/1‑70 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:3603268/1‑70 (MQ=255)
                                          ttGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTc                                                         <  1:2111259/41‑1 (MQ=255)
                                           tGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACtt                            >  1:3802192/1‑69 (MQ=255)
                                            gCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTAt                                                  >  1:2010108/1‑46 (MQ=255)
                                            gCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                             >  1:485373/1‑35 (MQ=255)
                                            gCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                             >  1:1464709/1‑35 (MQ=255)
                                               aaGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATc                                                 <  1:2346269/44‑1 (MQ=255)
                                                 gACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGa                                          >  1:1470635/1‑49 (MQ=255)
                                                    cccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCAc                                               >  1:2083138/1‑41 (MQ=255)
                                                    cccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaaaa                 >  1:3616581/1‑71 (MQ=255)
                                                     ccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGt                                            <  1:2057776/43‑1 (MQ=255)
                                                      cGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt               <  1:2780623/71‑1 (MQ=255)
                                                      cGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGATATCACCGTGATATGACACTCACTTTAAACATAAAAtt               <  1:3394056/71‑1 (MQ=255)
                                                      cGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt               <  1:469065/71‑1 (MQ=255)
                                                      cGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt               <  1:3383899/71‑1 (MQ=255)
                                                       gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCAGCGTGATATGACACTCACTTTAAACATAAAATTa              <  1:3792592/71‑1 (MQ=255)
                                                       gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt               >  1:3131546/1‑70 (MQ=255)
                                                       gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTa              <  1:1236299/71‑1 (MQ=255)
                                                       gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTa              <  1:3141727/71‑1 (MQ=255)
                                                       gTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt               >  1:1260076/1‑70 (MQ=255)
                                                       gTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAt                >  1:2129382/1‑69 (MQ=255)
                                                        tAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa             >  1:2699043/1‑71 (MQ=255)
                                                         aaGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa             >  1:589204/1‑70 (MQ=255)
                                                          aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGac                                    >  1:1473974/1‑46 (MQ=255)
                                                          aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa           >  1:96367/1‑71 (MQ=255)
                                                          aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa           >  1:1991821/1‑71 (MQ=255)
                                                             gtgtATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACttt                           >  1:2870420/1‑52 (MQ=255)
                                                             gtgtATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaaaa                 <  1:292139/62‑1 (MQ=255)
                                                              tgtATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATCACACTCACTTTAAACATAAAATTAACATTAg       <  1:2445451/71‑1 (MQ=255)
                                                                tATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAt          >  1:680340/1‑66 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGAcactcac                              >  1:218358/1‑42 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:1410271/1‑39 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                 >  1:1298536/1‑39 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATc    >  1:2727979/1‑68 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCTa  >  1:2846874/1‑70 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCTa  >  1:2364418/1‑70 (MQ=255)
                                                                    tttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCTa  >  1:1280091/1‑70 (MQ=255)
                                                                    |                                                                     
CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCTA  >  minE/2877641‑2877778

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: