Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	618,346	C→T	9.4%	V75M (GTG→ATG)	PP_0532 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	618,346	0	C	T	9.4%	176.1 / 7.7	64	V75M (GTG→ATG)	PP_0532	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base T (2/4);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGT  >  NC_002947/618201‑618492                                                                                                                                                   |                                                                                                                                                   tCCAGCGTGCGGCCGGGCCTGTCGATGCGCAGTTCGTTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCg                                                                                                                                                  <  5:21871/148‑1 (MQ=255) tCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCgg                                                                                                                                                 >  8:259842/1‑149 (MQ=255)   cAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCggg                                                                                                                                                >  8:61142/1‑148 (MQ=255)   cAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCggg                                                                                                                                                <  3:78358/148‑1 (MQ=255)   cAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGt                                                                                                                                               <  1:220787/149‑1 (MQ=255)       gTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGTTCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTAccc                                                                                                                                           >  1:252317/1‑149 (MQ=255)            gccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATg                                                                                                                                      >  5:54739/1‑149 (MQ=255)                 cccTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgtg                                                                                                                                 <  8:236946/149‑3 (MQ=255)                   cTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgtgga                                                                                                                               <  8:260380/149‑5 (MQ=255)                      tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgtggatg                                                                                                                             <  7:259842/148‑7 (MQ=255)                      tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                             <  2:289926/148‑1 (MQ=255)                      tCGATGCGCAGTTCGGTCCGGTAACCCACCCCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATGCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                             <  5:247097/148‑1 (MQ=255)                       ccaggacctggcctatcacttcacccagcttgatctggactcggccgatggcctgcgtcactaccgcctgtgggttggcaagccgaatcggCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTgg                                                                                                                            >  8:270576‑M1/92‑149 (MQ=255)                         aTGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTgg                                                                                                                            >  2:283176/1‑146 (MQ=255)                         aTGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTgg                                                                                                                            <  1:283176/146‑1 (MQ=255)                           gcgcAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGAt                                                                                                                       >  5:237164/1‑149 (MQ=255)                             gcAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTc                                                                                                                     <  2:229081/149‑1 (MQ=255)                                 ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgt                                                                                                                                  <  8:292973/132‑2 (MQ=255)                                 ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgt                                                                                                                                  >  7:292973/1‑131 (MQ=255)                                   cGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCtt                                                                                                                <  2:252317/148‑1 (MQ=255)                                     ggccGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGcc                                                                                                             <  4:197309/147‑1 (MQ=255)                                      tCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTg                                                                                                                                   <  2:248944/126‑1 (MQ=255)                                      tCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTg                                                                                                                                   <  2:252513/126‑1 (MQ=255)                                         ggTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAcc                                                                                                         >  8:19994/1‑149 (MQ=255)                                         ggTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAcc                                                                                                         >  2:12740/1‑149 (MQ=255)                                               ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTg                                                                                                               <  5:110440/137‑1 (MQ=255)                                               ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTg                                                                                                               >  6:110440/1‑137 (MQ=255)                                                cACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCGCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGAt                                                                                                                                       >  1:59016/1‑113 (MQ=255)                                                cACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGAt                                                                                                                                       <  2:59016/112‑1 (MQ=255)                                                        ccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGa                                                                                          <  1:12740/149‑1 (MQ=255)                                                               gggCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGc                                                                                   <  7:61142/149‑1 (MQ=255)                                                                  gccgatggcctgcgtcactaccgcctgtgggttggcaagccgaatcgtCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                             >  6:37287‑M1/49‑105 (MQ=255)                                                                  gccgatggcctgcgtcactaccgcctgtgggttggcaagccgaatcgtCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                             <  5:37287‑M1/57‑1 (MQ=255)                                                                    ccTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCg                                                                              <  6:237164/149‑1 (MQ=255)                                                                          ccACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGGATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAg                                                                        <  4:288370/149‑1 (MQ=255)                                                                          ccACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAg                                                                        >  4:57039/1‑149 (MQ=255)                                                                          ccACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAg                                                                        <  6:301155/149‑1 (MQ=255)                                                                           cACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                 >  5:50952/1‑123 (MQ=255)                                                                           cACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                 <  6:50952/123‑1 (MQ=255)                                                                               gCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTc                                                                      >  3:272054/1‑146 (MQ=255)                                                                                cGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGa                                                                  <  6:54739/149‑1 (MQ=255)                                                                                       gTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                           >  8:268975/1‑85 (MQ=255)                                                                                       gTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                           <  7:268975/85‑1 (MQ=255)                                                                                       gTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCt                                                           >  6:150383/1‑149 (MQ=255)                                                                                            aGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGt                                                       <  7:19994/149‑1 (MQ=255)                                                                                                 ctgctgATCCAGCGTCACCAGCGCAGCATTCCCATCTAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                   <  7:159951/148‑1 (MQ=255)                                                                                                    ctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                                                >  8:50053/1‑148 (MQ=255)                                                                                                    ctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGAt                                               <  3:57039/149‑1 (MQ=255)                                                                                                       aTCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAgg                                            >  5:4991/1‑149 (MQ=255)                                                                                                         ccAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGc                                           >  7:35175/1‑148 (MQ=255)                                                                                                           aGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAg                                        >  2:20787/1‑149 (MQ=255)                                                                                                           aGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAg                                        >  4:230417/1‑149 (MQ=255)                                                                                                            gCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAgg                                       >  5:225033/1‑149 (MQ=255)                                                                                                                  CCAGCACAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                 <  7:260570‑M1/149‑1 (MQ=255)                                                                                                                     GCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGgcgc                              >  1:20787‑M1/1‑149 (MQ=255)                                                                                                                        CAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGt                           >  4:272054‑M1/1‑149 (MQ=255)                                                                                                                        CAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGt                           >  6:4991‑M1/1‑149 (MQ=255)                                                                                                                               cccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCg                    >  1:131159/1‑149 (MQ=255)                                                                                                                                ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGa                   >  6:242801/1‑149 (MQ=255)                                                                                                                                     gAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACgtgt              <  5:242801/149‑1 (MQ=255)                                                                                                                                           ccACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCa                                         >  2:258235/1‑116 (MQ=255)                                                                                                                                           ccACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCa                                         <  1:258235/116‑1 (MQ=255)                                                                                                                                              caGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttct     <  7:50053/149‑1 (MQ=255)                                                                                                                                                 cACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGt  >  8:232109/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTC‑CAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGT  >  NC_002947/618201‑618492

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