Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	598,658	A→C	5.5%	F96C (TTT→TGT)	trx ←	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	598,658	0	A	C	5.5%	235.6 / 4.1	73	F96C (TTT→TGT)	trx	thioredoxin
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base C (2/2);  total (37/36)
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.79e-01

CCTGTAGCAGCGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCTTTGC  >  NC_002947/598512‑598800                                                                                                                                                   |                                                                                                                                               ccTGTAGCAGCGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCaaaa                                                                                                                                              <  3:347500/149‑1 (MQ=255) ccTGTAGCAGCGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCaaaa                                                                                                                                              >  1:173832/1‑149 (MQ=255)     tAGCAGCGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCg                                                                                                                                           <  1:194446/148‑1 (MQ=255)       gcagcGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGt                                                                                                                                          >  2:306191/1‑147 (MQ=255)       gcagcGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGt                                                                                                                                          <  1:306191/147‑1 (MQ=255)             tttcggtttcggCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTTTGCCGGGGCACAACCg                                                                                                                                           <  8:234395/140‑1 (MQ=255)             tttcggtttcggCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTTTGCCGGGGCACAACCg                                                                                                                                           >  7:234395/1‑140 (MQ=255)             tttcggtttcggCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGc                                                                                                                                  >  8:206278/1‑149 (MQ=255)               tcggtttcggCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTg                                                                                                                                >  5:136285/1‑149 (MQ=255)                       ggCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCt                                                                                                                        >  6:270694/1‑149 (MQ=255)                         cAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGAGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATAGCCGATTCCGGTAGTGCACCGGCAACACCGTCCCCCGGCTGGCCGTCCtt                                                                                                                       >  5:100386/1‑148 (MQ=255)                              cGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAg                                                                                                                 <  6:136285/149‑1 (MQ=255)                                gCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGca                                                                                                               <  4:374576/149‑1 (MQ=255)                                 cTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGcac                                                                                                              >  4:198215/1‑149 (MQ=255)                                  tttCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAAc                                                                                                                   <  1:178210/143‑1 (MQ=255)                                  tttCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAAc                                                                                                                   >  2:178210/1‑143 (MQ=255)                                    tCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGcaccac                                                                                                           >  8:33098/1‑149 (MQ=255)                                    tCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGcaccac                                                                                                           >  1:155539/1‑149 (MQ=255)                                       gCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGt                                                                                                        >  7:73465/1‑149 (MQ=255)                                       gCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAAAAGCACCACTGt                                                                                                        >  3:343007/1‑149 (MQ=255)                                        cAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTc                                                                                                       <  5:270694/149‑1 (MQ=255)                                          aaCAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCgg                                                                                                     >  7:37229/1‑149 (MQ=255)                                           aCAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGc                                                                                                    <  2:347550/149‑1 (MQ=255)                                              gTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAgg                                                                                                 <  8:73465/149‑1 (MQ=255)                                                   tttGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACg                                                                                            >  2:344099/1‑149 (MQ=255)                                                    ttGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACg                                                                                            >  3:334175/1‑148 (MQ=255)                                                     tGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACgg                                                                                           <  6:208355/148‑1 (MQ=255)                                                       ccTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTc                                                                                                       <  4:36995/134‑1 (MQ=255)                                                       ccTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTc                                                                                                       >  3:36995/1‑134 (MQ=255)                                                       ccTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATg                                                                                        >  3:376462/1‑149 (MQ=255)                                                        cTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTgg                                                                                                                               <  3:175911/109‑1 (MQ=255)                                                        cTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTgg                                                                                                                               >  4:175911/1‑109 (MQ=255)                                                         tGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGcc                                                                                      >  3:137702/1‑149 (MQ=255)                                                           cTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCg                                                                                     <  2:173832/148‑1 (MQ=255)                                                             ccAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAAc                                                                                  >  3:283596/1‑149 (MQ=255)                                                                gTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTgg                                                                               <  1:255331/149‑1 (MQ=255)                                                                gTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTgg                                                                               <  5:175438/149‑1 (MQ=255)                                                                gTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTgg                                                                               <  4:283596/149‑1 (MQ=255)                                                                      cggcggTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAac                                                                          >  6:344756/1‑148 (MQ=255)                                                                        gcggTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAaccac                                                                       >  4:278772/1‑149 (MQ=255)                                                                          ggTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAATGCGTCGACGGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCt                                                                     <  4:23543/149‑1 (MQ=255)                                                                             gggCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGtt                                                                  >  5:188113/1‑149 (MQ=255)                                                                             gggCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGtt                                                                  <  2:221510/149‑1 (MQ=255)                                                                               gcgcGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGc                                                                <  4:137702/149‑1 (MQ=255)                                                                               gcgcGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGc                                                                <  2:155539/149‑1 (MQ=255)                                                                                 gcggcgGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTc                                                              <  4:376462/149‑1 (MQ=255)                                                                                              gTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGtt                                                  <  8:37229/148‑1 (MQ=255)                                                                                                 cATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCa                                                            <  5:363792/135‑1 (MQ=255)                                                                                                 cATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCa                                                            >  6:363792/1‑135 (MQ=255)                                                                                                     cGGCTCGATCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAACACGTCCACCGGCTTGCCGTCCTTGAACAGAACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTg                                           <  6:100386/148‑1 (MQ=255)                                                                                                     cGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTg                                           <  5:344756/148‑1 (MQ=255)                                                                                                      ggCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGt                                                                                                        <  4:316964/86‑1 (MQ=255)                                                                                                      ggCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGt                                                                                                        >  3:316964/1‑86 (MQ=255)                                                                                                              cATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCcagcagca                                 <  6:188113/149‑1 (MQ=255)                                                                                                               aTGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCcagcagcag                                >  4:233668/1‑149 (MQ=255)                                                                                                                tGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGt                               <  7:109673/149‑1 (MQ=255)                                                                                                                  gcgcGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTc                             <  3:233668/149‑1 (MQ=255)                                                                                                                    gcgGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGc                           <  3:198215/149‑1 (MQ=255)                                                                                                                     cgGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGcc                          <  4:334175/149‑1 (MQ=255)                                                                                                                        aTCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTg                       >  3:64718/1‑149 (MQ=255)                                                                                                                          cGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTgg                      >  8:360847/1‑148 (MQ=255)                                                                                                                            ccGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTaa                   >  7:285926/1‑149 (MQ=255)                                                                                                                                   cGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAg                                                    <  2:12605/109‑1 (MQ=255)                                                                                                                                   cGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAg                                                    >  1:12605/1‑109 (MQ=255)                                                                                                                                    ggTTGTCCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGACAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCgg            >  1:172111/1‑148 (MQ=255)                                                                                                                                     gTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGc           >  5:372427/1‑148 (MQ=255)                                                                                                                                        gtgCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAAt        <  3:278772/148‑1 (MQ=255)                                                                                                                                          gCCCCGGCAAAACCGTCCACCGCCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCCCCAGTTGCTCCACGACGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCtt     >  3:91187/1‑149 (MQ=255)                                                                                                                                           ccggggCACAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTg                                                 <  1:93373/100‑1 (MQ=255)                                                                                                                                           ccggggCACAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTg                                                 >  2:93373/5‑104 (MQ=255)                                                                                                                                           ccccGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCg                            >  2:137765/1‑125 (MQ=255)                                                                                                                                           ccccGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCg                            <  1:137765/125‑1 (MQ=255)                                                                                                                                           ccccGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCttt    >  8:212202/1‑149 (MQ=255)                                                                                                                                            cccGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCTTTg   >  1:171480/1‑149 (MQ=255)                                                                                                                                             ccGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCTTTGc  <  6:262273/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCTGTAGCAGCGTTTCGGTTTCGGCAAAACGGCTTTCGGCAAACAGTGCCTTTGCCTGCTCCAGTGGCGCGGCGGTGGGCGCGGCGGGCATTTGTACATGCGGCTCGAGCATGGCGCGGATCGCCGATTCCGGTTGTGCCCCGGCAAAACCGTCCACCGGCTGGCCGTCCTTGAACAGCACCACTGTCGGCAGGCTACGGATGCCGAACTGGGCAACCACCTGTTGCTCCACGTCGCAGTTGATCTTGGCCAGCAGCAGTTCGCCCTGGTAACCCTCGGCAATCTTTGC  >  NC_002947/598512‑598800

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