Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,981,882	T→C	6.0%	W843R (TGG→CGG)	PP_5242 →	GAF domain‑containing sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,981,882	0	T	C	6.0%	200.9 / 2.6	67	W843R (TGG→CGG)	PP_5242	GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CAGCGTGCTCGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTC  >  NC_002947/5981735‑5982028                                                                                                                                                    |                                                                                                                                                     cAGCGTGCTCGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTg                                                                                                                                                     <  8:159286/149‑1 (MQ=255)     gTGCTCGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGgcag                                                                                                                                                 <  4:257751/149‑1 (MQ=255)      tGCTCGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGgcag                                                                                                                                                 >  7:278237/1‑148 (MQ=255)       gctcgctcGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGgcag                                                                                                                                                 >  3:306872/1‑147 (MQ=255)          cgctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGt                                                                                                                                            >  8:58287/1‑149 (MQ=255)           gctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTc                                                                                                                                           <  8:49249/149‑1 (MQ=255)           gctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTc                                                                                                                                           >  8:209552/1‑149 (MQ=255)               ggtggtGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACCTGCTGCCAGTCGAgg                                                                                                                                       <  6:327133/149‑1 (MQ=255)               ggTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAgg                                                                                                                                       >  1:121107/1‑149 (MQ=255)                gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGc                                                                                                                                      >  4:218653/1‑149 (MQ=255)                gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACCTGCTGCCAGTCGAGGc                                                                                                                                      <  2:184209/149‑1 (MQ=255)                gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGCTAGCAGCCCGGTGCGGGGCCTCGACTCGCAGCAGGTCGAGGc                                                                                                                                      <  8:119891/149‑1 (MQ=255)                gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGAACAACTTGCGCGGTGACTATCAGACGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGCGCCTCGACTGGCACCAGGTCGAGGc                                                                                                                                      >  3:152180/1‑149 (MQ=255)                  gtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGc                                                                                                                                      >  5:146339/1‑147 (MQ=255)                      gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAg                                                                                                                                        >  6:65310/1‑141 (MQ=255)                      gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAg                                                                                                                                        <  5:65310/141‑1 (MQ=255)                      gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCGTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCa                                                                                                                                >  1:61611/1‑149 (MQ=255)                        ccAGGCACGGGCGTGCCGGCTGCAGGCGATGCTGTACACCTTGCGCGGTGACTATCAGACGGGCATCGACAGCGCGCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATg                                                                                                                              <  4:152180/149‑1 (MQ=255)                         cAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGc                                                                                                                             <  1:137144/149‑1 (MQ=255)                           ggCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGcc                                                                                                                            >  3:147974/1‑148 (MQ=255)                                   cGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑g                                                                                                                 <  4:131074/149‑1 (MQ=255)                                           cTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGAt                                                                                                         <  3:218653/149‑1 (MQ=255)                                            agCAGGCCATGCTGTACACCATGCGCCTAGACTATCAGGCGGCCATCGACAGCGCCCTGTCCGGCCTCGACTTACTGGGTATCAGCCCGGTGCGGGCCCTCGACTGCCACCAGGTCTAGGCCAGTCATGCA‑CATGTTCA‑GGCGTTGCTc                                                                                                        >  8:154997/2‑149 (MQ=255)                                              cAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATc                                                                                                        <  6:146339/147‑1 (MQ=255)                                               aGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGaa                                                                                                     <  4:188217/149‑1 (MQ=255)                                                ggCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAAc                                                                                                    <  2:264078/149‑1 (MQ=255)                                                 gcgcTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACa                                                                                                   <  4:147974/149‑1 (MQ=255)                                                        tgTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGAc                                                                                            >  5:58909/1‑149 (MQ=255)                                                         gTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGAcc                                                                                           >  8:255868/1‑149 (MQ=255)                                                          tACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCg                                                                                          >  5:18976/1‑149 (MQ=255)                                                          tACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCg                                                                                          <  6:58909/149‑1 (MQ=255)                                                            cacCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCgg                                                                                         <  7:209552/148‑1 (MQ=255)                                                             acCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGc                                                                                                               <  1:307795/125‑1 (MQ=255)                                                             acCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGc                                                                                                               >  2:307795/1‑125 (MQ=255)                                                               cTTGCGCGGTGACTATCAGGCTGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCAtt                                                                                     <  7:154997/149‑1 (MQ=255)                                                               cTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCAtt                                                                                     >  4:93006/1‑149 (MQ=255)                                                                    gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                >  8:195573/1‑149 (MQ=255)                                                                    gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                <  7:197/149‑1 (MQ=255)                                                                    gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                >  8:10529/1‑149 (MQ=255)                                                                      ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                              <  6:18976/149‑1 (MQ=255)                                                                      ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                              <  7:255868/149‑1 (MQ=255)                                                                      ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                              <  8:278237/149‑1 (MQ=255)                                                                                      ccATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                              <  2:47584/133‑1 (MQ=255)                                                                                      ccATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                              >  1:47584/1‑133 (MQ=255)                                                                                         tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGTGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGa                                                            >  6:83596/1‑148 (MQ=255)                                                                                         tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGAc                                                           <  7:195573/149‑1 (MQ=255)                                                                                            aCAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCg                                                        >  2:66684/1‑149 (MQ=255)                                                                                              aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCgga                                                      <  2:121107/149‑1 (MQ=255)                                                                                              aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCgga                                                      <  4:207956/149‑1 (MQ=255)                                                                                               gcgcCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCggag                                                     <  5:43494/149‑1 (MQ=255)                                                                                                         cggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACgg                                           >  7:8066/1‑149 (MQ=255)                                                                                                           gccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                         <  2:61611/149‑1 (MQ=255)                                                                                                             cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                       >  1:119112/1‑149 (MQ=255)                                                                                                             cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                       <  2:119112/149‑1 (MQ=255)                                                                                                                 aCTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGTGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATc                                   <  5:83596/149‑1 (MQ=255)                                                                                                                     aCTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGa                                                                                                          >  6:260980/1‑74 (MQ=255)                                                                                                                     aCTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGa                                                                                                          <  5:260980/74‑1 (MQ=255)                                                                                                                       tGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCAGGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCct                              >  1:253297/1‑149 (MQ=255)                                                                                                                         ggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCGCATGTGAA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGACACTGCCGCGCCCCGAAACGGAGGAAGTGACGTTTGCCATCACCctgc                            >  1:82700/1‑149 (MQ=255)                                                                                                                                 cccGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCAc                     >  6:162041/1‑147 (MQ=255)                                                                                                                                 cccGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCAc                     <  5:162041/147‑1 (MQ=255)                                                                                                                                     gTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGtc               <  3:93006/149‑1 (MQ=255)                                                                                                                                     gTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGt                >  1:247381/1‑148 (MQ=255)                                                                                                                                         ggggCCTCGACTGGGAGGAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGATGGTCAGGCATTGCCTGGATTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGtcatca           <  2:111032/149‑1 (MQ=255)                                                                                                                                             ccTCGACCTGCTGCCAGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCct       >  7:100915/1‑149 (MQ=255)                                                                                                                                             ccTCGACCTGCTGCCAGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCct       >  8:43287/1‑149 (MQ=255)                                                                                                                                                  aCTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCttcttc  <  4:306872/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CAGCGTGCTCGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCC‑CATGTGCA‑GGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTC  >  NC_002947/5981735‑5982028

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