Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,573,161	C→G	7.3%	P85R (CCC→CGC)	PP_4054 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,573,161	0	C	G	7.3%	170.9 / 4.0	55	P85R (CCC→CGC)	PP_4054	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/29);  new base G (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

CACTGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAGCACAATTCCCCTTCCATC  >  NC_002947/4573013‑4573308                                                                                                                                                      |                                                                                                                                                      cACTGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGcc                                                                                                                                                       >  7:170773/1‑149 (MQ=255) cACTGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGcc                                                                                                                                                       >  6:3755/1‑149 (MQ=255)  aCTGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGccc                                                                                                                                                      <  8:339660/149‑1 (MQ=255)    tGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGACCGCCGCCAAGGCGCAAGCGCCAGCCGAACCCAAGGCCGCAGCCCAGCCCGGGGGTAAACCCAAGGCCGTGA‑AAAAGCCCCg                                                                                                                                                    >  1:295357/1‑149 (MQ=255)      ggAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCggg                                                                                                                                                  <  6:35723/149‑1 (MQ=255)          aTGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGt                                                                                                                                              <  6:354548/149‑1 (MQ=255)               ccGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGaaa                                                                                                                                          >  6:75140/1‑148 (MQ=255)                     gCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTa                                                                                                                                   <  6:84643/149‑1 (MQ=255)                      cTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTAc                                                                                                                                  >  3:226182/1‑149 (MQ=255)                      cTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTAc                                                                                                                                  >  7:83111/1‑149 (MQ=255)                           cgaggccgaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCGCCGGGCGGTGAAAAAGCCTACTGCCg                                                                                                                             >  7:135980/1‑149 (MQ=255)                           cgaggccgaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGGTAAACCCACTGCGGTGA‑AAAAGCGCCGGGCGGTGAAAAAGCCTACTGCCg                                                                                                                             >  7:141103/1‑149 (MQ=255)                             aggccgaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGt                                                                                                                           <  1:322317/149‑1 (MQ=255)                              ggccgaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGtt                                                                                                                          <  1:1379/149‑1 (MQ=255)                                 cgaAGCGCTGGTGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCAACGGGCGGTGAAAAAGCCTACTGCCGGTTAAg                                                                                                                       >  4:115055/1‑149 (MQ=255)                                 cgaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAg                                                                                                                       >  3:332053/1‑149 (MQ=255)                                  gaAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGc                                                                                                                      <  6:335151/149‑1 (MQ=255)                                   aaGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGcc                                                                                                                     <  7:179266/149‑1 (MQ=255)                                   aaGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTCCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGcc                                                                                                                     >  4:231115/1‑149 (MQ=255)                                    aGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCGACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCt                                                                                                                    <  6:150515/149‑1 (MQ=255)                                     gcgcTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTg                                                                                                                   <  5:3755/149‑1 (MQ=255)                                             gcCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGcc                                                                                                         >  6:106560/1‑149 (MQ=255)                                                      aGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCgg                                                                                                <  4:332053/149‑1 (MQ=255)                                                      aGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCgg                                                                                                <  5:75140/149‑1 (MQ=255)                                                        ccccGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCggg                                                                                               >  2:106156/1‑148 (MQ=255)                                                            gcgcAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCGCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAAc                                                                                          <  8:135980/149‑1 (MQ=255)                                                            gcgcAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCGCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAAc                                                                                          <  8:141103/149‑1 (MQ=255)                                                            gcgcAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAAc                                                                                          >  5:315866/1‑149 (MQ=255)                                                              gcAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAAccc                                                                                        <  4:226182/149‑1 (MQ=255)                                                                  gaagGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTc                                                                                    <  4:775/149‑1 (MQ=255)                                                                          cGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGt                                                                            <  6:315866/149‑1 (MQ=255)                                                                             ccAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGtca                                                                         >  4:310276/1‑149 (MQ=255)                                                                                 gCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGc                                                                                      >  3:185387/1‑132 (MQ=255)                                                                                 gCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGc                                                                                      <  4:185387/132‑1 (MQ=255)                                                                                     aaGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGgc                                                                 <  6:322635/149‑1 (MQ=255)                                                                                           ccTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGt                                                           <  3:115055/149‑1 (MQ=255)                                                                                            cTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATg                                                            >  3:153411/1‑147 (MQ=255)                                                                                            cTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTa                                                          <  5:106560/149‑1 (MQ=255)                                                                                                   aCCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGTCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGTCGCATGTAGGAGCgg                                                   <  8:170773/149‑1 (MQ=255)                                                                                                     ccAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGt                                                 <  4:153411/149‑1 (MQ=255)                                                                                                      cAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCCTGGGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCggg                                                  <  1:106156/149‑1 (MQ=255)                                                                                                       aaGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTTAACCCGCTGCTACAGGGTGATGACAGGGGGATGTAGGGGCGGGgtt                                               >  2:134894/1‑149 (MQ=255)                                                                                                       aaGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                  >  7:339516/1‑114 (MQ=255)                                                                                                       aaGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                  <  8:339516/114‑1 (MQ=255)                                                                                                            ggcTGCCCTGCCGGGGGTTAAACCCACTGCGGTGAGAAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTAccc                                           >  1:55016/1‑149 (MQ=255)                                                                                                                   ctgccGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACCGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTg                                   <  2:295357/149‑1 (MQ=255)                                                                                                                   ctgccGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGGCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTg                                   <  1:134894/149‑1 (MQ=255)                                                                                                                       cGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGcc                                >  1:54247/1‑148 (MQ=255)                                                                                                                                aaCCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGAc                      <  8:83111/149‑1 (MQ=255)                                                                                                                                  cccACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGc                                                     <  6:134165/116‑1 (MQ=255)                                                                                                                                  cccACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGc                                                     <  6:128250/116‑1 (MQ=255)                                                                                                                                   ccACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTTAAACCCCGGCCACAGGGTGATCACCGGCGCCAGTATGCGCCGGGGTACCCGCGCAGGGGGCCAAGCACGGCAGc                   >  5:40580/1‑149 (MQ=255)                                                                                                                                   ccACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATAACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAg                    >  5:321048/1‑148 (MQ=255)                                                                                                                                     aCTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACACGCGAAGTGGCCAATGCAGGACAGca                  >  7:258013/1‑148 (MQ=255)                                                                                                                                                aaaaGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAGCACAATTCCCCtt       <  3:231115/149‑1 (MQ=255)                                                                                                                                                   aGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATcaca                                                                    >  2:323443/1‑85 (MQ=255)                                                                                                                                                   aGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATcaca                                                                    <  1:323443/85‑1 (MQ=255)                                                                                                                                                    gCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAGCACAATTCCCCTTCCa    >  2:315859/1‑148 (MQ=255)                                                                                                                                                     ccccGGGCGGTGAAAAAGCCTACTGCCGGTTACGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAGCACAATTCCCCTTCCATc  <  5:21012/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CACTGGGAAATGGCCCGCAAGCTGGCCGAGGCCGAAGCGCTGGCGCCCAAGGCAGCCCCGCGCAAGAAGGCAGCGGCCAAGCCGAAGGCGCCTGCCGAACCCAAGGCGGCTGCCCTGCCGGGGGTTAAACCCACTGCGGTGA‑AAAAGCCCCGGGCGGTGAAAAAGCCTACTGCCGGTTAAGCCTGTCC‑‑GGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTTCATCACAGGCGCATGTAGGAGCGGGTTTACCCGCGAAGTGGCCAATGCAGGACAGCACAATTCCCCTTCCATC  >  NC_002947/4573013‑4573308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑