Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,984	C→A	9.7%	R163R (CGC→CGA) ‡	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,984	0	C	A	9.7%	270.1 / 13.7	90	R163R (CGC→CGA) ‡	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (4/4);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAA  >  NC_002947/5187836‑5188128                                                                                                                                                     |                                                                                                                                                          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                                                                                                                      >  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                                                                                                           <  1:240321/149‑1 (MQ=255)                                        aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTc                                                                                                           >  5:270510/1‑149 (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:25818/148‑1 (MQ=255)                                                                          gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                          >  8:103066/1‑148 (MQ=255)                                                                          gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAc                                                                         <  4:23264/149‑1 (MQ=255)                                                                              gATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagca                                                                     <  2:316099/149‑1 (MQ=255)                                                                                 aCCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                  >  8:239826/1‑149 (MQ=255)                                                                                  cccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGTCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                  >  6:334814/1‑148 (MQ=255)                                                                                   ccTGCTCTACTACCTGCTCACCCCGCGCGAAGCCGGCCTGCTGCACTTCGAGGAAGGCGTGATGCGCCAAGCCATCCCCCAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGccc                                                                 >  8:134989/1‑148 (MQ=255)                                                                                   ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGCACAGCAGccc                                                                 >  8:175108/1‑148 (MQ=255)                                                                                   ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTAGGCCAGGCGCGCAGCGGGCTTGCGAAACACGAACAGCAGccc                                                                 >  7:364493/1‑148 (MQ=255)                                                                                    cTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTa                                                               >  2:232905/1‑149 (MQ=255)                                                                                          tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGa                                                         >  7:347651/1‑149 (MQ=255)                                                                                             tacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                        >  4:373389/1‑147 (MQ=255)                                                                                               ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCTGGTGGTAGATGAGCGGGCAACGCATCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTAc                                                              >  1:359613/1‑148 (MQ=255)                                                                                                   cTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCa                                                >  5:166428/1‑149 (MQ=255)                                                                                                    tCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCAt                                               >  1:252026/1‑149 (MQ=255)                                                                                                      aCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGc                                             >  4:424126/1‑149 (MQ=255)                                                                                                       cccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                            <  7:119576/149‑1 (MQ=255)                                                                                                        ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCGAAATCAGCCCCATGcc                                            >  1:150933/1‑148 (MQ=255)                                                                                                        ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                            >  1:231747/1‑148 (MQ=255)                                                                                                        ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                            >  6:239710/1‑148 (MQ=255)                                                                                                         cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                            >  7:189735/1‑147 (MQ=255)                                                                                                         cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                            >  3:412921/1‑147 (MQ=255)                                                                                                         cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCca                                          <  7:236902/149‑1 (MQ=255)                                                                                                          aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                         <  2:103286/149‑1 (MQ=255)                                                                                                            cgcgCAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagca                                       <  2:350528/149‑1 (MQ=255)                                                                                                              cgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGTCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                     <  5:334814/149‑1 (MQ=255)                                                                                                              cgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                     >  8:170124/1‑149 (MQ=255)                                                                                                              cgAAGACGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAgg                                     <  4:17736/149‑1 (MQ=255)                                                                                                                      gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCgg                             >  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                    <  7:103066/149‑1 (MQ=255)                                                                                                                               aCTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                    <  2:55915/149‑1 (MQ=255)                                                                                                                                cTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGcc                         >  3:203133/1‑143 (MQ=38)                                                                                                                                cTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGcc                         <  4:203133/143‑1 (MQ=38)                                                                                                                                 ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                  >  7:59245/1‑149 (MQ=38)                                                                                                                                 ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                  >  8:224811/1‑149 (MQ=38)                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                  >  5:129325/1‑149 (MQ=255)                                                                                                                                     aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc              <  6:393118/149‑1 (MQ=255)                                                                                                                                     aGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc              >  1:24480/1‑149 (MQ=255)                                                                                                                                      ggAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc              >  6:144736/1‑148 (MQ=255)                                                                                                                                       gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCg             >  3:386555/1‑148 (MQ=255)                                                                                                                                       gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑CGCCTACGCTACGGCCAGGCGCGCGGCTGGCTTCCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCACCAGTCTGAGCGGTGCGAGGCGGTTCCCg             >  4:312773/1‑148 (MQ=255)                                                                                                                                        aaGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa           <  5:144736/149‑1 (MQ=255)                                                                                                                                          ggCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCg             >  6:405515/1‑145 (MQ=255)                                                                                                                                            cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATc       <  2:231747/149‑1 (MQ=255)                                                                                                                                                  gcgcCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCaa  <  6:166428/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                          TGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCT‑‑‑‑‑‑‑‑AG‑GGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAA  >  NC_002947/5187836‑5188128

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