Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,995	C→T	7.4%	P167L (CCC→CTC)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,995	0	C	T	7.4%	276.7 / 7.5	94	P167L (CCC→CTC)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/40);  new base T (4/3);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

CGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCG  >  NC_002947/5187849‑5188139                                                                                                                                                   |                                                                                                                                                          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                                                                                                                                 <  2:399207/142‑1 (MQ=255)                       gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑Ag‑g                                                                                                                                 >  1:399207/1‑142 (MQ=255)                       gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTACTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAagtcc                                                                                                                                                   <  6:226179/133‑4 (MQ=37)                       gTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTACTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAagtcc                                                                                                                                                   >  5:226179/1‑130 (MQ=37)                          aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCtt                                                                                                                       >  1:316099/1‑149 (MQ=255)                           aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                      >  5:63292/1‑149 (MQ=255)                           aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                      <  1:240321/149‑1 (MQ=255)                           aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                                      >  5:270510/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                                                                                             >  5:344624/1‑138 (MQ=255)                                                            cgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAAGCCCGCCGCGCGCCTggcc                                                                                             <  6:344624/141‑4 (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                                                       >  3:412921/1‑147 (MQ=255)                                                                                            cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                       >  7:189735/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                                        >  6:430596/1‑149 (MQ=255)                                                                                                         gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCgg                                        >  5:393118/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                               <  2:55915/149‑1 (MQ=255)                                                                                                                  aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                               <  7:103066/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                             >  8:224811/1‑149 (MQ=38)                                                                                                                    ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                             >  7:59245/1‑149 (MQ=38)                                                                                                                    ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                             >  5:129325/1‑149 (MQ=255)                                                                                                                        aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                         >  1:24480/1‑149 (MQ=255)                                                                                                                        aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                         <  6:393118/149‑1 (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)                                                                                                                                                tCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCg  >  6:297726/1‑148 (MQ=255)                                                                                                                                                tCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCg  >  2:434175/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                          CGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCG  >  NC_002947/5187849‑5188139

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