Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,464,992	C→T	5.9%	A216T (GCC→ACC)	PP_2157 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,464,992	0	C	T	5.9%	203.9 / 2.1	68	A216T (GCC→ACC)	PP_2157	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGGCTGCGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGTGATTAGC  >  NC_002947/2464844‑2465140                                                                                                                                                     |                                                                                                                                                     gcacTGCATGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGc                                                                                                                                                      <  5:245748/146‑1 (MQ=255) aGGCTGCGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGc                                                                                                                                                      <  8:186683/149‑1 (MQ=255)  ggCTGCGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGcc                                                                                                                                                     <  5:374229/149‑1 (MQ=255)   gCTGCGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCa                                                                                                                                                    <  3:135443/149‑1 (MQ=255)       cGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTc                                                                                                                                                <  8:207375/149‑1 (MQ=255)           ggCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGc                                                                                                                                             <  2:393705/148‑1 (MQ=255)               ggggagggtgaCCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAAct                                                                                                                                        <  2:378596/146‑1 (MQ=255)                     ggtgaCCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCgg                                                                                                                                  >  1:361792/1‑149 (MQ=255)                       tgaCCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGTCAGCTCGGCCAACTCTTCggg                                                                                                                                 >  6:267883/1‑148 (MQ=255)                          ccTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATc                                                                                                                             <  6:4624/149‑1 (MQ=255)                                 gtgCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGcc                                                                                                                       <  8:345299/148‑1 (MQ=255)                                 gtgCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGcc                                                                                                                       >  7:345299/1‑148 (MQ=255)                                    cGTCAGCAGGTTGGACAGCGGCGTAGGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCg                                                                                                                   >  5:233806/1‑149 (MQ=255)                                     gTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCgg                                                                                                                  >  5:190416/1‑149 (MQ=255)                                        agcagGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCgg                                                                                                                  >  1:313698/1‑146 (MQ=255)                                            ggTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGTCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACgg                                                                                                           <  5:267883/149‑1 (MQ=255)                                                ggACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGggtg                                                                                                        <  2:361792/148‑1 (MQ=255)                                                 gACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGggtggt                                                                                                      <  5:71022/149‑1 (MQ=255)                                                 gACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGggtggt                                                                                                      >  1:114209/1‑149 (MQ=255)                                                    aGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCAc                                                                                                                     <  6:79060/131‑1 (MQ=255)                                                    aGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCAc                                                                                                                     >  5:79060/1‑131 (MQ=255)                                                    aGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAg                                                                                                   <  2:313698/149‑1 (MQ=255)                                                    aGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAg                                                                                                   >  4:222799/1‑149 (MQ=255)                                                       ggcgTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCTGGCATCTGCGCCAcc                                                                                                                    <  6:30420/129‑1 (MQ=255)                                                       ggcgTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATAGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCAcc                                                                                                                    >  5:30420/1‑129 (MQ=255)                                                       ggcgTACGGAACTCGTGGGCGTTTTCGGCCGTAAACGGTGACAGGCGCGGGAAGGCATCGTCCAGGCGCTGCAGGATCGCGTTCATGGTGCTGGCCAGCTGGGCCAAGTGTTCGGGCATCTGCGGGACGGGCAGACGGTTGGTGAGGGa                                                                                                <  4:111538/149‑1 (MQ=255)                                                         cgtacgtaACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAgc                                                                                              >  4:104828/1‑149 (MQ=255)                                                              gtaACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGTCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCt                                                                                         >  6:110441/1‑149 (MQ=255)                                                              gtaACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCt                                                                                         >  6:203485/1‑149 (MQ=255)                                                                   tCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAgg                                                                                                  >  1:193440/1‑135 (MQ=255)                                                                   tCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAgg                                                                                                  <  2:193440/135‑1 (MQ=255)                                                                   tCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGacac                                                                                    >  7:319181/1‑149 (MQ=255)                                                                   tCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGacac                                                                                    >  4:91966/1‑149 (MQ=255)                                                                          cGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTggcg                                                                             >  3:64678/1‑149 (MQ=255)                                                                               tCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCTCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCAc                                                                        <  4:249423/149‑1 (MQ=255)                                                                                 ggCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGTCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCAcc                                                                       <  5:110441/148‑1 (MQ=255)                                                                                 ggCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGacac                                                                                    <  6:100975/135‑1 (MQ=255)                                                                                 ggCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGacac                                                                                    >  5:100975/1‑135 (MQ=255)                                                                                 ggCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCt                                                                      >  6:25016/1‑149 (MQ=255)                                                                                    cGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGc                                                                   <  3:104828/149‑1 (MQ=255)                                                                                                   gcgcTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCAcc                                                                       <  1:396316/130‑1 (MQ=255)                                                                                                   gcgcTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCAcc                                                                       >  2:396316/1‑130 (MQ=255)                                                                                                     gcTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATggg                                                  <  5:203485/149‑1 (MQ=255)                                                                                                         gAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGgcg                                                >  8:134715/1‑147 (MQ=255)                                                                                                              cATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGc                                                                   <  4:357712/123‑1 (MQ=255)                                                                                                              cATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGc                                                                   >  3:357712/1‑123 (MQ=255)                                                                                                               aTCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGc                                        >  1:369377/1‑149 (MQ=255)                                                                                                                    ccAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgg                                   >  4:207195/1‑149 (MQ=255)                                                                                                                     cAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACTGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcggg                                  <  6:190416/149‑1 (MQ=255)                                                                                                                      aGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTCAGCGAGCGCGCTGCCACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgggc                                 >  7:261997/1‑149 (MQ=255)                                                                                                                       ggCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCa                                                       >  2:10701/1‑126 (MQ=255)                                                                                                                       ggCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgggcg                                <  5:215432/149‑1 (MQ=255)                                                                                                                       ggCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgggcg                                >  7:127093/1‑149 (MQ=255)                                                                                                                       ggCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCCGGCGGCCACCTGGCCCATCTGCCGCa                                                       <  1:10701/126‑1 (MQ=255)                                                                                                                        gcgcTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgg                                   >  2:139034/1‑145 (MQ=255)                                                                                                                        gcgcTGCAAGATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTgcgg                                   <  1:139034/145‑1 (MQ=255)                                                                                                                          gcTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGcc                             >  2:173116/1‑149 (MQ=255)                                                                                                                              cagcaTCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCa                           <  1:173116/147‑1 (MQ=255)                                                                                                                              cagcaTCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAgg                         <  4:64678/149‑1 (MQ=255)                                                                                                                              cagcaTCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCATCGGGCAGCTGCGCCACCGGCAGACGGGGGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCa                           <  3:222799/147‑1 (MQ=255)                                                                                                                                gcaTCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCa                       <  2:369377/149‑1 (MQ=255)                                                                                                                                     gcgTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCaa                  >  4:370939/1‑149 (MQ=255)                                                                                                                                      cgTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCaa                  >  1:321387/1‑148 (MQ=255)                                                                                                                                       gTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATa                <  2:114209/149‑1 (MQ=255)                                                                                                                                           aTGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTg             <  2:374647/148‑1 (MQ=255)                                                                                                                                           aTGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTg             <  8:319181/148‑1 (MQ=255)                                                                                                                                           aTGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGc            >  4:312281/1‑149 (MQ=255)                                                                                                                                              ctgctgGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGt         <  3:91966/149‑1 (MQ=255)                                                                                                                                               tgctgGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGt         >  1:42512/1‑148 (MQ=255)                                                                                                                                               tgctgGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGt         >  8:272720/1‑148 (MQ=255)                                                                                                                                                 ctgGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGTg        <  7:134715/147‑1 (MQ=255)                                                                                                                                                     ccAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGTGATTAgc  >  4:268098/1‑149 (MQ=255)                                                                                                                                                     ccAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGTGATTAgc  >  8:36677/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGCTGCGTGGGCGGGTGAGGGTGACCTGGGTGTGCGTCAGCAGGTTGGACAGCGGCGTACGTAACTCGTGGGCGATATCGGCCGAAAACGCTGACAGGCGCTGGAAGGCATCGTCCAGGCGCTGCAGCATCGCGTTCATGCTGCTGGCCAGCTCGGCCAACTCTTCGGGCATCTGCGCCACCGGCAGACGGGTGGTGAGGGAGCGCGCTGACACGCTGGCGGCCACCTGGCCCATCTGCCGCAATGGGCGCAGGCCGCTGCGGGCGGCCCAGGCACCCAATAGTGCGGTGATTAGC  >  NC_002947/2464844‑2465140

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: ‑