Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,991	A→T	8.8%	I166F (ATC→TTC)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,991	0	A	T	8.8%	249.2 / 12.6	91	I166F (ATC→TTC)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/38);  new base T (4/4);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACC  >  NC_002947/5187845‑5188138                                                                                                                                                   |                                                                                                                                                             gCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGTCTGCTGGACTTCGAGGAAGGCGTGGCGCAACACGCCTTc                                                                                                                                                            <  7:284829/149‑1 (MQ=255)     cgcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                                          <  4:145506/147‑1 (MQ=255)     cgcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                                          >  3:145506/1‑147 (MQ=255)       ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCg                                                                                                                                                        >  3:283051/1‑147 (MQ=255)       ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                                      <  7:364715/149‑1 (MQ=255)       ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                                      <  5:62134/149‑1 (MQ=255)       ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                                      <  1:421154/149‑1 (MQ=255)        cTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCg                                                                                                                                                        >  3:6810/1‑146 (MQ=255)           tgAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAagtc                                                                                                                                                   <  5:9368/148‑3 (MQ=38)              aaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTTCTACCTGCTCACCCAGCGTGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTcc                                                                                                                                               <  5:172623/149‑1 (MQ=255)               aCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCg                                                                                                                                              <  8:347371/149‑1 (MQ=255)               aCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCg                                                                                                                                              <  6:252292/149‑1 (MQ=255)                ccGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGc                                                                                                                                             >  1:102210/1‑149 (MQ=255)                        gCAGTGAAGATGTTGGCGGTCGAGGTCACGCCGCAACTCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGATGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGcc                                                                                                                             >  1:48495/1‑148 (MQ=255)                           gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑Ag‑g                                                                                                                                >  1:399207/1‑142 (MQ=255)                           gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑Ag‑g                                                                                                                                <  2:399207/142‑1 (MQ=255)                           gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTACTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAagtcc                                                                                                                                                  >  5:226179/1‑130 (MQ=37)                           gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTACTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAagtcc                                                                                                                                                  <  6:226179/133‑4 (MQ=37)                              aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCtt                                                                                                                      >  1:316099/1‑149 (MQ=255)                               aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                     >  5:270510/1‑149 (MQ=255)                               aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                     <  1:240321/149‑1 (MQ=255)                               aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                     >  5:63292/1‑149 (MQ=255)                                  tGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGc                                                                                                                  <  3:224916/149‑1 (MQ=255)                                   gTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGc                                                                                                                  >  8:119576/1‑148 (MQ=255)                                    ttGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCa                                                                                                                <  4:283051/149‑1 (MQ=255)                                     tGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAg                                                                                                               >  1:103286/1‑149 (MQ=255)                                      ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAg                                                                                                               <  1:247318/148‑1 (MQ=255)                                             gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCgcgg                                                                                                       <  6:341501/149‑1 (MQ=255)                                              aGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                      >  1:55915/1‑149 (MQ=255)                                                 tCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                                   <  1:106623/149‑1 (MQ=255)                                                  cccGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                                   <  2:366645/148‑1 (MQ=255)                                                  cccGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                                   >  1:366645/1‑148 (MQ=255)                                                  cccGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGc                                                                                                  <  6:290463/149‑1 (MQ=255)                                                   ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCt                                                                                                 >  5:404844/1‑149 (MQ=255)                                                      ccgcAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGc                                                                                              <  7:107834/149‑1 (MQ=255)                                                       cgcAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGc                                                                                              >  3:17736/1‑148 (MQ=255)                                                            cgctcgctCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAaca                                                                                        <  2:102210/149‑1 (MQ=255)                                                                cgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAAGCCCGCCGCGCGCCTggcc                                                                                            <  6:344624/141‑4 (MQ=255)                                                                cgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAAGCCCGCCGCGCGCCTggcc                                                                                            >  5:344624/1‑138 (MQ=255)                                                                 gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                    >  8:103066/1‑148 (MQ=255)                                                                 gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                    <  8:25818/148‑1 (MQ=255)                                                                 gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAc                                                                                   <  4:23264/149‑1 (MQ=255)                                                                     gATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagca                                                                               <  2:316099/149‑1 (MQ=255)                                                                        aCCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                            >  8:239826/1‑149 (MQ=255)                                                                         cccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGTCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                            >  6:334814/1‑148 (MQ=255)                                                                          ccTGCTCTACTACCTGCTCACCCCGCGCGAAGCCGGCCTGCTGCACTTCGAGGAAGGCGTGATGCGCCAAGCCATCCCCCAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGccc                                                                           >  8:134989/1‑148 (MQ=255)                                                                          ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGCACAGCAGccc                                                                           >  8:175108/1‑148 (MQ=255)                                                                          ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTAGGCCAGGCGCGCAGCGGGCTTGCGAAACACGAACAGCAGccc                                                                           >  7:364493/1‑148 (MQ=255)                                                                           cTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTa                                                                         >  2:232905/1‑149 (MQ=255)                                                                                 tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGa                                                                   >  7:347651/1‑149 (MQ=255)                                                                                    tacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                                  >  4:373389/1‑147 (MQ=255)                                                                                      ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCTGGTGGTAGATGAGCGGGCAACGCATCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTAc                                                                        >  1:359613/1‑148 (MQ=255)                                                                                          cTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCa                                                          >  5:166428/1‑149 (MQ=255)                                                                                           tCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCAt                                                         >  1:252026/1‑149 (MQ=255)                                                                                             aCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGc                                                       >  4:424126/1‑149 (MQ=255)                                                                                              cccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                      <  7:119576/149‑1 (MQ=255)                                                                                               ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCGAAATCAGCCCCATGcc                                                      >  1:150933/1‑148 (MQ=255)                                                                                               ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                      >  1:231747/1‑148 (MQ=255)                                                                                               ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                      >  6:239710/1‑148 (MQ=255)                                                                                                cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                      >  7:189735/1‑147 (MQ=255)                                                                                                cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                      >  3:412921/1‑147 (MQ=255)                                                                                                cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCca                                                    <  7:236902/149‑1 (MQ=255)                                                                                                 aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                                   <  2:103286/149‑1 (MQ=255)                                                                                                   cgcgCAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagca                                                 <  2:350528/149‑1 (MQ=255)                                                                                                     cgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGTCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                               <  5:334814/149‑1 (MQ=255)                                                                                                     cgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                               >  8:170124/1‑149 (MQ=255)                                                                                                     cgAAGACGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                               <  4:17736/149‑1 (MQ=255)                                                                                                             gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCgg                                       >  5:393118/1‑149 (MQ=255)                                                                                                             gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCgg                                       >  6:430596/1‑149 (MQ=255)                                                                                                              ccTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                      <  1:137814/149‑1 (MQ=255)                                                                                                                  ctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGc                                    >  8:64803/1‑147 (MQ=255)                                                                                                                    ggACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAgg                                >  1:286897/1‑149 (MQ=255)                                                                                                                     gACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAgg                                <  3:424126/148‑1 (MQ=255)                                                                                                                     gACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGc                               <  6:404844/149‑1 (MQ=255)                                                                                                                      aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                              <  6:270510/149‑1 (MQ=255)                                                                                                                      aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                              <  7:103066/149‑1 (MQ=255)                                                                                                                      aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                              <  2:55915/149‑1 (MQ=255)                                                                                                                       cTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGcc                                   <  4:203133/143‑1 (MQ=38)                                                                                                                       cTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGcc                                   >  3:203133/1‑143 (MQ=38)                                                                                                                        ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                            >  7:59245/1‑149 (MQ=38)                                                                                                                        ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                            >  8:224811/1‑149 (MQ=38)                                                                                                                        ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                            >  5:129325/1‑149 (MQ=255)                                                                                                                            aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                        <  6:393118/149‑1 (MQ=255)                                                                                                                            aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                        >  1:24480/1‑149 (MQ=255)                                                                                                                             ggAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                        >  6:144736/1‑148 (MQ=255)                                                                                                                              gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCg                       >  3:386555/1‑148 (MQ=255)                                                                                                                              gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑CGCCTACGCTACGGCCAGGCGCGCGGCTGGCTTCCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCACCAGTCTGAGCGGTGCGAGGCGGTTCCCg                       >  4:312773/1‑148 (MQ=255)                                                                                                                               aaGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa                     <  5:144736/149‑1 (MQ=255)                                                                                                                                 ggCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCg                       >  6:405515/1‑145 (MQ=255)                                                                                                                                   cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATc                 <  2:231747/149‑1 (MQ=255)                                                                                                                                         gcgcCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCaa            <  6:166428/148‑1 (MQ=255)                                                                                                                                             cacgccaTCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCa       >  2:117845/1‑149 (MQ=255)                                                                                                                                              acgccaTCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCAt      <  4:386555/149‑1 (MQ=255)                                                                                                                                               cgccaTCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATc     <  5:239710/149‑1 (MQ=255)                                                                                                                                                  caTCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCAcc  <  2:359613/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                             GCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACC  >  NC_002947/5187845‑5188138

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