Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,304,452	G→A	6.1%	R363H (CGC→CAC)	cysS →	cysteine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,304,452	0	G	A	6.1%	195.7 / 2.0	65	R363H (CGC→CAC)	cysS	cysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (2/2);  total (30/35)
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.17e-01

GCTGGAGCGCTTCTACCACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTGGTGCCGAAGGCAAGGTC  >  NC_002947/3304309‑3304591                                                                                                                                                |                                                                                                                                            gcTGGAGCGCTTCTACCACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCgaga                                                                                                                                        <  7:280834/149‑1 (MQ=255)       gcgcTTCTACCACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAAcc                                                                                                                                  >  3:143997/1‑149 (MQ=255)       gcgcTTCTACAACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAAcc                                                                                                                                  >  6:50787/1‑149 (MQ=255)           ttCTACCACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCt                                                                                                                              >  5:230546/1‑149 (MQ=255)                 cACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGa                                                                                                                        >  7:135893/1‑149 (MQ=255)                   cGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACa                                                                                                                      <  5:232725/149‑1 (MQ=255)                   cGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACa                                                                                                                      <  6:230546/149‑1 (MQ=255)                    gCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAg                                                                                                                     >  7:58701/1‑149 (MQ=255)                           gcgGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTg                                                                                                              <  7:266098/149‑1 (MQ=255)                            cgGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTgg                                                                                                             <  8:135893/149‑1 (MQ=255)                              gCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAg                                                                                                           <  6:325973/149‑1 (MQ=255)                               ccTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAgg                                                                                                          <  1:327775/149‑1 (MQ=255)                                 tGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGcc                                                                                                        >  4:144571/1‑149 (MQ=255)                                  gCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCg                                                                                                       >  5:1930/1‑149 (MQ=255)                                      cGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCgctgc                                                                                                   <  2:80427/149‑1 (MQ=255)                                        ggTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGctgctg                                                                                                 <  7:62896/149‑1 (MQ=255)                                         gTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTCGCGCACCAGATcgaac                                                                                                                                  >  3:168459/1‑111 (MQ=255)                                         gTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTCGCGCACCAGATcgaac                                                                                                                                  <  4:168459/115‑5 (MQ=255)                                          tGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGctggc                                                                                               <  8:58701/149‑1 (MQ=255)                                             cggcCAATGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGctggcctg                                                                                            >  7:237803/1‑149 (MQ=255)                                               gcCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGctggcctggc                                                                                          >  8:19097/1‑149 (MQ=255)                                               gcCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGctggcctggc                                                                                          <  5:35249/149‑1 (MQ=255)                                                ccAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTCGCGCACCAGATcgaac                                                                                                                                  >  2:360108/1‑104 (MQ=255)                                                ccAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTCGCGCACCAGATcgaac                                                                                                                                  <  1:360108/108‑5 (MQ=255)                                                   aGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGt                                                                                      <  5:50787/149‑1 (MQ=255)                                                    gggTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTc                                                                                     <  6:103518/149‑1 (MQ=255)                                                       tGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGAc                                                                                  >  1:292431/1‑149 (MQ=255)                                                         gCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTg                                                                                >  6:362287/1‑149 (MQ=255)                                                          cGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGTGGCCGCTGCTGGCCTGGCCGGTCGACTgc                                                                               >  5:155583/1‑149 (MQ=255)                                                           gAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTgcg                                                                              >  4:65732/1‑149 (MQ=255)                                                            aaGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGTGGCCGCTGCTGGCCTGGCCGGTCGACTgcgc                                                                             <  6:155583/149‑1 (MQ=255)                                                               cGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAg                                                                          >  7:46537/1‑149 (MQ=255)                                                               cGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAg                                                                          <  7:203152/149‑1 (MQ=255)                                                                       gAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGa                                                                  >  4:318475/1‑149 (MQ=255)                                                                        aGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAg                                                                          <  2:312798/140‑1 (MQ=255)                                                                        aGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAg                                                                          >  1:312798/1‑140 (MQ=255)                                                                        aGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGAt                                                                 >  4:56158/1‑149 (MQ=255)                                                                          cgcTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGt                                                               <  4:143997/149‑1 (MQ=255)                                                                             ttCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCg                                                                                                                   >  5:39527/1‑94 (MQ=255)                                                                             ttCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCg                                                                                                                   <  6:39527/94‑1 (MQ=255)                                                                               cAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTg                                                           <  2:13246/148‑1 (MQ=255)                                                                                aGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTgg                                                          <  3:318475/148‑1 (MQ=255)                                                                                aGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTgg                                                          <  7:19097/148‑1 (MQ=255)                                                                                     ggCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCgc                                                                                                      >  1:329597/1‑99 (MQ=255)                                                                                     ggCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCgc                                                                                                      <  2:329597/99‑1 (MQ=255)                                                                                      gCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGAt                                                                 <  2:33432/135‑1 (MQ=255)                                                                                      gCGATGAACGACGACTTCGGCACCACCGAACCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGAt                                                                 >  1:33432/1‑135 (MQ=255)                                                                                        gATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGAt                                                                 <  1:209426/133‑1 (MQ=255)                                                                                        gATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGAt                                                                 >  2:209426/1‑133 (MQ=255)                                                                                        gATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGtgctg                                                  >  3:233272/1‑148 (MQ=255)                                                                                             acgacgACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTg                                            <  6:1930/149‑1 (MQ=255)                                                                                               gacgacTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTgc                                                                               <  2:193788/112‑1 (MQ=255)                                                                                               gacgacTTCGGCACCCCCGAAGCCTGGGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTgc                                                                               >  1:193788/1‑112 (MQ=255)                                                                                                       cGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGAt                                   <  3:144571/148‑1 (MQ=255)                                                                                                            gccccGACGCCGGCGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTc                             <  8:237803/148‑1 (MQ=255)                                                                                                             ccccGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTgg                                           >  2:56062/1‑134 (MQ=255)                                                                                                             ccccGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTgg                                           <  1:56062/134‑1 (MQ=255)                                                                                                                      cTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTg                   <  7:224972/149‑1 (MQ=255)                                                                                                                      cTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTg                   <  3:65732/149‑1 (MQ=255)                                                                                                                      cTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATAACTTCCTGCGTGCTg                   <  2:292431/149‑1 (MQ=255)                                                                                                                            ggTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTGGTGcc              >  3:169488/1‑148 (MQ=255)                                                                                                                               ccTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTgg                                                          >  2:263897/1‑101 (MQ=255)                                                                                                                               ccTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTgg                                                          <  1:263897/101‑1 (MQ=255)                                                                                                                                  gTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTGGTGCCGAAGGc        >  5:122593/1‑148 (MQ=255)                                                                                                                                       aTCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTGGTGCCGAAGGCAAGGTc  >  5:239699/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                            GCTGGAGCGCTTCTACCACGCCCTGCGCGGCCTGCCACGGGTGGCGGCCAAGGGTGGCGAAGCGTTTGTCGAGCGCTTCAGCGTGGCGATGAACGACGACTTCGGCACCCCCGAAGCCTGTGCGGTCCTGTTCGATCTGGTGCGCGAGATCAACCGCCTGCGCGACAGCGACGTGGAGGCCGCTGCTGGCCTGGCCGGTCGACTGCGCGAGCTGGGTGATGTACTGGGTGTGCTGCAACTGGAGGCCGATGACTTCCTGCGTGCTGGTGCCGAAGGCAAGGTC  >  NC_002947/3304309‑3304591

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