Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260480	2,308,233	(C)5→4	intergenic (+176/‑161)	h16_B2050 → / → h16_B2051	response regulator containing a GGDEF effector domain/conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	2,308,229	0	C	.	100.0%	152.0 / NA	34	intergenic (+172/‑165)	h16_B2050/h16_B2051	response regulator containing a GGDEF effector domain/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (16/18);  total (16/18)

CGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGCC  >  AM260480/2308100‑2308355                                                                                                                                  |                                                                                                                               cggcCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccg                                                                                                                            <  1:600266/134‑2 (MQ=255)  ggcCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccgc                                                                                                                           >  1:947521/1‑132 (MQ=255)    ccACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccgcga                                                                                                                         >  1:739109/1‑130 (MQ=255)            cAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCg                                                                                                                <  2:739109/134‑1 (MQ=255)            cAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCg                                                                                                                <  1:543183/134‑1 (MQ=255)             aaCCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGc                                                                                                               <  1:781560/134‑1 (MQ=255)                      accaGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGc                                                                                                               <  2:879669/125‑1 (MQ=255)                      accaGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGc                                                                                                               >  1:879669/1‑125 (MQ=255)                       ccaGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTg                                                                                                     <  2:62599/134‑1 (MQ=255)                                            cAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGa                                                                                <  1:903839/134‑1 (MQ=255)                                              ggCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGa                                                                                >  2:903839/1‑132 (MQ=255)                                                   cGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTg                                                                          >  2:994295/1‑133 (MQ=255)                                                           cGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTg                                                                 <  2:962306/134‑1 (MQ=255)                                                              gCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTgcgc                                                              >  1:953089/1‑134 (MQ=255)                                                               ccTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTggg                                                                                                        <  1:897453/91‑1 (MQ=255)                                                               ccTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTggg                                                                                                        >  2:897453/1‑91 (MQ=255)                                                                      ggccCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGcaac                                                       >  2:865786/1‑133 (MQ=255)                                                                        ccccGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCaacaac                                                    <  1:865786/134‑1 (MQ=255)                                                                            gCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGc                                                                                                               <  1:743973/71‑1 (MQ=255)                                                                            gCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGc                                                                                                               >  2:743973/1‑71 (MQ=255)                                                                                 gcacgcacgcGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGcc                                            >  2:919771/1‑133 (MQ=255)                                                                                       acgcGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGc                                     <  1:994295/134‑1 (MQ=255)                                                                                          cgTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCt                                                                                                                      >  1:830245/1‑50 (MQ=39)                                                                                          cgTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCt                                                                                                                      <  2:830245/50‑1 (MQ=39)                                                                                               ccGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTa                                                                                              >  2:440982/1‑69 (MQ=255)                                                                                               ccGCTTGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTa                                                                                              <  1:440982/69‑1 (MQ=255)                                                                                                   ttGTCGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGgc                         >  2:156924/1‑134 (MQ=255)                                                                                                       cGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATc                                                                                    <  2:1063475/71‑1 (MQ=255)                                                                                                       cGGAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATc                                                                                    >  1:1063475/1‑71 (MQ=255)                                                                                                         gAATTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGccc                   >  1:796522/1‑134 (MQ=255)                                                                                                           aTTGTCCGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGccccc                 <  2:796522/134‑1 (MQ=255)                                                                                                                 cGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCgg           >  2:587364/1‑134 (MQ=255)                                                                                                                 cGACATGCGATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCgg           >  1:355875/1‑134 (MQ=255)                                                                                                                         gATTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGc   <  2:355875/134‑1 (MQ=255)                                                                                                                          aTTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGcc  <  2:947521/134‑1 (MQ=255)                                                                                                                          aTTCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGcc  <  1:156924/134‑1 (MQ=255)                                                                                                                           ttCGGAG‑CCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGc   >  2:503690/1‑132 (MQ=255)                                                                                                                                  |                                                                                                                               CGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGCC  >  AM260480/2308100‑2308355

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