Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,094 C→A 12.5% A416E (GCA→GAA) ‡ sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,0940CA12.5% 233.2 / 11.5 78A416E (GCA→GAA) ‡sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base A (3/7);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATA  >  minE/1202024‑1202163
                                                                      |                                                                     
aTTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGc                                                                       <  1:4469686/71‑1 (MQ=255)
   gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:1055199/1‑69 (MQ=255)
   gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:5573887/1‑69 (MQ=255)
   gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCAcc                                                                    >  1:580349/1‑71 (MQ=255)
       gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                                <  1:58668/71‑1 (MQ=255)
       gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                                <  1:5164087/71‑1 (MQ=255)
        tAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgcg                                                               >  1:880853/1‑71 (MQ=255)
            gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:3333967/1‑71 (MQ=255)
            gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:3739672/1‑71 (MQ=255)
             ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:4229551/1‑70 (MQ=255)
             ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                          >  1:2211402/1‑71 (MQ=255)
               aTTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                          <  1:845554/69‑1 (MQ=255)
                 tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                       >  1:3480006/1‑70 (MQ=255)
                 tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                       >  1:611908/1‑70 (MQ=255)
                      ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:4869360/1‑50 (MQ=255)
                      ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                                 >  1:4823271/1‑71 (MQ=255)
                      ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                                 >  1:4865052/1‑71 (MQ=255)
                       tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:4900608/1‑71 (MQ=255)
                       tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:3969631/1‑71 (MQ=255)
                       tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:1298988/1‑71 (MQ=255)
                         tgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTc                                              <  1:3333494/71‑1 (MQ=255)
                           aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCa                                             >  1:1598500/1‑70 (MQ=255)
                           aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCAt                                            <  1:1483214/71‑1 (MQ=255)
                            aGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATc                                           <  1:1486707/71‑1 (MQ=255)
                               cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                        <  1:1098307/71‑1 (MQ=255)
                               cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                        <  1:1944295/71‑1 (MQ=255)
                                  ttAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                                      >  1:3092243/1‑70 (MQ=255)
                                   tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    >  1:573100/1‑71 (MQ=255)
                                    aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    <  1:2635090/70‑1 (MQ=255)
                                    aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   >  1:4061821/1‑71 (MQ=255)
                                    aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    <  1:2331851/70‑3 (MQ=9)
                                    aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    <  1:2338470/70‑3 (MQ=9)
                                    aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    >  1:1675178/1‑68 (MQ=9)
                                     aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   <  1:1416027/70‑1 (MQ=255)
                                     aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   <  1:5636342/70‑1 (MQ=255)
                                     aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                  >  1:1591409/1‑71 (MQ=255)
                                     aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                  <  1:214850/71‑1 (MQ=255)
                                      cgccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                                      <  1:1318944/66‑1 (MQ=255)
                                        ccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    >  1:2896731/1‑66 (MQ=255)
                                         cgcacgTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACg                              <  1:2094673/71‑1 (MQ=12)
                                          gcgcgGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                             <  1:3272937/71‑1 (MQ=255)
                                             cgtaTGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                          <  1:3611663/68‑1 (MQ=12)
                                               gATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCg                        >  1:3287383/1‑71 (MQ=255)
                                                  ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                      >  1:1703545/1‑70 (MQ=255)
                                                   gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:602788/70‑1 (MQ=255)
                                                   gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:5028452/70‑1 (MQ=255)
                                                   gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:374468/70‑1 (MQ=255)
                                                    ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCg                                          <  1:1429092/48‑1 (MQ=255)
                                                     gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                  <  1:3752269/70‑3 (MQ=12)
                                                     gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                  <  1:155195/70‑3 (MQ=12)
                                                      ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa                    >  1:1338565/1‑68 (MQ=16)
                                                      ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa                    >  1:4658923/1‑68 (MQ=16)
                                                      ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                 >  1:5130055/1‑71 (MQ=255)
                                                       tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                               <  1:5138410/56‑1 (MQ=255)
                                                       tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                               <  1:279462/56‑1 (MQ=255)
                                                         cgccgtACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata               <  1:1741280/70‑3 (MQ=9)
                                                         cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                             >  1:4681059/1‑56 (MQ=255)
                                                         cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt              <  1:5233609/71‑1 (MQ=255)
                                                          gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:1053086/71‑1 (MQ=255)
                                                          gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:3386637/71‑1 (MQ=255)
                                                          gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:2945933/71‑1 (MQ=255)
                                                           ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:485031/71‑1 (MQ=255)
                                                           ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:740483/71‑1 (MQ=255)
                                                           ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:2237289/71‑1 (MQ=255)
                                                           ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:3133687/71‑1 (MQ=255)
                                                            cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                       >  1:4772866/1‑43 (MQ=255)
                                                            cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                       >  1:1721249/1‑43 (MQ=255)
                                                            cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:122387/69‑1 (MQ=255)
                                                            cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa           >  1:2482802/1‑71 (MQ=255)
                                                              cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg         >  1:3152718/1‑71 (MQ=255)
                                                              cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg         >  1:2128450/1‑71 (MQ=255)
                                                                ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       <  1:2282120/71‑1 (MQ=255)
                                                                ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       <  1:1873247/71‑1 (MQ=255)
                                                                 cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       >  1:4645873/1‑70 (MQ=255)
                                                                  aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTa                                >  1:3994514/1‑44 (MQ=255)
                                                                   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa    <  1:5386880/71‑1 (MQ=255)
                                                                   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa    <  1:2388045/71‑1 (MQ=255)
                                                                    tGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAt   <  1:3202003/71‑1 (MQ=255)
                                                                     gCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa  >  1:1085806/1‑71 (MQ=255)
                                                                      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa  <  1:588897/70‑1 (MQ=255)
                                                                      |                                                                     
ATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATA  >  minE/1202024‑1202163

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: ‑