Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	226,159	C→A	P163T (CCT→ACT)	gmhB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	226,159	0	C	A	100.0%	181.1 / NA	54	P163T (CCT→ACT)	gmhB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/31);  total (23/31)

GCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGA  >  REL606/226111‑226205                                                 |                                               gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACt                                              <  1:1108340/51‑1 (MQ=255) gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACt                                              <  1:6415888/51‑1 (MQ=255)   gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGa                                            <  1:328493/51‑1 (MQ=255)    aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGaa                                           <  1:5799142/51‑1 (MQ=255)     aCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                          >  1:7282673/1‑51 (MQ=255)      cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGc                                         <  1:4920655/51‑1 (MQ=255)      cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGc                                         <  1:3625115/51‑1 (MQ=255)       gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCa                                        <  1:4365542/51‑1 (MQ=255)        tGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAg                                       >  1:3853927/1‑51 (MQ=255)         gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGa                                      <  1:4092680/51‑1 (MQ=255)          ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                     <  1:561351/51‑1 (MQ=255)          ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                     >  1:4201095/1‑51 (MQ=255)            aaCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaaaa                                   <  1:7254487/51‑1 (MQ=255)             aCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAc                                  <  1:1939852/51‑1 (MQ=255)             aCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAc                                  <  1:6748187/51‑1 (MQ=255)               aaaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgc                                >  1:896699/1‑51 (MQ=255)                aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                               <  1:6773983/51‑1 (MQ=255)                aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                               >  1:7058187/1‑51 (MQ=255)                aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                               >  1:5235902/1‑51 (MQ=255)                   gTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                            <  1:5027922/51‑1 (MQ=255)                   gTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                            <  1:3616034/51‑1 (MQ=255)                    tGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGcggcgg                           <  1:2209481/51‑1 (MQ=255)                    tGCTGGTGCGTACGGGCAAACCTATTACGACTGAAGCAGAAAACGcggcgg                           >  1:4016315/1‑51 (MQ=255)                       tGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAtt                        >  1:2175915/1‑51 (MQ=255)                       tGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAtt                        <  1:2381747/51‑1 (MQ=255)                        ggTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTg                       <  1:3690026/51‑1 (MQ=255)                         gTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTgg                      >  1:3043958/1‑51 (MQ=255)                          tGCGTACGGGAAAACCTATTACGACTGAAGCAGAAAACGcggcgg                           >  1:796270/1‑45 (MQ=38)                             gTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTggg                     >  1:1469603/1‑48 (MQ=255)                              tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                 >  1:7171823/1‑51 (MQ=255)                              tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                 <  1:3621892/51‑1 (MQ=255)                              tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                 <  1:1893633/51‑1 (MQ=255)                              tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                 >  1:2680817/1‑51 (MQ=255)                              tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                 >  1:3006609/1‑51 (MQ=255)                               aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTggg                     >  1:5022179/1‑46 (MQ=255)                               aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTa                >  1:6507423/1‑51 (MQ=255)                               aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTa                >  1:2837205/1‑51 (MQ=255)                                 gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa              >  1:5723656/1‑51 (MQ=255)                                 gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa              <  1:7033026/51‑1 (MQ=255)                                   gTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATa            <  1:4993723/51‑1 (MQ=255)                                    tAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAg           <  1:119403/51‑1 (MQ=255)                                      aaCCTATTACGACTGAAGCAGAAAACgcgg                              >  1:1301711/1‑30 (MQ=38)                                      aaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAAt             >  1:6829306/1‑47 (MQ=255)                                       aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt        <  1:671791/51‑1 (MQ=255)                                       aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt        <  1:5624080/51‑1 (MQ=255)                                        ccTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTg       <  1:5883/51‑1 (MQ=255)                                           aTTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa    <  1:5191097/51‑1 (MQ=255)                                           aTTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa    <  1:6060314/51‑1 (MQ=255)                                            ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg   <  1:6316799/51‑1 (MQ=255)                                            ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg   >  1:2906367/1‑51 (MQ=255)                                            ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg   <  1:4316866/51‑1 (MQ=255)                                            ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg   >  1:3349489/1‑51 (MQ=255)                                             tACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGa  <  1:6690845/51‑1 (MQ=255)                                             tACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGa  >  1:1951912/1‑51 (MQ=255)                                                 |                                               GCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGA  >  REL606/226111‑226205

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