Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,319	A→G	8.6%	A224A (GCA→GCG)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,319	0	A	G	8.6%	156.7 / 4.9	59	A224A (GCA→GCG)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (3/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645173‑5645467                                                                                                                                                   |                                                                                                                                                      cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                                       <  5:9241/147‑1 (MQ=255) cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                                       >  6:9241/1‑147 (MQ=255) cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                                     >  5:17126/1‑149 (MQ=255)    agCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCTGCGTTCAGCGCGGCC‑ac                                                                                                                                                 <  5:106341/149‑2 (MQ=255)      ccTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGCGGTGGCCGCGCTGAACGCAGCC‑AAgg                                                                                                                                               <  1:8638/149‑1 (MQ=255)          gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                            <  6:172980/148‑1 (MQ=255)          gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                            <  4:283318/148‑1 (MQ=255)              gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACTGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGatcatc                                                                                                                                       <  8:287798/149‑1 (MQ=255)               aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                            <  4:263112/143‑1 (MQ=255)               aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                            >  3:263112/1‑143 (MQ=255)                cTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATcc                                                                                                                                      <  7:276147/148‑1 (MQ=255)                      ttcGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCg                                                                                                                               >  5:90306/1‑149 (MQ=255)                          aGCTGCTCAGCTACGTGAAGCTCAAGTAAATCCTGTTTTTTACGGTCCGCTCGCTCGACGTCTAGGCAAGGCGCCTGGAACTGGCCCAGCAGAAGGGGAAGGTGGTGGAGGGGCTGATCGGTGCC‑AAGGAGATCTTCCTGGAGCGGGTg                                                                                                                           <  8:236219/149‑1 (MQ=255)                                         tGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacc                                                                                                             >  6:60923/1‑148 (MQ=255)                                          gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacca                                                                                                            >  6:195988/1‑148 (MQ=255)                                             gttCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                                        >  6:90464/3‑149 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                       >  8:252379/1‑148 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                       >  2:43673/1‑148 (MQ=255)                                                 aaGGCCAGCCTGTTTGACCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCg                                                                                                      >  1:129732/1‑147 (MQ=255)                                                 aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGa                                                                                                     <  5:195988/148‑1 (MQ=255)                                                 aaGGCCAGCCTGTCTGATCCCGTCCGGCCGCTCGACCACGACGCTCGGCGCCTGGAACTGGCCCAGCAGTACGGCAAGGTGGTGGACGCGCTGAACGCAGCC‑AAGCAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                                    <  6:164/149‑1 (MQ=255)                                                    gCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                                 <  5:138587/149‑1 (MQ=255)                                                        gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                             <  5:60923/149‑1 (MQ=255)                                                                gAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                                     >  3:165729/1‑149 (MQ=255)                                                                  aCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGCACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                                    >  4:172981/1‑148 (MQ=255)                                                                   ccGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                                    >  3:13846/1‑147 (MQ=255)                                                                    cGGTCCGGACGCTCGACGTCGAGGCACGGCGCCTGGAACTGTCCCAGGAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGTTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                                 <  6:90306/149‑1 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCGAGCCCACCTCGAAGGTTAGCCGCTGCCTCa                                                                                   >  8:211064/1‑145 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCa                                                                                   <  7:211064/145‑1 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCt                                                                               >  8:206521/1‑149 (MQ=255)                                                                         cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                             <  1:43673/148‑1 (MQ=255)                                                                          gCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTAc                                                                           <  3:294666/149‑1 (MQ=255)                                                                              gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                               >  6:201418/1‑109 (MQ=255)                                                                              gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                               <  5:201418/109‑1 (MQ=255)                                                                                tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTg                                                                     >  1:312403/1‑149 (MQ=255)                                                                                       cGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTTGCCCAgg                                                              <  7:206521/149‑1 (MQ=255)                                                                                          ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                             >  5:204582/1‑131 (MQ=255)                                                                                          ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                             <  6:204582/131‑1 (MQ=255)                                                                                             aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCTAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATc                                                         >  1:249052/1‑149 (MQ=255)                                                                                                  gcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAg                                                   >  3:88431/1‑149 (MQ=255)                                                                                                      ggAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgtg                                               >  5:141565/1‑149 (MQ=255)                                                                                                       gAACTGGCCCAGCAGAACGGCAAGGTTGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                                <  2:129732/147‑1 (MQ=255)                                                                                                               ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                                      <  4:165729/149‑1 (MQ=255)                                                                                                                    agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                                 <  5:90464/149‑1 (MQ=255)                                                                                                                    agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                                 <  4:13846/149‑1 (MQ=255)                                                                                                                              ggtggtGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTAc                        >  2:54891/1‑148 (MQ=255)                                                                                                                               tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                              <  7:27213/137‑1 (MQ=255)                                                                                                                               tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                              >  8:27213/5‑141 (MQ=255)                                                                                                                                tggtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc                     >  7:9925/1‑149 (MQ=255)                                                                                                                                  ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg                    >  4:136751/2‑148 (MQ=255)                                                                                                                                  ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg                    >  2:97939/2‑148 (MQ=255)                                                                                                                                  gtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg                    >  2:310520/1‑148 (MQ=255)                                                                                                                                    ggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCAc                 >  5:71656/1‑149 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                            >  4:11081/1‑119 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                            <  3:11081/119‑1 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGa              <  4:88431/149‑1 (MQ=255)                                                                                                                                              aaCGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGc       <  6:141565/149‑1 (MQ=255)                                                                                                                                              aaCGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGc       <  2:249052/149‑1 (MQ=255)                                                                                                                                                 gCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGa                                                    >  6:204013/1‑101 (MQ=255)                                                                                                                                                 gCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGa                                                    <  5:204013/101‑1 (MQ=255)                                                                                                                                                 gCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtt    >  8:81838/1‑149 (MQ=255)                                                                                                                                                   aGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGGCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtttt  <  2:312403/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645173‑5645467

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