Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,128 G→C 9.2% T427T (ACG→ACC sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1280GC 9.2% 230.5 / 10.4 76T427T (ACG→ACCsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (40/29);  new base C (3/4);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAA  >  minE/1202058‑1202196
                                                                      |                                                                    
ttAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACg                                                                      <  1:3742065/71‑1 (MQ=255)
 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                     >  1:2325584/1‑71 (MQ=255)
 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                     >  1:3721231/1‑71 (MQ=255)
 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                     >  1:3726082/1‑71 (MQ=255)
    cgccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGt                                                                  <  1:1091951/71‑1 (MQ=255)
                ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                                                       >  1:2007158/1‑70 (MQ=255)
                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                    <  1:1463917/71‑1 (MQ=255)
                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                    <  1:3473038/71‑1 (MQ=255)
                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                    >  1:3828835/1‑71 (MQ=255)
                    ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                                  <  1:712785/71‑1 (MQ=12)
                    ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                  >  1:2208185/1‑71 (MQ=255)
                     tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                   <  1:1377931/53‑1 (MQ=255)
                               cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                       <  1:2452681/71‑1 (MQ=255)
                                aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                       <  1:1034642/70‑1 (MQ=255)
                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                               >  1:2380527/1‑61 (MQ=255)
                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                               >  1:2898195/1‑61 (MQ=255)
                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                     <  1:2126291/71‑1 (MQ=255)
                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                      >  1:489676/1‑70 (MQ=255)
                                  tGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCgg                                                        <  1:3799768/51‑1 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                                    >  1:667878/1‑37 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                                    >  1:584241/1‑37 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                                    >  1:3454234/1‑37 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                  >  1:2019384/1‑71 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                  >  1:1360214/1‑71 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                  <  1:2214042/71‑1 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAAAc                                  >  1:2187312/1‑71 (MQ=255)
                                      ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                   >  1:2645408/1‑68 (MQ=255)
                                      ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                 <  1:3132007/70‑1 (MQ=255)
                                      ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                                >  1:1763813/1‑71 (MQ=255)
                                       cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATg                                                                   <  1:2738225/35‑1 (MQ=255)
                                        gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                      >  1:375368/1‑63 (MQ=255)
                                         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                        >  1:326770/1‑60 (MQ=255)
                                         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                             <  1:2960531/71‑1 (MQ=255)
                                              ctcGCTGGATAAGGTCATCGAATCGATGTACGAAACCGGTAAGGa                                                  <  1:2101549/45‑1 (MQ=255)
                                              ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                  <  1:3929155/45‑1 (MQ=255)
                                              ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                                >  1:2464782/1‑47 (MQ=255)
                                                cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                      >  1:3246103/1‑71 (MQ=255)
                                                 gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAc                     <  1:2278689/71‑1 (MQ=255)
                                                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                    <  1:901792/71‑1 (MQ=255)
                                                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                    <  1:1667805/71‑1 (MQ=255)
                                                    ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                  <  1:4065938/39‑1 (MQ=255)
                                                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                            >  1:2266177/1‑44 (MQ=255)
                                                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                  <  1:190977/70‑1 (MQ=255)
                                                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                    <  1:3972827/68‑1 (MQ=255)
                                                       tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg               >  1:461654/1‑69 (MQ=9)
                                                       tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg               <  1:382465/71‑3 (MQ=9)
                                                        aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                              <  1:355144/39‑1 (MQ=255)
                                                        aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                          >  1:1684957/1‑59 (MQ=255)
                                                        aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg              >  1:2299572/1‑71 (MQ=255)
                                                         aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                            <  1:1213629/56‑1 (MQ=255)
                                                         aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg              <  1:969779/70‑1 (MQ=255)
                                                           gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                       <  1:4137919/43‑1 (MQ=255)
                                                            ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                       <  1:2694751/40‑1 (MQ=17)
                                                            ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                       <  1:310864/40‑1 (MQ=17)
                                                            tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTg          <  1:2297868/71‑1 (MQ=255)
                                                              aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt           >  1:3159503/1‑68 (MQ=255)
                                                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                       <  1:3073080/39‑1 (MQ=255)
                                                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                        >  1:3004455/1‑54 (MQ=255)
                                                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                   <  1:142465/59‑1 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                     >  1:298657/1‑39 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                           >  1:2227986/1‑49 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                           >  1:2223297/1‑49 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                          >  1:2673881/1‑50 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt     >  1:3157823/1‑71 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt     >  1:4121969/1‑71 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt     >  1:3602236/1‑71 (MQ=255)
                                                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACAGCGAAAc                          >  1:3177867/1‑50 (MQ=255)
                                                                 gAAACCATGTACGAAACAGGTAAAGATATGAACGcc                                        >  1:1676812/1‑36 (MQ=12)
                                                                 gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt     >  1:743552/1‑71 (MQ=17)
                                                                   aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa   <  1:2023539/71‑1 (MQ=255)
                                                                    aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                          >  1:2836474/1‑47 (MQ=255)
                                                                    aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc        >  1:562780/1‑65 (MQ=255)
                                                                    aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaa  <  1:1724198/71‑1 (MQ=255)
                                                                    aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa   >  1:3562572/1‑70 (MQ=255)
                                                                    aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa   >  1:1993506/1‑70 (MQ=255)
                                                                     cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa   >  1:4007118/1‑69 (MQ=255)
                                                                     cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa   >  1:2792915/1‑69 (MQ=255)
                                                                      gATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                        >  1:3024897/1‑47 (MQ=255)
                                                                      |                                                                    
TTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAA  >  minE/1202058‑1202196

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: