Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,426,626	Δ1 bp	7.1%	coding (268/942 nt)	PP_2127 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,426,626	0	C	.	7.1%	222.1 / 9.8	70	coding (268/942 nt)	PP_2127	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base . (2/3);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

GCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCA  >  NC_002947/2426485‑2426769                                                                                                                                               |                                                                                                                                                gCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCggg                                                                                                                                          >  4:163009/1‑149 (MQ=255)  cAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCgg                                                                                                                                           <  2:189275/147‑1 (MQ=255)        gtcgGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCCGAGCT‑GGCCGAGCGCat                                                                                                                                     <  1:25376/147‑2 (MQ=255)            ggCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGAt                                                                                                                               <  3:278329/149‑1 (MQ=255)              cATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATc                                                                                                                              >  5:337541/1‑148 (MQ=255)                  cgcCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCt                                                                                                                         <  4:71322/149‑1 (MQ=255)                  cgcCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCt                                                                                                                         <  3:163009/149‑1 (MQ=255)                         gCAGCCACTGCAGTTGCAGCTTGGCAGCTGTGAGCATTTGACCGATCTCGTCTTGCTGTTCCCGGCCAATATGCTGGCGCTCGTCTTCTTGCACCTTGATCATGCGCTCGTCC‑ATCTCGGCCGGTCTCAGGTGTATCGACCTGTCgcc                                                                                                                   <  6:198779/148‑1 (MQ=255)                             ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                              >  5:42782/1‑149 (MQ=255)                             ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                              <  2:367125/149‑1 (MQ=255)                             ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGAGCTCGTCTTCTTGCACCTAGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                              >  7:149955/1‑149 (MQ=255)                              cACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                             >  1:194770/1‑149 (MQ=255)                              cACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                             >  3:68572/1‑149 (MQ=255)                                cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggccca                                                                                                           <  2:378388/149‑1 (MQ=255)                                cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                             >  2:107515/1‑147 (MQ=255)                                 tGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcccag                                                                                                          >  5:257582/1‑149 (MQ=255)                                   cAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGtt                                                                                                        <  3:20887/149‑1 (MQ=255)                                    aGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTg                                                                                                       >  7:322712/1‑149 (MQ=255)                                        gCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTAGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTac                                                                                                    >  1:24631/1‑148 (MQ=255)                                         cAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCAGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaa                                                                                                  >  4:295267/1‑149 (MQ=255)                                           gCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaaca                                                                                                <  7:164486/149‑1 (MQ=255)                                           gCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaaca                                                                                                <  1:107515/149‑1 (MQ=255)                                                  aGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCa                                                                                         <  7:309798/149‑1 (MQ=255)                                                    cGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGc                                                                                       >  1:260243/1‑149 (MQ=255)                                                      gTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGt                                                                                     >  4:104465/1‑149 (MQ=255)                                                       tGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTg                                                                                    <  5:136156/149‑1 (MQ=255)                                                         aGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                                  >  3:359714/1‑149 (MQ=255)                                                           cAGTTGCCCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGa                                                                                <  8:346873/149‑1 (MQ=255)                                                           cAGGTGACCGATGGCGTCGTGCAGTTCCCGGCCCAGATGCGGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGa                                                                                <  4:27256/149‑1 (MQ=255)                                                            aGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                               >  8:256664/1‑149 (MQ=255)                                                            aGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                               <  4:306264/149‑1 (MQ=255)                                                              ttGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAcc                                                                              >  4:14147/1‑148 (MQ=255)                                                                gACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCaa                                                                            <  6:257582/148‑1 (MQ=255)                                                                    gATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTCCAACCGCACCCCGCGGGGCGCCCAACTgcc                                                                       >  6:258495/1‑148 (MQ=255)                                                                    gATGTCGTCGTGCAGTTCCCCGCCCAGATGCTGGCGCTCGTCTTCTTTCACC‑TGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTgcag                                                                      <  6:249418/149‑1 (MQ=255)                                                                     aTGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTgcag                                                                      <  8:322712/149‑1 (MQ=255)                                                                        tcgtcgTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCAGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCAc                                                                   <  3:295267/149‑1 (MQ=255)                                                                           tcgtGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCCGCACcag                                                                >  2:189611/1‑149 (MQ=255)                                                                                      ccGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGcc                                                      <  4:359714/148‑1 (MQ=255)                                                                                       cGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCATCGTGGCGACCAACTGCAGCACCAGCAGTCTGGCCgg                                                    <  6:326532/149‑1 (MQ=255)                                                                                       cGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCgg                                                    <  7:256664/149‑1 (MQ=255)                                                                                            cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTg                                               >  8:122311/1‑149 (MQ=255)                                                                                                     gcgcTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                        <  7:41733/147‑1 (MQ=255)                                                                                                     gcgcTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTg                                       <  6:257028/148‑1 (MQ=255)                                                                                                           gTCTTCTTGCACCTTGCGCATGCGCGCGGCC‑CGCTCGGCCGGGCGCAGGGGGCTCGACCTGGCGCCCAGGCCCAGTTGTACACCAGCACCCCGCGTGGCGACCAACTGCAGCACCCGCAGGCTGGCCGGGATTGCCTGGTTGGCCAgg                                 <  5:258495/148‑1 (MQ=255)                                                                                                           gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAgg                                 >  6:280233/1‑148 (MQ=255)                                                                                                           gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                  <  2:260243/147‑1 (MQ=255)                                                                                                           gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                  <  3:14147/147‑1 (MQ=255)                                                                                                                 ttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgcg                          >  3:133543/1‑149 (MQ=255)                                                                                                                   gCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgcg                          <  6:334970/147‑1 (MQ=255)                                                                                                                   gCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCCGGATTGCCTGGTTGGCCAGGTACAgcg                          >  5:334970/1‑147 (MQ=255)                                                                                                                    cACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGcc                        <  6:337541/148‑1 (MQ=255)                                                                                                                     aCCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCAt                      >  6:312435/1‑149 (MQ=255)                                                                                                                      ccTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCa                                   <  5:6359/135‑1 (MQ=255)                                                                                                                      ccTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCa                                   >  6:6359/1‑135 (MQ=255)                                                                                                                       cTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGt                    >  2:321490/1‑149 (MQ=255)                                                                                                                         tccacctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTg                                           <  8:94005/118‑1 (MQ=38)                                                                                                                          gAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGtt                   <  6:42782/147‑1 (MQ=255)                                                                                                                            acctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGTAGGCTGGCCGGGa                                                  <  4:87398/111‑1 (MQ=38)                                                                                                                            acctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGTAGGCTGGCCGGGa                                                  >  3:87398/4‑114 (MQ=38)                                                                                                                              ctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAAt             >  7:271365/2‑149 (MQ=255)                                                                                                                              aTGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGcc                        >  8:157681/1‑138 (MQ=255)                                                                                                                                gcgcTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAAt             >  6:106742/1‑147 (MQ=255)                                                                                                                                   cTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGt         >  2:238715/1‑148 (MQ=255)                                                                                                                                   cTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGt         >  4:20804/1‑148 (MQ=255)                                                                                                                                      ggCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTgg                 >  3:22696/1‑137 (MQ=255)                                                                                                                                      ggCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTgg                 <  4:22696/137‑1 (MQ=255)                                                                                                                                      ggCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGa       <  1:238715/147‑1 (MQ=255)                                                                                                                                      ggCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGa       <  1:321490/147‑1 (MQ=255)                                                                                                                                         c‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCa  >  3:33493/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                GCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCA  >  NC_002947/2426485‑2426769

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