Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,947,541 G→A 15.0% F89F (TTC→TTT dnaT ← DNA biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,947,5410GA15.0% 96.3 / 8.4 40F89F (TTC→TTTdnaTDNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (21/13);  new base A (2/4);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACA  >  minE/2947472‑2947610
                                                                     |                                                                     
tAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCg                                                                       >  1:1027677/1‑70 (MQ=255)
 aGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCaaa                                                                     >  1:3752166/1‑71 (MQ=255)
           cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCAtt                                                           >  1:990787/1‑70 (MQ=255)
           cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATg                                                           >  1:3236815/1‑71 (MQ=255)
                    ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGc                                                   >  1:2488774/1‑70 (MQ=255)
                       ccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACg                                               >  1:1346603/1‑71 (MQ=255)
                         aGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACggg                                             <  1:555904/71‑1 (MQ=255)
                          ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCa                                                  <  1:1840051/65‑1 (MQ=255)
                           gCGGATAAATTCGGCATCGGGTTGCCAGTCCGGATACATGGCAAATTTCCCCATGGGTACAGCAACGGTAg                                           <  1:3563904/71‑1 (MQ=255)
                            cGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGc                                          >  1:2492069/1‑71 (MQ=255)
                                 aaaaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCAt                                                            >  1:155886/1‑48 (MQ=255)
                                   aaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCa                                                             <  1:1294936/45‑1 (MQ=255)
                                    aTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCa                                                             >  1:1840092/1‑44 (MQ=255)
                                     tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTa                                                       >  1:2693671/1‑49 (MQ=255)
                                     tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTa                                                       >  1:3781657/1‑49 (MQ=255)
                                     tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTa                                                       >  1:125885/1‑49 (MQ=255)
                                     tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTa                                                       >  1:1292226/1‑49 (MQ=255)
                                       ggCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGtt                               <  1:2579747/71‑1 (MQ=255)
                                       ggCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTACGGtt                               <  1:2573320/71‑1 (MQ=255)
                                          aTCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCgg                                 <  1:3848978/66‑1 (MQ=255)
                                           tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCt                            >  1:2045869/1‑70 (MQ=255)
                                           tCGGGTTGCCAGTCCGGATACATGGCAAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGc                                   <  1:3739293/63‑1 (MQ=255)
                                             gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                          >  1:758188/1‑70 (MQ=255)
                                             gggTTGCCAGTCCGGATACATGGCAAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCgg                                 >  1:895553/1‑63 (MQ=255)
                                              ggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGc                                   <  1:3649385/60‑1 (MQ=255)
                                                ttGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATaga                      <  1:1871676/71‑1 (MQ=255)
                                                ttGCCAGTCCGGATACATGGCAAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATaga                      <  1:2465640/71‑1 (MQ=255)
                                                 tGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCgg                                 >  1:644068/1‑59 (MQ=255)
                                                 tGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCg                                  >  1:930045/1‑58 (MQ=255)
                                                 tGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATagag                     >  1:3080843/1‑71 (MQ=255)
                                                     aGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATagag                     >  1:2399594/1‑67 (MQ=255)
                                                     aGTCCGGATACATGGCAAATTTCCCCATGGGTACAGCAACGGGa                                            >  1:788113/1‑44 (MQ=255)
                                                      gTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGAt                >  1:2107637/1‑71 (MQ=255)
                                                       tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTcc                             >  1:2955875/1‑57 (MQ=255)
                                                          ggATACATGGCAAATTTCCCCATGGGTACAGCAACGGGAGCgg                                        <  1:2527920/43‑1 (MQ=255)
                                                            aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCt                            <  1:25000/53‑1 (MQ=255)
                                                            aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTc            >  1:3258203/1‑69 (MQ=255)
                                                              aCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGt             <  1:917525/66‑1 (MQ=255)
                                                               cATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCt                                      <  1:1948553/40‑1 (MQ=255)
                                                                 tGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTg                                    <  1:3825032/40‑1 (MQ=255)
                                                                  ggCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTc                              >  1:604028/1‑45 (MQ=255)
                                                                    cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACa  <  1:1181797/71‑1 (MQ=255)
                                                                     |                                                                     
TAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACA  >  minE/2947472‑2947610

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: