Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,561,943	C→T	E380K (GAA→AAA)	digH ←	glycoside hydrolase family 10 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,561,943	0	C	T	100.0%	146.1 / NA	38	E380K (GAA→AAA)	digH	glycoside hydrolase family 10 protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/22);  total (16/22)

TGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGC  >  NZ_CP009273/1561807‑1562072                                                                                                                                         |                                                                                                                                  tGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 >  1:301031/1‑139 (MQ=255) tGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 >  2:746241/1‑139 (MQ=255)   cGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                <  1:1244339/138‑1 (MQ=255)   cGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                <  1:740520/138‑1 (MQ=255)   cGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                >  2:801582/1‑138 (MQ=255)   cGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                <  1:381838/138‑1 (MQ=255)   cGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 <  1:1000968/137‑1 (MQ=255)        tATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACtt                                                                                                                           <  2:916883/138‑1 (MQ=255)            cAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTtata                                                                                                                       <  2:805021/138‑1 (MQ=255)             aGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGa                                                                                                                     <  2:1279792/139‑1 (MQ=255)                tAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAgg                                                                                                                  >  2:968626/1‑139 (MQ=255)                   tCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGa                                                                                                               <  2:1112051/139‑1 (MQ=255)                   tCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGa                                                                                                               <  2:592939/139‑1 (MQ=255)                     cTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 >  2:441090/1‑119 (MQ=255)                     cTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 <  1:441090/119‑1 (MQ=255)                        aCGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATAcc                                                                                                           >  1:866569/1‑138 (MQ=255)                             aCAAGATGGTGCCGCTAATTTCGGGCACATCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGAtata                                                                                                     <  2:301031/139‑1 (MQ=255)                             aCAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAc                                                                                                                                 <  2:357975/111‑1 (MQ=255)                             aCAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGAtata                                                                                                     <  2:472123/139‑1 (MQ=255)                                 gATGGTGCCGGTAATTTCGGGCACAGCATCGTTAAGATCGGGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAAACCACCAGTCTGGCTCTAACTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAgg                                                                                                 <  2:682476/139‑1 (MQ=255)                                 gATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAgg                                                                                                 >  1:20272/1‑139 (MQ=255)                                  aTGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGc                                                                                                <  1:629495/139‑1 (MQ=255)                                      tGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGt                                                                                            >  2:1245354/1‑139 (MQ=255)                                        ccGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTc                                                                                           <  2:465818/138‑1 (MQ=255)                                              aTTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGt                                                                                      <  1:801582/137‑1 (MQ=255)                                                               gTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGcc                                                                   <  1:504904/139‑1 (MQ=255)                                                                ttAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                >  2:252317/1‑77 (MQ=255)                                                                ttAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTAcc                                                                                                                                <  1:252317/77‑1 (MQ=255)                                                                        gAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGACCTGGTGGGTTTATAGACATTCGCCCCCCAttc                                                           >  1:494601/1‑137 (MQ=255)                                                                          gctgctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGaa                                                                                                                    <  2:1040211/79‑1 (MQ=255)                                                                          gctgctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGaa                                                                                                                    >  1:1040211/1‑79 (MQ=255)                                                                              ctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAac                                                    <  1:748881/139‑1 (MQ=255)                                                                              ctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAac                                                    <  1:1245354/139‑1 (MQ=255)                                                                                            gTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTGCTGGGCGGTTTCAAGCCCTCCACCCCCCCATTTTTCACCACGTCATAAagccc                                      >  2:649940/1‑137 (MQ=255)                                                                                                          atcatcCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGGCGGGTTAAAGAAATACGCCCAACATTTTTACCACACGTCATAAAGGGCAGGAAGCCGTgag                         >  1:269589/1‑138 (MQ=255)                                                                                                              tcCAGTCTGGCTCTATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCCTGTCGGTTTAACGACATCCCCCCCCCATTTTTCCAAAACCTCATAAAGCGCCGCACTCCGTcagcaacg                    >  1:34628/1‑134 (MQ=255)                                                                                                                        ctctATCTTTGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTGCTGGTCGGGTTTAAGAAATCCGCCCACCATTTTTCTCACAACTCATAACACGCGGGTATCCCTGAGAAACGCCAGCGAAtt          >  2:677355/1‑139 (MQ=255)                                                                                                                                 tGAAGGTTTACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCTGGGCGGGTTTACGACATACGCCCCCCATTTTGCCAACACGCCATAAAGCCCGGCACTCCGCGAGAACGGCCCGGAAATTTTGGGAgg  >  1:368470/1‑137 (MQ=255)                                                                                                                                         |                                                                                                                                  TGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGC  >  NZ_CP009273/1561807‑1562072

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