Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,972	A→T	7.4%	E159D (GAA→GAT)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,972	0	A	T	7.4%	273.6 / 9.6	95	E159D (GAA→GAT)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base T (4/3);  total (52/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 = 5.66e-01

CGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAA  >  NC_002947/5187825‑5188119                                                                                                                                                    |                                                                                                                                                             cgcAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAg                                                                                                                                                             <  2:48309/149‑1 (MQ=255)  gcagcaTTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTTCTACCTGCTCACCCAGCGTGAAGCCGGCCTGCTGGACTTCGAGGAAgg                                                                                                                                                            >  6:172623/1‑149 (MQ=255)      caTTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCgtg                                                                                                                                                        >  4:224916/1‑149 (MQ=255)           cTGGATTTCGGCTGCGCCTGTGAAACCGAACCAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgc                                                                                                                                                    <  1:325059/148‑1 (MQ=255)            tGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgcgc                                                                                                                                                  >  2:240321/1‑149 (MQ=255)             ggACTTCGGCTGCGCCTGCGAAACCGAACCAGCAGTGAAGATGTTGGCGGTCGAGGGCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGcgcc                                                                                                                                                 <  4:20887/149‑1 (MQ=255)              gATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgcgc                                                                                                                                                  <  4:97300/147‑1 (MQ=255)                tttCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGcgccacg                                                                                                                                              >  2:421154/1‑149 (MQ=255)                     gCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGTCTGCTGGACTTCGAGGAAGGCGTGGCGCAACACGCCTTc                                                                                                                                         <  7:284829/149‑1 (MQ=255)                         cgcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                       >  3:145506/1‑147 (MQ=255)                         cgcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                       <  4:145506/147‑1 (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                                                                                                                           <  6:252292/149‑1 (MQ=255)                                   aCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCg                                                                                                                           <  8:347371/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                                                                                >  1:366645/1‑148 (MQ=255)                                                                      cccGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                <  2:366645/148‑1 (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                                   >  6:239710/1‑148 (MQ=255)                                                                                                                   ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                   >  1:231747/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                    >  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           <  2:55915/149‑1 (MQ=255)                                                                                                                                          aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg           <  6:270510/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         >  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)                                                                                                                                                    |                                                                                                                                                             CGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAA  >  NC_002947/5187825‑5188119

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