Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,501,500	C→T	L81L (CTG→TTG)	yjgZ →	uncharacterized protein YjgZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,501,500	0	C	T	96.4%	154.0 / ‑5.7	55	L81L (CTG→TTG)	yjgZ	uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/32);  minor base G (0/2);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAC  >  NC_000913/4501260‑4501742                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   ttGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATTAGCCTTCTGGGAACCGTAGTCTGGAGCCACATGCTTGGTTgacggacg                                                                                                                                                                                                                                           <  1:333818/250‑1 (MQ=255) ttGTTACCCCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTgacggacg                                                                                                                                                                                                                                           <  1:268967/250‑1 (MQ=255)       cctcctGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGGCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGAc                                                                                                                                                                                                                                    <  2:265149/251‑1 (MQ=255)       cctcctGGACCGCTACATGTCCTACCAGGCGGAGAGAGGGTTAGCCGGGGCGGACTTTCACCGACGGCAAGATCGTGCGTTTAAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTACGGTCGTTATGGGCCGAAAGGCCTGGTTGTGCGCAGGCTCAATGGTGGCCGGCAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGACTGGCGCCACATGCGTGGTTGACGGCCGACCTGAc                                                                                                                                                                                                                                    >  2:51223/1‑251 (MQ=255)           cgggACCGCTTCTTGTGCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAAGAATATCTGCGAGCGTGCTAGCCGTCCGGGCGGTATGGGACCAAATGGCTGGTTGGTCGCTGGTTCACTGGGGGCCGGGAACCGGGCGTCATAGATAATGAGCCTTCTGGGATGCGCAGTCTGGAGCCACATGCTTTGTTGACGGACGTCCTGACGCGt                                                                                                                                                                                                                                <  1:57593/249‑1 (MQ=255)              gACCGCTTCTTGTCCTTCCAGGCGGATCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACGCACAACAATATGTGCGAGCGTGCTATCCGTCCGGTCGTTATTGGACGAAAGGCCTGGTTGTGCGCTGGTTCACTGGTGGCCGGGAGCCGCGCTGCTCGTATAATGGGCCTTCTGGGAGCCGCAGTCTGGAGCCATGTGCGTGGGTGACGGACGTTCTGAGGCGTCTg                                                                                                                                                                                                                             <  2:403418/251‑1 (MQ=255)                         gTCCTTCCAGTCGGAGCGAGAGTTAGCCGTGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTTGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCGGGAGCCTCATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGcc                                                                                                                                                                                                                  <  1:28049/251‑1 (MQ=255)                             ttCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTATGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGCGTGGCCgg                                                                                                                                                                                                                >  2:240290/1‑249 (MQ=255)                                      gAGCGAGAGTTAGACGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACAGCGCGGCACAGATAATGAGCCTTCTGGGACCCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCATGACGCGTCTGCCGGCGTGGCCGGGGGAGCGGTTg                                                                                                                                                                                                     >  2:334742/1‑251 (MQ=255)                                                     ggggCAGACTTTCACCGGCGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGGGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGGGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCGCATGCTTGGGTGACGGACGGCCTGACGCGGCTGCCGGAGTGGCCGGAGGAGAGGTTGGGTGATTTACTGCct                                                                                                                                                                                      <  2:183257/251‑1 (MQ=255)                                                                 cACCGACGGAAAGATCGTGCATTTCAAGGCACACATAACAATATCTGCGTGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGGCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCtt                                                                                                                                                                           <  2:285308/250‑1 (MQ=255)                                                                    cGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTAcc                                                                                                                                                                       >  2:289064/1‑251 (MQ=255)                                                                        ggAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTACGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGCGAGGTTGGCTGCGTTACTGCCTCTTGAGGGCTTTACcttct                                                                                                                                                                   >  1:75184/1‑251 (MQ=255)                                                                               tcgtgactttCAAGGCACACACAACAATATCTGCGAGCGTGTTATCCGACCGGTCGTTCTGGGACGAAAGGCATGGTTGTTCGCTGGTTAACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGCGCCACATGCTTGGTTGACGGACGTCCTGCCGCGTCTGCCGGAGTGGCCGGAGGCGAGGTGGGCTGAGTTACTGCCTCTTGAGGTCTTTACATTCTTCGGGt                                                                                                                                                             >  2:326407/8‑250 (MQ=255)                                                                                cGTGCATTTCAAGGCACACTCAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCTCTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCGGAGTTACTGCCACTTGAGGGTTTTACCTTCTTCggggaa                                                                                                                                                           <  2:263728/251‑4 (MQ=255)                                                                                 gTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAg                                                                                                                                                          <  1:203443/251‑1 (MQ=255)                                                                                    cATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACt                                                                                                                                                                                          >  1:23256/1‑216 (MQ=255)                                                                                    cATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACt                                                                                                                                                                                          <  2:23256/216‑1 (MQ=255)                                                                                       ttCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATAc                                                                                                                                                    >  1:203002/1‑251 (MQ=255)                                                                                         cAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGCTAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCt                                                                                                                                                  >  2:374373/1‑251 (MQ=255)                                                                                         cAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGGTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGAGAATTAGCCTTCTGGTAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTTGCCGGAGGAGATGGTGGCTTATTTACTGCCTGTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCt                                                                                                                                                  <  2:56608/251‑1 (MQ=255)                                                                                          aaGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTg                                                                                                                                                 >  1:31142/1‑251 (MQ=255)                                                                                                 acacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGTGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTTCCGCTTGATGGGTTTCCCATCTGCGGGTAAGTGATACCTTCCGTCAg                                                                                                                                          <  2:247884/251‑1 (MQ=255)                                                                                                   acaacaATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAGGTGATACCTGACGTCAggg                                                                                                                                        >  2:12381/1‑250 (MQ=255)                                                                                                     aacaaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGCCGCAAGACCTGGTTGTTCGCTGGTTCACTGATGGCCGGGAACCGCGCGGCACGGCTAATGAGCCTTCGGGGAACCGCAGTCTGGAGCCACATGCTTGGTTTTCGGACGTCCTGACGCGTCTGCCGGCGTCGCCGGCGGAGAGGTTGTCTGAGGTACTGCATCTTGAGGTCTGTACCTTCTTCGTTTAGGTGATACCTGCCGTCATgtgt                                                                                                                                      >  2:117697/1‑251 (MQ=255)                                                                                                      acaaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGtt                                                                                                                                     <  1:250583/251‑1 (MQ=255)                                                                                                               gCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGTTGTTTGTGCAcc                                                                                                                             <  1:173040/250‑1 (MQ=255)                                                                                                                 gAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACATATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCggg                                                                                                                          <  1:229964/251‑1 (MQ=255)                                                                                                                 gAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCggg                                                                                                                          <  1:143900/251‑1 (MQ=255)                                                                                                                       gCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATa                                                                                                                     >  1:316393/1‑250 (MQ=255)                                                                                                                                 cgGTCGTTATGGGACGAAAGCCCTGGTTAATCGCTCGCTCCCTGGTGGCCGGGCACCGCGCGGCCCAGATAATGAGCCTTCTGGGAACCGCCGGCTGGAGCCACCTGCCTGGTTGCCGGCCGTCCTGACGCGTCTGCCGGAGTGACCGGAGGAGAGGATGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGCCAATGATACCTGCCGTAAGGGGTTTGTGCACCGGGACCTAAACTGCAgatg                                                                                                          >  2:402467/1‑248 (MQ=255)                                                                                                                                         aTGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTg                                                                                                   <  1:34100/250‑1 (MQ=255)                                                                                                                                                 aaaGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAAACGCAGTCTGGAGCCACATGCTTGGTTGACGGCCGTCCTGACGCGTCTGCCGGAGTGGACGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTg                                                                                           >  1:356736/1‑250 (MQ=255)                                                                                                                                                  aaGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGac                                                                                         <  1:374373/251‑1 (MQ=255)                                                                                                                                                             gttCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAAt                                                                               <  1:119983/250‑1 (MQ=255)                                                                                                                                                                cGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGGGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGc                                                                           <  1:397923/251‑1 (MQ=255)                                                                                                                                                                 gCTGTTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGGCTGGAGCCACATGCTTGGTTGACGGTCGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACATGAAATGAGCt                                                                          <  2:316393/251‑1 (MQ=255)                                                                                                                                                                            ggtggCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTaa                                                                <  1:103227/250‑1 (MQ=255)                                                                                                                                                                             gtggCCGGGAACCGCGCGGCACAGATAATGAGCCTGCTGGGAACCGTAGTCTGGAGCCACATGCTTGGTTGGCGGACGTCCTGACGCGTCGGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAac                                                              <  2:255786/251‑1 (MQ=255)                                                                                                                                                                                    gtaaCCGCGCGGCACAGATAATGAGCCGTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTGGGCTGAGTTACGGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGGTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCggg                                                       <  2:175083/249‑1 (MQ=255)                                                                                                                                                                                         cgcgcgGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGAc                                                  <  1:289064/251‑1 (MQ=255)                                                                                                                                                                                             cgGCACAGATAATGAGCCTTCTGGGAACCGCAGTATGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATg                                              >  2:19588/1‑251 (MQ=255)                                                                                                                                                                                              ggCACAGATAAGGAGCATGATGGGAACCGCAGTCTGGAGCCACATGCTTGGTGGACGGACGTCCTGACGCGTCGTCCGGAGTGGGCGGAGGAGAGGTTGGTTGAGTTACTGCCTTTTGAGGGCTTTACCTTCTTCGGGGCAGTGATTCCTGCCGTAAGGTGGTTGTGCACCGGGCCATAACTTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATgg                                             <  1:117697/251‑1 (MQ=255)                                                                                                                                                                                               gCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATggg                                            >  1:119123/1‑251 (MQ=255)                                                                                                                                                                                                  aaGATAATGAGCCTTCTGTGACCCGAAGTCTGGAGCCACAGGCTTGGTTGACGGACGTCATGACGCTTTTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTGAATGCCTCTTGAGGGCTTTACCTTCTTCTGGTAAGTGATACCTGCCGTCAGTGGTTTTTGCTCCGGGCCATAACCTGCAGTTGGGTATTGATCGCCTGAAAACAGGAAATGAGCTAGAGCAGTAAACACACGGGCTGACCATGGGTCa                                         <  2:271674/250‑1 (MQ=255)                                                                                                                                                                                                   aGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTc                                          >  1:318710/1‑249 (MQ=255)                                                                                                                                                                                                        aTGAGCCTTCTGGTAACCGCAGTCGGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCTTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTTAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAatat                                   <  2:165571/251‑1 (MQ=255)                                                                                                                                                                                                        aTGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAatat                                   >  2:339113/1‑251 (MQ=255)                                                                                                                                                                                                           aGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGa                                                                                                                                                       <  2:360929/132‑1 (MQ=255)                                                                                                                                                                                                           aGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGa                                                                                                                                                       >  1:360929/1‑132 (MQ=255)                                                                                                                                                                                                                tCTGGGAACCGCAGGCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTa                             <  2:350972/249‑1 (MQ=255)                                                                                                                                                                                                                    ggAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGTGCTGACAATGGGTCAGAATATGAGTTAGCCCgg                       <  1:19491/251‑1 (MQ=255)                                                                                                                                                                                                                               caggAGCCACACGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCGCTTGAGGGCTTTCCCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTg            <  2:327492/249‑1 (MQ=255)                                                                                                                                                                                                                                  gAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGCGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACCATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTgggg         >  2:255430/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGGGTGGCCGGAGGAGAGGTTGGCTGGGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGtta  <  2:107112/250‑3 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   TTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAC  >  NC_000913/4501260‑4501742

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