Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,008,960 G→A 11.3% L98L (CTG→CTA folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,008,9600GA11.3% 204.9 / 14.7 71L98L (CTG→CTAfolMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base A (4/4);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCATGCGCTGGA  >  minE/1008895‑1009028
                                                                 |                                                                    
cATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGa                                                                 >  1:2022046/1‑71 (MQ=255)
cATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGa                                                                 >  1:2505838/1‑71 (MQ=255)
cATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGa                                                                 >  1:1538924/1‑71 (MQ=255)
 aTGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTAGCCGAc                                                                >  1:928713/1‑71 (MQ=255)
 aTGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTAGCCGAc                                                                >  1:1733008/1‑71 (MQ=255)
  tGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                               <  1:3442281/71‑1 (MQ=255)
         cGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGc                                                        >  1:3096909/1‑71 (MQ=255)
         cGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGc                                                        >  1:2098179/1‑71 (MQ=255)
             cTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg                                                     <  1:3140984/70‑1 (MQ=255)
             cTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg                                                     <  1:2749661/70‑1 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                                  <  1:3735903/70‑1 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                                  >  1:2752781/1‑70 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                                  >  1:1701222/1‑70 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                                  <  1:1041558/70‑1 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                                  <  1:629509/70‑1 (MQ=255)
                ttttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCa                                                   <  1:2174603/69‑1 (MQ=255)
                 tttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTa                                                             >  1:892153/1‑58 (MQ=255)
                  ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGa                                                                 >  1:3296841/1‑53 (MQ=255)
                  ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGa                                                                 >  1:2676694/1‑53 (MQ=255)
                  ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgat                                               >  1:792688/1‑71 (MQ=255)
                  ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                                <  1:1717110/70‑1 (MQ=255)
                    gCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:2105339/1‑70 (MQ=255)
                     cATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:653997/1‑69 (MQ=255)
                     cATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGc                                             <  1:3246801/70‑1 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:3115346/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:10702/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:577385/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:3478766/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:2365689/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:1963250/1‑66 (MQ=255)
                        aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatgatg                                              >  1:1875320/1‑66 (MQ=255)
                         aCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc                                        >  1:3701897/1‑71 (MQ=255)
                               gTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                                  <  1:2319177/71‑1 (MQ=255)
                                tGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                                  >  1:1630902/1‑70 (MQ=255)
                                  cGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                                <  1:3233342/70‑1 (MQ=255)
                                   gTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                                  <  1:2473100/67‑1 (MQ=255)
                                    tGGATGGCGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGCTTGCATGATGCAGa                                          <  1:760720/58‑1 (MQ=255)
                                      gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAccc                           >  1:2256862/1‑71 (MQ=255)
                                      gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAccc                           >  1:3355078/1‑71 (MQ=255)
                                      gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAccc                           >  1:1974251/1‑71 (MQ=255)
                                      gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATACACGTTAATAccc                           >  1:212901/1‑71 (MQ=255)
                                          gCGGAAAAACCGGGTGCGCCACTTGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCAt                        <  1:3260704/70‑1 (MQ=255)
                                          gCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCAt                        <  1:3055029/70‑1 (MQ=255)
                                          gCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCAt                        <  1:1757784/70‑1 (MQ=255)
                                          gCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCAt                        <  1:460858/70‑1 (MQ=255)
                                          gCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCAt                        <  1:2922072/70‑1 (MQ=255)
                                           cGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAc                      <  1:881431/71‑1 (MQ=255)
                                           cGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAc                      <  1:1508201/71‑1 (MQ=255)
                                           cGGAAAAACCGGGTGCGCCACTAGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAc                      <  1:2685411/71‑1 (MQ=255)
                                            ggAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAcc                     >  1:3502031/1‑71 (MQ=255)
                                            ggAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAcc                     >  1:2369780/1‑71 (MQ=255)
                                             gAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGc                                             <  1:2535201/46‑1 (MQ=255)
                                             gAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc                                        <  1:1911871/51‑1 (MQ=255)
                                             gAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc                                        <  1:1294290/51‑1 (MQ=255)
                                               aaaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCa                                      >  1:1110766/1‑51 (MQ=255)
                                               aaaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGc                  >  1:2066756/1‑71 (MQ=255)
                                                aaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAcc                            >  1:1515652/1‑60 (MQ=255)
                                                aaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAcc                            >  1:2923605/1‑60 (MQ=255)
                                                aaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCt                 >  1:2614677/1‑71 (MQ=255)
                                                 aaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTc                <  1:1628381/71‑1 (MQ=255)
                                                 aaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTc                <  1:1680472/71‑1 (MQ=255)
                                                  aCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTc                >  1:320394/1‑70 (MQ=255)
                                                  aCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCa               <  1:1617249/71‑1 (MQ=255)
                                                      ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                                <  1:88400/50‑1 (MQ=255)
                                                      ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCa           >  1:2528083/1‑71 (MQ=255)
                                                      ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCa           >  1:2913333/1‑71 (MQ=255)
                                                      ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCa           >  1:2222451/1‑71 (MQ=255)
                                                      ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCa           <  1:888050/71‑1 (MQ=255)
                                                       gTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCAt          >  1:1108117/1‑71 (MQ=255)
                                                               cTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCATGCGCTGGa  <  1:377234/71‑1 (MQ=255)
                                                               cTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCATGCGCTGGa  <  1:3414497/71‑1 (MQ=255)
                                                                 ggCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTc                <  1:2796979/55‑1 (MQ=255)
                                                                 ggCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCATGCGCTGGa  >  1:770773/1‑69 (MQ=255)
                                                                 |                                                                    
CATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAACCATGCGCTGGA  >  minE/1008895‑1009028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: