Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	406,027	A→C	5.6%	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	5.6%	281.7 / 6.5	90	R10R (CGT→CGG)	PP_0337	sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (43/42);  new base C (3/2);  total (46/44)
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.85e-01

ACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCT  >  NC_002947/405880‑406174                                                                                                                                                    |                                                                                                                                                    acacaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGTATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCacac                                                                                                                                                    <  8:205616/149‑1 (MQ=255)  cacaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACg                                                                                                                                                   >  7:305414/1‑149 (MQ=255)   acaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACgg                                                                                                                                                  >  2:288034/1‑149 (MQ=255)        tCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTggc                                                                                                                                             >  5:114016/1‑149 (MQ=255)          gCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCAGCCGTGTggcgg                                                                                                                                           >  8:230533/1‑149 (MQ=255)          gCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTggcgg                                                                                                                                           >  8:45585/1‑149 (MQ=255)              gTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATc                                                                                                                                       <  3:180603/149‑1 (MQ=255)              gTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATc                                                                                                                                       >  3:285048/1‑149 (MQ=255)               ttttCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCAGCCGTGTggcggccgag                                                                                                                                      <  7:230533/149‑5 (MQ=255)                  tCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGt                                                                                                                                   >  1:172121/1‑149 (MQ=255)                      gCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCa                                                                                                                                         <  5:269356/139‑1 (MQ=255)                      gCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCa                                                                                                                                         >  6:269356/1‑139 (MQ=255)                      gCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGc                                                                                                                               >  5:2563/1‑149 (MQ=255)                      gCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGc                                                                                                                               >  6:319719/1‑149 (MQ=255)                          aaCTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCtttt                                                                                                                           <  1:288034/149‑1 (MQ=255)                            cTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCa                                                                                                                         <  1:220156/149‑1 (MQ=255)                            cTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCa                                                                                                                         <  5:289486/149‑1 (MQ=255)                               gcgcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATc                                                                                                                       >  8:220870/1‑148 (MQ=255)                                 gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  6:2563/149‑1 (MQ=255)                                 gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  6:114016/149‑1 (MQ=255)                                  cTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCg                                                                                                                      >  3:51040/1‑146 (MQ=255)                                  cTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  4:51040/148‑1 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGt                                                                                                                                   <  4:231008/130‑1 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGt                                                                                                                                   >  3:231008/1‑130 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATc                                                                                                                       <  6:58903/142‑1 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATc                                                                                                                       >  5:58903/1‑142 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGa                                                                                                                  <  2:226657/147‑1 (MQ=255)                                     acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACt                                                                                                                >  3:236546/1‑149 (MQ=255)                                          cTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCAt                                                                                                                        <  5:297426/136‑1 (MQ=255)                                          cTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCAt                                                                                                                        >  6:297426/1‑136 (MQ=255)                                           tGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCtt                                                                                                           >  4:106409/1‑148 (MQ=255)                                             ggTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTg                                                                                                          <  8:305414/147‑1 (MQ=255)                                             ggTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGaa                                                                                                        <  6:302094/149‑1 (MQ=255)                                             ggTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGaa                                                                                                        <  5:319719/149‑1 (MQ=255)                                               tCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAaga                                                                                                      <  3:106409/149‑1 (MQ=255)                                                 aaGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAgg                                                                                                    >  1:93163/1‑149 (MQ=255)                                                   ggTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCAGCCGTGTggcgg                                                                                                                                           >  6:258788/1‑108 (MQ=255)                                                    gTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCa                                                                                                                         <  3:177044/125‑1 (MQ=255)                                                    gTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCa                                                                                                                         >  4:177044/1‑125 (MQ=255)                                                              gCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTtg                                                                                       >  2:309362/1‑149 (MQ=255)                                                               cccAGGTGGGTTCGTTGAGCCTTTCGAAACGCAGCTTGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTtgt                                                                                      <  2:172121/149‑1 (MQ=255)                                                               cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCtttt                                                                                                                           >  2:276494/1‑112 (MQ=255)                                                               cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCtttt                                                                                                                           <  1:276494/112‑1 (MQ=255)                                                               cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTtgt                                                                                      <  5:327038/149‑1 (MQ=255)                                                                     tGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTgc                                                                                <  2:93163/149‑1 (MQ=255)                                                                       gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTgcgg                                                                              <  1:309362/149‑1 (MQ=255)                                                                         ttCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    >  5:16570/1‑109 (MQ=255)                                                                         ttCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  6:16570/109‑1 (MQ=255)                                                                               gAGCCTTTCGAAACGCAGCATGACGAGCCGCGAGGGAACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTCGGCTTTTCATCGAAGACTCGCTTGCAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTg                                                                      >  4:230634/1‑149 (MQ=255)                                                                                aGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTg                                                                      <  7:45585/148‑1 (MQ=255)                                                                                 gCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTgg                                                                     >  1:137699/1‑148 (MQ=255)                                                                                 gCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTAATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTgg                                                                     >  6:67773/1‑148 (MQ=255)                                                                                    tttCGAAACGCAGCATCCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGCGAACGCGGCCTTGACGCTGCGGCGGATTTGCCAAg                                                                 >  8:310598/1‑149 (MQ=255)                                                                                         aaaCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTg                                                                      >  4:296311/1‑139 (MQ=255)                                                                                         aaaCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTg                                                                      <  3:296311/139‑1 (MQ=255)                                                                                         aaaCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTg                                                             >  2:297574/1‑148 (MQ=255)                                                                                           aCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCat                                                          <  4:285048/149‑1 (MQ=255)                                                                                           aCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCat                                                          <  4:236546/149‑1 (MQ=255)                                                                                            cGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCat                                                          >  2:349291/1‑148 (MQ=255)                                                                                               agcaTGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCatgatg                                                      >  7:99497/1‑149 (MQ=255)                                                                                                   tGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCATACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTccc                                                  <  5:321735/149‑1 (MQ=255)                                                                                                     ccGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTccc                                                  >  5:629/1‑147 (MQ=255)                                                                                                     ccGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                                <  8:99497/149‑1 (MQ=255)                                                                                                      cGAGCCGCGAGGGCACGGGTAACTGTGTCACGCCCTCGGCCGCCACACGGCTGGCGTCATCGTGTTGGCTTTTCATCGAATACTCGCTTGAATATGGAACGCGGCCTTGTGGCTGCGGCGGATTTCGCAAGGTTGCATGATGTCCCGa                                                <  7:220870/148‑1 (MQ=255)                                                                                                          cctcgAGGGCACGGGTAACTGGGTCACGACCTCGGCCGCCTCACGGCTGGCGGCTTCGGGTTGGCGTTTCATCGAAGACTCGCTTGAAGACGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGTGGGcg                                           <  2:243936/146‑2 (MQ=255)                                                                                                          ccGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCt                                           >  7:312330/1‑149 (MQ=255)                                                                                                           cgcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTg                                          <  5:67773/149‑1 (MQ=255)                                                                                                           cgcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTg                                          <  8:312330/149‑1 (MQ=255)                                                                                                            gcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTgg                                         >  6:341614/1‑149 (MQ=255)                                                                                                             cgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTgg                                         >  4:145115/1‑148 (MQ=255)                                                                                                              gAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGc                                        >  5:277113/1‑148 (MQ=255)                                                                                                                      gggTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTgg                                                                     >  3:148704/1‑111 (MQ=255)                                                                                                                      gggTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTgg                                                                     <  4:148704/111‑1 (MQ=255)                                                                                                                        gTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCGTCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATg                             >  1:1040/1‑149 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTg                                                             <  5:247471/108‑1 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTg                                                             >  6:247471/1‑108 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAggg                                             >  6:138454/1‑124 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAggg                                             <  5:138454/124‑1 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                    >  3:333920/1‑149 (MQ=255)                                                                                                                                 tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                    >  7:203822/1‑149 (MQ=255)                                                                                                                                  cACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGa                                                                          >  1:277405/1‑94 (MQ=255)                                                                                                                                  cACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGa                                                                          <  2:277405/94‑1 (MQ=255)                                                                                                                                    cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg                 <  3:41186/149‑1 (MQ=255)                                                                                                                                    cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg                 <  1:297574/149‑1 (MQ=255)                                                                                                                                        cTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa               >  3:260066/1‑147 (MQ=255)                                                                                                                                         gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                      >  3:344865/4‑139 (MQ=255)                                                                                                                                         gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                      <  4:344865/136‑1 (MQ=255)                                                                                                                                           ggCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg                 <  2:46337/142‑1 (MQ=255)                                                                                                                                           ggCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg                 >  1:46337/1‑142 (MQ=255)                                                                                                                                           ggCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg            <  6:277113/147‑1 (MQ=255)                                                                                                                                           ggCCGCCACACGGCTGGCGGCATCGGGTCGGCTTTTCATCGAAGACTCGCTTGAAGAGTGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGTTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg            <  3:230634/147‑1 (MQ=255)                                                                                                                                                  acacGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGc   <  1:349291/149‑1 (MQ=255)                                                                                                                                                   cacGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCt  <  4:333920/149‑1 (MQ=255)                                                                                                                                                    acGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAggg                                             >  2:263994/1‑105 (MQ=255)                                                                                                                                                    acGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCt  >  7:92028/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCT  >  NC_002947/405880‑406174

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