Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	782,109	T→C	Q167R (CAG→CGG)	hlyE ←	hemolysin E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	782,109	0	T	C	100.0%	179.5 / NA	52	Q167R (CAG→CGG)	hlyE	hemolysin E
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/27);  total (25/27)

TGGACCGGCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGTTTCCCGG  >  minE/782042‑782175                                                                    |                                                                   tGGACCGGCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCgg                                                                  >  1:1394068/1‑70 (MQ=255) tGGACCGGCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCgg                                                                  >  1:2575497/1‑70 (MQ=255) tGGACCGGCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGa                                                                 <  1:2027761/71‑1 (MQ=255)        gCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGaaa                                                               <  1:2747051/66‑1 (MQ=255)        gCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAgct                                                          <  1:2129490/71‑1 (MQ=255)             gacACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCttt                                                     <  1:2629204/71‑1 (MQ=255)              acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                                                    <  1:2252259/71‑1 (MQ=255)                 ccGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTc                                                 >  1:1201772/1‑71 (MQ=255)                 ccGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTc                                                 <  1:329942/71‑1 (MQ=255)                  cGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGaaa                                                               <  1:1679789/56‑1 (MQ=255)                    gCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttttt                                                  >  1:56152/1‑67 (MQ=255)                     cTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCTGAAATAGCTGCTTTTTTCTGa                                              <  1:2614921/70‑1 (MQ=255)                     cTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGaa                                             <  1:2700055/71‑1 (MQ=255)                      tGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCt                                                >  1:2926008/1‑67 (MQ=255)                        cGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAgctgc                                                        <  1:2649087/57‑1 (MQ=255)                         ggCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAAt                                         >  1:661080/1‑71 (MQ=255)                          gCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGaaa                                                               <  1:2058327/48‑1 (MQ=255)                           cACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCa                                       >  1:2371674/1‑71 (MQ=255)                            aCCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGaaa                                            <  1:1533907/65‑1 (MQ=255)                             ccGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAAt                                                              <  1:919909/46‑1 (MQ=255)                             ccGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGaaaaa                                          <  1:2740209/66‑1 (MQ=255)                             ccGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCAtt                                     >  1:2241636/1‑71 (MQ=255)                              cGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGaa                                                                >  1:294010/1‑43 (MQ=255)                                gCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGt                                  >  1:948736/1‑71 (MQ=255)                                 cATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGtt                                 <  1:2501971/71‑1 (MQ=255)                                   tatGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTaa                               <  1:1753112/71‑1 (MQ=255)                                     tGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATa                                                             >  1:1585921/1‑39 (MQ=255)                                     tGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACt                             <  1:1516048/71‑1 (MQ=255)                                      gCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTg                            >  1:1301394/1‑71 (MQ=255)                                       cttccttcCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCAt                                      <  1:1396689/60‑1 (MQ=255)                                       cttccttcCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTgg                           <  1:676083/71‑1 (MQ=255)                                       cttccttcCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTg                            <  1:2191966/70‑1 (MQ=255)                                         tccttcctGATTTTATCTACCTGTGTCCGGAAATAgctg                                                         >  1:389586/1‑39 (MQ=38)                                          ccttcctGATTTTATCTACCTGTGACCGGAAATAgctg                                                         >  1:959058/1‑38 (MQ=255)                                          ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc                          <  1:3031502/69‑1 (MQ=255)                                           cttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAAc                              >  1:331834/1‑64 (MQ=255)                                            ttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTa                                <  1:2345859/61‑1 (MQ=255)                                            ttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTa                                <  1:2961109/61‑1 (MQ=255)                                            ttcctGATTTTATCTACCTGTCACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTa                                <  1:1970854/61‑1 (MQ=255)                                               ctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTa                    <  1:142543/70‑1 (MQ=255)                                                  aTTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGc                <  1:192117/71‑1 (MQ=255)                                                  aTTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGc                >  1:2514297/1‑71 (MQ=255)                                                    tttATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGc                >  1:2641499/1‑69 (MQ=255)                                                    tttATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGc                >  1:1808911/1‑69 (MQ=255)                                                         cTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGt         >  1:259105/1‑71 (MQ=255)                                                         cTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGt         >  1:1489399/1‑71 (MQ=255)                                                           aCCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCa                                       >  1:2745309/1‑39 (MQ=255)                                                             cTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGTTTcc     >  1:1731778/1‑71 (MQ=255)                                                              tgtgACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATc                      >  1:2695840/1‑53 (MQ=255)                                                                tgACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGTTTCCCgg  <  1:1042780/71‑1 (MQ=255)                                                                 gACCGGAAATAGCTGCTTTTTTCTGAAAAATCAtt                                     >  1:2231194/1‑35 (MQ=255)                                                                 gACCGGAAATAGCTGCTTTTTTCTGAAAAATCAtt                                     >  1:752477/1‑35 (MQ=255)                                                                    |                                                                   TGGACCGGCGACGACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAGCAGTTTCCCGG  >  minE/782042‑782175

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