Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	598,428	T→C	19.1%	L317P (CTC→CCC)	cusB →	copper/silver efflux system, membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	598,428	0	T	C	19.1%	100.1 / 2.1	56	L317P (CTC→CCC)	cusB	copper/silver efflux system, membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (34/10);  new base C (4/8);  total (38/18)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGAAGCCAATATTTCTGGCGCAC  >  NC_000913/598201‑598661                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           ggACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCACCTGCGTCTGGAAGTCGACAACCCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTCCAACTCAACA‑CCGcc                                                                                                                                                                                                                                 >  1:481115/1‑238 (MQ=255)   gcccGGTTTGGGTCACTGCTGGGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAg                                                                                                                                                                                                                               <  1:1161133/237‑1 (MQ=255)        gtgtGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGCCGAGGCGCTCAAACCGGGCATGAACCCCTGGTTGCAACTCAACA‑CCGCCCGCGAAcc                                                                                                                                                                                                                         >  1:1046302/1‑239 (MQ=255)        gtgtGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAAcc                                                                                                                                                                                                                         >  1:217309/1‑239 (MQ=255)        gtgtGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCCCCCCCTGCAGCTGCGCCTGGCAGTCGACAACGCCAGCCAGCCCCTCAAACCGGGCCTGAACCCCTGGTTCCACCTCAAAC‑ACGCCCGCAAAcc                                                                                                                                                                                                                         >  2:340982/1‑239 (MQ=255)                   tgctgcGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGAtgcggccaa                                                                                                                                                                                                              >  1:487623/1‑236 (MQ=255)                   tgctgcGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCa                                                                                                                                                                                                              >  1:500889/1‑239 (MQ=255)                   tgctgcGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAAAGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACCAAACA‑CCGCCAGCGAACCGATGCTGCCCa                                                                                                                                                                                                              >  1:641498/1‑239 (MQ=255)                        cGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGTCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCCCGGAGGCGCTAAAACCGGGCATTAACGCCGGGTTTCAAACCACAA‑CGCCCAGCGAACCGAtgcggcccctcccg                                                                                                                                                                                                         >  2:1116070/1‑231 (MQ=255)                            ccgggAGTTTAACGCCTGGCTGGGGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACCCTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCGGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACCCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTcac                                                                                                                                                                                                     <  1:932152/236‑1 (MQ=255)                            cccgggatcattCGCCGGGCTGGTGAAAGATGCCGCGCCGTTTACCCCCCCCGTTCCGGCGGGACCGGATAAAACACTCCCCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTcac                                                                                                                                                                                                     <  1:131787/228‑1 (MQ=255)                             ccGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACACCGCCGACGAGGCGCTAAAACCGGGCCTGCCGGCCTGGTTGCAAACCAACA‑CCGCCAACGAACCGACGCTGCCAACTCCGTcaca                                                                                                                                                                                                    >  2:1142714/1‑239 (MQ=255)                                     tATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCtga                                                                                                                                                                                            >  1:1213990/1‑239 (MQ=255)                                             ggctggTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGCCGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCCACAAACCGATGCGCCCCATTCCGGCACAAACCCTGATAGATAcg                                                                                                                                                                                    >  2:75352/1‑238 (MQ=255)                                              gctggTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGCAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCCGCGAACCGATGCTGCCCATTCCGTCACAAGCGCTCATTGATACCg                                                                                                                                                                                   >  1:942904/1‑239 (MQ=255)                                                       aGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACCCTCACCATCCGCAAATGGACGCTTCTCCCTGGCGGGGATGCCGCGACCCGCACCCGGCAGCTCGGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCCGGCGGGAACGCCGTGTTGCACACCAACC‑ACGCCAGCGAACCCATGCGCCCCATTCCGTCCCAACCGCTGACTGATAccccccgccaca                                                                                                                                                                          >  2:714645/1‑229 (MQ=255)                                                        gATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCCTTCCGTCCCAACCCCTGATGAATACCGGCAGGCAACa                                                                                                                                                                         >  1:1005113/1‑239 (MQ=255)                                                            ccTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCACCGAACCGATGCTGCACATTCCGTAACAAGCGCTGATTGAAACGGCCAGGGAACAgggg                                                                                                                                                                     >  2:598752/1‑239 (MQ=255)                                                            ccTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCCCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGCCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCCCATTCCGTCACAAGCGCCGATTGATACCGGCAGCGAAACGCgg                                                                                                                                                                     >  2:378314/1‑239 (MQ=255)                                                             cTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGACACGATGCTGCTCATTCCGTCACAAGGGCCGATTGATACCGGCACCGAACAACggg                                                                                                                                                                    >  1:1230658/1‑239 (MQ=255)                                                                            cACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCCCAAGCGCTGATTGATACCGGCCGCGACACGCGGGGGATAACCGtttgt                                                                                                                                                      >  1:74142/1‑235 (MQ=255)                                                                              ccGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGTTACCGGCAGCGAACAGCGGGGGCTTACCGTTGTTGc                                                                                                                                                    >  1:1208112/1‑238 (MQ=255)                                                                                  tCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGCCAGCGAACAACGGGGGATTACCGCTGATCCCGACg                                                                                                                                               >  1:749429/1‑239 (MQ=255)                                                                                           aCCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCTTTGATGCCGCCGGGCGCTTTg                                                                                                                                      >  1:512759/1‑239 (MQ=255)                                                                                                         tCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCGGCCCCTTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGGTTACCGTTTGTGCCGGCGGGGGCTTTTAACCGAAAcggggg                                                                                                                        >  2:824777/1‑235 (MQ=255)                                                                                                          cACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAAAAGCGGGTGATTACCGTTGATGCCGACGGGAGCTTTGTACCGAAACGCGTTg                                                                                                                       >  1:530110/1‑239 (MQ=255)                                                                                                          cACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAATCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCAAACCGCTGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCCAACAGCGGGGCATTCCCGTTGTTGCCAACGGGCCCCTTGTACCACAAcgcagtg                                                                                                                       >  2:486374/1‑235 (MQ=255)                                                                                                            ccATCCGAAAAGGGACGCGGCTACCTGGCGTGGATGCCGCGACCCGGACGCTGCAGCTGCGTCTGGAAGTCGAAAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCt                                                                                                                     <  1:340982/239‑1 (MQ=255)                                                                                                            ccAACCGCAAATGGACGCTGCTACCTGGGGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCt                                                                                                                     <  1:378314/239‑1 (MQ=255)                                                                                                              aTCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAACGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGt                                                                                                                   >  1:369143/1‑239 (MQ=255)                                                                                                                ccGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGttt                                                                                                                 >  1:261521/1‑239 (MQ=255)                                                                                                                ccGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCACCGAACAGCGGGGACTTACCGTTGTTGCCGACGGGGGCTTTGTACCGAAAAGCGGGTCTGttt                                                                                                                 >  2:129768/1‑239 (MQ=255)                                                                                                                       tGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTAGAAGTCGACAACGCCGCCGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCGCATTCCTCTACAAGCGCTGATTGCTACCGGCAGCGAACAGCGGGGGATTACCGTTGATTCCCACGGGCGCTTTGAACACAAAAGGGTTCAGGGTTTCCAgga                                                                                                          >  2:663166/1‑238 (MQ=255)                                                                                                                                        cGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCAcc                                                                                          >  1:986220/1‑238 (MQ=255)                                                                                                                                        cGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGGGATTACCGTTGATGCCGCCGGGCGCCTTGTACCGAAACCGGTGGCTTGTTTCCAGGGAACCCAAGGGGtaccc                                                                                          >  2:730118/1‑234 (MQ=255)                                                                                                                                          tGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑ACGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAAAAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCCCAAGGGGTCTCCGCa                                                                                       >  1:889695/1‑239 (MQ=255)                                                                                                                                            gATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACCAGCGCTGATTGATACCGGCAGCGAACAGCGGGTTGTTACCGTTGATGCCGACGCGCGCTTTGTACCGAAACGCGTTGCCGTTTTCCAGGCATCGCACGGCGTCACCGCCtt                                                                                     >  1:429034/1‑239 (MQ=255)                                                                                                                                              tcccGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCCGCGAACCGATGATGCTCATTCCGTCACAAGCGCTGATTGATACCGCCAGCGAACAGCGGGTGATTACCGTTATCGCCACCGTGCGCTCGTTACCCAAACGCGTTGCTGCTTTACAGGCCTTGCAAAGAGTTAccccgtcac                                                                                   >  2:263055/3‑232 (MQ=255)                                                                                                                                                              cgcaaGCTGCGTCTGGAAGCCGAAAACGCCGACGGGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACCCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGa                                                                   <  1:730118/235‑1 (MQ=255)                                                                                                                                                                       cGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTATGGCGGAAGGTGgaaa                                                          >  1:620025/1‑239 (MQ=255)                                                                                                                                                                         tCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAgg                                                        >  1:278325/1‑239 (MQ=255)                                                                                                                                                                         tCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGGGAAAAgg                                                        >  1:88393/1‑239 (MQ=255)                                                                                                                                                                         tCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGCTTGGCGGAAGGTGCAAAgg                                                        >  1:243034/1‑239 (MQ=255)                                                                                                                                                                         ccTGGAAGCCGACAACGCCGAGGGGGCCCTAAAACCGGGCATGAACCCCTGGTTCCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATA‑CGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAgg                                                        <  1:75352/237‑1 (MQ=255)                                                                                                                                                                                   gacgacGCCGACGAGGCCATAAAACGGGGCAGAACCGCCTGGTTGCAACCCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCa                                               <  1:163914/238‑1 (MQ=255)                                                                                                                                                                                   gACACCGCCGACGGGGGCGTAAACCCGGGAAGGAACGCCTGGGTGCAACCCAACA‑CCGCCAGCGAACCGATGCGGCTCATTCCGGCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCa                                               <  1:186307/238‑1 (MQ=255)                                                                                                                                                                                       acccggccggggCGCAAAAACCGGCACTGACCGCCTGGATCAAACCCACAC‑CCGCCAGCGAACCGATGCTGCTCATTCCGGCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGTGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGc                                          <  2:852543/229‑1 (MQ=255)                                                                                                                                                                                       aCGCCGAGGGGGCGCTAAAACCGGGCATGAACGCCTGGTGGCAACCCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGc                                          <  1:297038/239‑1 (MQ=255)                                                                                                                                                                                       aCGCCGACGGGGCCTTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGc                                          <  1:1303292/239‑1 (MQ=255)                                                                                                                                                                                       aCGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGc                                          <  1:78534/239‑1 (MQ=255)                                                                                                                                                                                                        aacccGGGAAGGAACGCCTGGTTGCAACCCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTc                         <  1:1152806/236‑1 (MQ=255)                                                                                                                                                                                                         acccGGGAAGGAACGCCTGGTTGCAACTAAACC‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTcc                        <  1:486374/237‑1 (MQ=255)                                                                                                                                                                                                            ggggCATGACCGCGGGGTGGCAACTCAACA‑CCGCCAGCGAACCGAGGCTGCTCATTCCGTCACAAGCGCTGATTGATCCCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGaa                     <  1:663166/238‑1 (MQ=255)                                                                                                                                                                                                                 cttaccctcttgTTGCAACCCCACCCCCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGTTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGAAGCCa                 <  2:641498/228‑1 (MQ=255)                                                                                                                                                                                                                    aaCGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGGGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGAAGCCCatat             >  1:1063780/1‑239 (MQ=255)                                                                                                                                                                                                                               gCAACTCACCA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGAAGCCAATATTTCTGGCGCAc  <  1:1280993/239‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           GGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACA‑CCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCCGAAGCCAATATTTCTGGCGCAC  >  NC_000913/598201‑598661

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑