Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,201,060	A→G	11.7%	A71A (GCA→GCG)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,201,060	0	A	G	11.7%	204.0 / 15.8	77	A71A (GCA→GCG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base A (27/41);  new base G (5/4);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAC  >  minE/1200993‑1201130                                                                    |                                                                       tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCCTGCAg                                                                       <  1:2544885/69‑1 (MQ=255)  ttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                                                     >  1:395702/1‑70 (MQ=255)    aTGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                                                  >  1:3449968/1‑71 (MQ=255)     tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                                                  <  1:1349149/70‑1 (MQ=255)     tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                 <  1:3669786/71‑1 (MQ=255)       gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                <  1:2893648/70‑1 (MQ=255)       gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                <  1:718191/70‑1 (MQ=255)        ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                              <  1:260318/71‑1 (MQ=255)        ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                              <  1:2552819/71‑1 (MQ=255)          cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                            >  1:3847740/1‑71 (MQ=255)             tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                          >  1:3577907/1‑70 (MQ=255)               tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                 <  1:3503918/61‑1 (MQ=255)               tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                       >  1:3246124/1‑71 (MQ=255)               tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                       >  1:1524997/1‑71 (MQ=255)                 gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                                                      <  1:1479701/70‑1 (MQ=255)                 gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                     <  1:2339332/71‑1 (MQ=255)                 gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                     <  1:191605/71‑1 (MQ=255)                      ccggTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGat                                                <  1:2709983/69‑1 (MQ=255)                      cGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGat                                                <  1:2131114/71‑1 (MQ=255)                        ggTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                              >  1:1541547/1‑71 (MQ=255)                         gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                              >  1:3463086/1‑70 (MQ=255)                         gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                             >  1:3458611/1‑71 (MQ=255)                         gTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                              >  1:3475882/1‑70 (MQ=255)                          tAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCg                                            <  1:3424658/71‑1 (MQ=255)                          tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                            <  1:1739371/55‑1 (MQ=255)                          tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATc                                             >  1:2437087/1‑70 (MQ=255)                           aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                            >  1:3160451/1‑70 (MQ=255)                           aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                           >  1:1698745/1‑71 (MQ=255)                           aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                           >  1:28292/1‑71 (MQ=255)                           aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                           >  1:2807012/1‑71 (MQ=255)                            aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                     <  1:2795167/60‑1 (MQ=255)                            aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                                   >  1:221149/1‑62 (MQ=255)                            aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                           >  1:2524524/1‑70 (MQ=255)                            aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                          <  1:1115253/71‑1 (MQ=255)                            aaGGCCACCACACCGATATCCCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                     <  1:1499202/60‑1 (MQ=255)                             aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                         <  1:2122150/71‑3 (MQ=255)                             aGGCCACCACACCGATATCGCCA‑TATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                        <  1:1121595/71‑1 (MQ=255)                              ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                          >  1:2891181/1‑53 (MQ=255)                              ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                            <  1:693914/67‑1 (MQ=255)                              ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                        <  1:242644/71‑1 (MQ=255)                               gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                        >  1:36065/1‑67 (MQ=255)                               gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                        >  1:2613467/1‑67 (MQ=255)                                 caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                                                      <  1:1450100/54‑1 (MQ=255)                                 caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                     >  1:658597/1‑71 (MQ=255)                                 caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                     <  1:1369986/71‑1 (MQ=255)                                 caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                     <  1:3356469/71‑1 (MQ=255)                                   ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                            <  1:1035295/62‑1 (MQ=255)                                     acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                              >  1:850225/1‑42 (MQ=255)                                      cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                      >  1:2320488/1‑65 (MQ=255)                                      cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAtt                                    <  1:762724/67‑1 (MQ=255)                                      cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATAGACAGt                                       >  1:1564110/1‑64 (MQ=255)                                         cGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                             <  1:2998993/71‑1 (MQ=255)                                                gCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcg                      >  1:2261783/1‑71 (MQ=255)                                                 ccATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTcag                     <  1:923037/71‑3 (MQ=255)                                                    ttattaTGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGAc                  <  1:2480922/71‑1 (MQ=255)                                                        tatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                    <  1:910709/65‑1 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              >  1:2038524/1‑70 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              <  1:1986121/70‑1 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              <  1:1646315/70‑1 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              <  1:765074/70‑1 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              <  1:2999549/70‑1 (MQ=255)                                                         atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag              <  1:3080664/70‑1 (MQ=255)                                                          tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga             <  1:2338361/70‑1 (MQ=255)                                                          tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga             <  1:1195346/70‑1 (MQ=255)                                                          tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga             <  1:2166205/70‑1 (MQ=255)                                                             gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag         <  1:24227/71‑1 (MQ=255)                                                              tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag         <  1:2159079/70‑1 (MQ=255)                                                                ttGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc      >  1:1185062/1‑71 (MQ=255)                                                                 tGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc      >  1:3573672/1‑70 (MQ=255)                                                                 tGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc      >  1:3569103/1‑70 (MQ=255)                                                                  gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                     >  1:3881240/1‑54 (MQ=255)                                                                  gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc      <  1:1174304/69‑1 (MQ=255)                                                                  gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg     <  1:1485302/70‑1 (MQ=255)                                                                  gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg     <  1:681754/70‑1 (MQ=255)                                                                  gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATT‑‑CGACGTAGAAGAGCGCGAAc  >  1:2948877/1‑71 (MQ=255)                                                                   cAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                 >  1:3610725/1‑57 (MQ=255)                                                                   cAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGa    <  1:453335/70‑1 (MQ=255)                                                                    aGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc  >  1:2086589/1‑71 (MQ=255)                                                                    |                                                                       TTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAC  >  minE/1200993‑1201130

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