Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,968	A→G	7.5%	E158G (GAG→GGG)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,968	0	A	G	7.5%	259.8 / 9.1	93	E158G (GAG→GGG)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/42);  new base G (4/3);  total (48/45)
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.74e-01

TACGCGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACC  >  NC_002947/5187821‑5188116                                                                                                                                                    |                                                                                                                                                              tACGCGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAg                                                                                                                                                              <  5:104302/149‑1 (MQ=255) tACGCGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAg                                                                                                                                                              <  3:390238/149‑1 (MQ=255)   cgcgCTGCTTTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCTCGAAGCCGGCCTGCTGGACTTCGAGGa                                                                                                                                                            <  3:354155/149‑1 (MQ=255)     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                                                                                                                                <  1:421154/149‑1 (MQ=255)                               ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                <  5:62134/149‑1 (MQ=255)                               ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                <  7:364715/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: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                                >  6:239710/1‑148 (MQ=255)                                                                                                                       ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                >  1:231747/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                 >  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        <  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  >  1:24480/1‑149 (MQ=255)                                                                                                                                                    aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc  <  6:393118/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                              TACGCGCAGCATTGCTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACC  >  NC_002947/5187821‑5188116

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