Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	406,027	A→C	8.7%	R10R (CGT→CGG)	PP_0337 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	406,027	0	A	C	8.7%	222.0 / 10.8	79	R10R (CGT→CGG)	PP_0337	sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base C (5/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGG  >  NC_002947/405880‑406172                                                                                                                                                      |                                                                                                                                                    acacaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTGGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCacac                                                                                                                                                    >  7:176709/1‑149 (MQ=255)  cacaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACg                                                                                                                                                   >  7:98507/1‑149 (MQ=255)        tCGCCGGTTTTCAGGCCGAACTGGCACTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTggc                                                                                                                                             >  4:109884/1‑149 (MQ=255)               ttttCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCAGCCGTGTggcggccgag                                                                                                                                      <  1:185101/149‑5 (MQ=255)               ttttCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCg                                                                                                                                      >  4:283471/1‑149 (MQ=255)                    aGGCCGAACTGGCGCTCACCGGCTGGGTCAAGGTAGCGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCTGCATCCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  1:276674/148‑1 (MQ=255)                    aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGATGATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTGACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  8:137104/148‑1 (MQ=255)                    aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  6:7493/148‑1 (MQ=255)                    aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  3:283471/148‑1 (MQ=255)                     ggCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTgg                                                                                                                                >  2:83204/1‑149 (MQ=255)                          aaCTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCtttt                                                                                                                           <  1:199637/149‑1 (MQ=255)                                 gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  1:125994/149‑1 (MQ=255)                                 gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  1:58133/149‑1 (MQ=255)                                 gcTCACAGGCTGGGTCAAGGTAGAGCATCGCCCAGGTGGCTTAGTTGAGACTTTCGAAACGCAGCATGCCGAGACGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACCCGGCTGGCGGCATCGGGTTCGCTTTTCATCGaa                                                                                                                    >  5:316284/1‑149 (MQ=255)                                   tCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑Ga                                                                                                                 >  4:228202/1‑149 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACt                                                                                                               >  3:149277/1‑149 (MQ=255)                                              gTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAag                                                                                                      <  8:98507/149‑1 (MQ=255)                                               tCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGc                                                                                                            <  7:233007/142‑1 (MQ=255)                                               tCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGc                                                                                                            >  8:233007/1‑142 (MQ=255)                                                 aaGGTAGAGCATGGCGCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAgg                                                                                                   <  3:5/149‑1 (MQ=255)                                                 aaGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAgg                                                                                                   <  4:234492/149‑1 (MQ=255)                                                     tAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑c                                                                                              <  1:83204/149‑1 (MQ=255)                                                       gagCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑cgc                                                                                            <  3:228202/149‑1 (MQ=255)                                                         gCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          >  8:32213/1‑149 (MQ=255)                                                         gCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          <  8:176709/149‑1 (MQ=255)                                                          cATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGc                                                                                         <  1:233660/149‑1 (MQ=255)                                                               cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGaa                                                                                                       >  6:333399/1‑131 (MQ=255)                                                               cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGaa                                                                                                       <  5:333399/131‑1 (MQ=255)                                                                ccAGGTGGGTTCGTTGAGCCTTTAGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCCCACGGCTGGCGGCATCGGGGTGGCTTTTAATCGAA‑GACTCGCTTGAAGAGGGCA‑CGCGGCCGtgtg                                                                                   >  2:110763/1‑149 (MQ=255)                                                                  aGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          >  6:36449/1‑140 (MQ=255)                                                                  aGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          <  5:36449/140‑1 (MQ=255)                                                                       gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    <  2:286565/141‑1 (MQ=255)                                                                       gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    >  1:286565/1‑141 (MQ=255)                                                                       gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTgcgg                                                                            <  7:32213/149‑1 (MQ=255)                                                                            gTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGAt                                                                       >  3:205179/1‑149 (MQ=255)                                                                              tGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    >  2:235004/1‑134 (MQ=255)                                                                              tGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    <  1:235004/134‑1 (MQ=255)                                                                                aGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    <  1:42663/148‑1 (MQ=255)                                                                                 gCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCAGGGGTAACTGTGTCACG‑‑ACCTCGGCCGCGACACTGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  <  2:93692/149‑1 (MQ=255)                                                                                  ccTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  >  2:253161/1‑148 (MQ=255)                                                                                   cTTTCGAAACGCAGCATGCCGAGCCGCGCGGGCACGTGTAACTGTGTCACGTCACATCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  <  1:110763/149‑1 (MQ=255)                                                                                    tttCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAg                                                               >  7:152422/1‑149 (MQ=255)                                                                                    tttCGAAACGCAGCATGCCGAGCCGCGAGGGCACGCGTAAGTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCGTTTAATGGAA‑GACGATCTTGAAAAGGGAA‑CGCTGCATTGTGCCTGAGGAGGATTTGCCAAg                                                               >  4:189327/1‑149 (MQ=255)                                                                                      tCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGt                                                             >  2:250470/1‑149 (MQ=255)                                                                                        gAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTg                                                           >  5:222882/1‑149 (MQ=255)                                                                                            cGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCatg                                                       <  1:253161/149‑1 (MQ=255)                                                                                                  aTGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTcc                                                 <  3:109884/149‑1 (MQ=255)                                                                                                    gCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAACGACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    <  1:245518/129‑1 (MQ=255)                                                                                                    gCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    >  2:245518/1‑128 (MQ=255)                                                                                                      cGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAACCGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              <  1:250470/149‑1 (MQ=255)                                                                                                      cGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                             >  8:315580/1‑149 (MQ=255)                                                                                                       gAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              >  3:285552/1‑147 (MQ=255)                                                                                                       gAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              <  4:285552/147‑1 (MQ=255)                                                                                                           cgcgAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTg                                        >  5:253721/1‑149 (MQ=255)                                                                                                                  gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  6:253721/149‑1 (MQ=255)                                                                                                                  gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  4:205179/149‑1 (MQ=255)                                                                                                                  gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  3:166770/149‑1 (MQ=255)                                                                                                                   cACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCgg                                >  1:161468/1‑149 (MQ=255)                                                                                                                    aCGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGc                               <  8:285092/149‑1 (MQ=255)                                                                                                                           aCTGTGTCACG‑‑ACCTCTGCCGCCACACGGCTGGCGGCATCGGGTTGTCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGTTTTGGCAAGGTTGCATTATGTCCCGAGGGCTGGCAAGCGGCCATGCt                         <  6:316284/148‑1 (MQ=255)                                                                                                                           aCTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        >  7:269990/1‑149 (MQ=255)                                                                                                                               tgtCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                    >  4:246680/1‑149 (MQ=255)                                                                                                                                 tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGTCGGATTTGGCAAGGTTGCATCATGTAACGAGGGCTGGCAGGCGGCCCTGGTTgcttct                  >  8:223417/1‑144 (MQ=255)                                                                                                                                 tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                  >  2:147347/1‑149 (MQ=255)                                                                                                                                 tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                  <  8:152422/149‑1 (MQ=255)                                                                                                                                  ccgctctccTCGGCCGCCACACGGCTGGCGGGATCGGGTTGGCTGTTCATCGAA‑GACTCGCTTGCAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGTCAAGGGTGCATGATGGCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                    <  7:241907/141‑1 (MQ=255)                                                                                                                                   aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        <  1:209937/141‑1 (MQ=255)                                                                                                                                   aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        >  2:209937/1‑141 (MQ=255)                                                                                                                                   aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                >  3:160547/1‑149 (MQ=255)                                                                                                                                   aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                >  7:72172/1‑149 (MQ=255)                                                                                                                                    cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  7:315580/149‑1 (MQ=255)                                                                                                                                    cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  6:222882/149‑1 (MQ=255)                                                                                                                                    cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  3:246680/149‑1 (MQ=255)                                                                                                                                         ccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa             <  5:58155/143‑1 (MQ=255)                                                                                                                                         ccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa             >  6:58155/6‑148 (MQ=255)                                                                                                                                          cTCGGCCGCCACACGGCTGGCGGCAGCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGATGGAA‑CGCGGCCTTGTGGCTGAGGCGGCTTTGGCAAGGTTGCATGATGGCCCGAGGGCTGGCACGAGGCCCTGCTAGCTGcagctcgg            >  8:241907/1‑141 (MQ=255)                                                                                                                                           gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GCCTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg          >  3:291221/4‑149 (MQ=255)                                                                                                                                           gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg          >  3:108598/4‑149 (MQ=255)                                                                                                                                            atgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAgg         >  3:192489/3‑149 (MQ=255)                                                                                                                                                  ccACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCg   >  7:196789/1‑149 (MQ=255)                                                                                                                                                   cacaCGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCgg  <  1:147347/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGG  >  NC_002947/405880‑406172

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