Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	798,949	C→T	M214I (ATG→ATA)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	798,949	0	C	T	100.0%	163.9 / NA	57	M214I (ATG→ATA)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/29);  total (28/29)

TCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA  >  minE/798881‑799015                                                                     |                                                                   tCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:2920064/1‑71 (MQ=255)  cTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  <  1:744440/70‑1 (MQ=255)  cTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:535859/1‑70 (MQ=255)   tttttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:1346239/1‑69 (MQ=255)   tttttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  <  1:1773269/69‑1 (MQ=255)     tttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATc                                                                 <  1:1143520/68‑1 (MQ=255)     tttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAgg                                                              >  1:2563085/1‑71 (MQ=255)          aTTTGCGGCAGGTATTGCTTGAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCaaa                                                         <  1:1779296/71‑1 (MQ=255)          aTTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCa                                                           >  1:689471/1‑69 (MQ=255)           tttGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:2872456/71‑1 (MQ=255)             tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:714812/1‑69 (MQ=255)             tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:2908026/1‑69 (MQ=255)               cGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAgg                                                    <  1:695571/71‑1 (MQ=255)               cGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAgg                                                    <  1:2862191/71‑1 (MQ=255)                      tATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaag                                                     <  1:28879/63‑1 (MQ=255)                       aTTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:913127/1‑59 (MQ=255)                       aTTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGc                                            >  1:1789109/1‑71 (MQ=255)                        ttGCTTCAATGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCg                                           <  1:2445896/71‑1 (MQ=255)                          gCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAg                                         >  1:1682802/1‑71 (MQ=255)                            ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaag                                                     <  1:691499/57‑1 (MQ=255)                            ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCg                                           <  1:1732338/67‑1 (MQ=255)                            ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCg                                       >  1:432470/1‑71 (MQ=255)                            ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCg                                       >  1:161083/1‑71 (MQ=255)                             tCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGtt                                              <  1:1596911/63‑1 (MQ=255)                              cAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGc                                     >  1:139052/1‑71 (MQ=255)                               aaCGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCt                                    >  1:1897303/1‑71 (MQ=255)                                aCGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:579679/57‑1 (MQ=255)                                  gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:233216/55‑1 (MQ=255)                                  gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:122578/55‑1 (MQ=255)                                  gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGCTTGCGAGCGGCTTa                                  >  1:662891/1‑70 (MQ=255)                                   cgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGTCTTaaa                                <  1:1020930/71‑1 (MQ=255)                                   cgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa                                <  1:2835425/71‑1 (MQ=255)                                    ggCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGAGCGGTTGCGAGCGGCTTaaa                                >  1:1371751/1‑70 (MQ=255)                                       aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTTAATAGt                            <  1:4293/71‑1 (MQ=255)                                       aTCCCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaa                                                      <  1:736753/45‑1 (MQ=255)                                        tCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:684828/42‑1 (MQ=255)                                         cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:208803/41‑1 (MQ=255)                                         cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:1514597/41‑1 (MQ=255)                                         cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          >  1:2216305/1‑71 (MQ=255)                                          gCCTTCATGCTGCGCCAGTAACTCTACTGTCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          <  1:2075664/70‑1 (MQ=255)                                          gCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAAc                         >  1:1998568/1‑71 (MQ=255)                                            cTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGtttc                                            >  1:2488433/1‑48 (MQ=255)                                               cATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                               >  1:319140/1‑44 (MQ=255)                                               cATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                               >  1:2305478/1‑44 (MQ=255)                                                   ctgcGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCgg                                      >  1:1721021/1‑49 (MQ=255)                                                   ctgcGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          >  1:1562831/1‑61 (MQ=255)                                                      cgcCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGt             <  1:2538757/71‑1 (MQ=255)                                                        ccAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATa                              >  1:2602138/1‑52 (MQ=255)                                                          aGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg         <  1:570269/71‑1 (MQ=255)                                                           gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg         >  1:1363671/1‑70 (MQ=255)                                                           gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:305344/71‑1 (MQ=255)                                                           gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:613395/71‑1 (MQ=255)                                                           gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:1904082/71‑1 (MQ=255)                                                            tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGc                <  1:2595292/62‑1 (MQ=255)                                                            tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa       >  1:1022819/1‑71 (MQ=255)                                                            tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa       >  1:895210/1‑71 (MQ=255)                                                                tctACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTgg          >  1:1412777/1‑64 (MQ=255)                                                                 ctactaTCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAAAAGTAACTGCGGTTGCCGTTGGGAATATGa  >  1:918833/1‑71 (MQ=255)                                                                     |                                                                   TCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA  >  minE/798881‑799015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑