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.3%	158.7 / NA	54	L81L (CTG→TTG)	yjgZ	uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/34);  minor base G (0/2);  total (18/36)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTG  >  NC_000913/4501259‑4501732                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          tttgggaccgcctGGACCGCTTCTTGGCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCTGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAAGATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTTGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTgacggac                                                                                                                                                                                                                                   <  2:151005/240‑1 (MQ=255)    gTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGCCGGACAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGACTGGAGCCACAGTCTTGGTTGACGGACGTcc                                                                                                                                                                                                                               >  2:41600/1‑251 (MQ=255)         ctcctGGACCGCTTGTTTTCCTTCCAGTCGGAGCGAGCGTTAGCCGGGGCGGACTTTCACCGACGAAAAGATCGCGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCCCTGGTGGCCGGGAACCGCGCGGCACCGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCAACCTGCTTGGTTGCCGGCGGTCCTGAcg                                                                                                                                                                                                                          >  1:75830/1‑251 (MQ=255)           cctGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGAcg                                                                                                                                                                                                                          <  1:133939/249‑1 (MQ=255)                               tCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggag                                                                                                                                                                                                     <  1:13120/250‑1 (MQ=255)                                ccAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggaggg                                                                                                                                                                                                   >  2:375991/1‑250 (MQ=255)                                      ggAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCTTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTGCTGGGAACCGCAGTCTGGAGCCACATGCTGGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGtt                                                                                                                                                                                             <  1:20494/251‑1 (MQ=255)                                        aGCGAGAGTTAACGGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACATAAATCGCTGCGAGCGTGCTATCCGTCCGGGCGTTATGGGACGAAAGGCCTGGTGGT‑TCGCGGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGAGAAGGAGCCTGCGGGGAACCGCAGTCTGGGGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTgg                                                                                                                                                                                           <  2:191431/251‑1 (MQ=255)                                             gagTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAg                                                                                                                                                                                      <  1:67988/251‑1 (MQ=255)                                                 tAGCCGGTGCGTACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGTGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTa                                                                                                                                                                                   <  2:210926/250‑1 (MQ=255)                                                          cGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTg                                                                                                                                                                         >  1:391222/1‑251 (MQ=255)                                                                 cgaCCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAAGATCTGCGAGCGGGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCt                                                                                                                                                                   <  1:413148/248‑1 (MQ=255)                                                                  cACCGACGGAAAGGTCGTGCATTTCAAGGCACACACAACAATATCTGCGTGCGTGCTATCCGTCCGGTCGTGATGGGACGTAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCGGGAGCTACAGGCTTGGTTGTCGGACGTTGTGTCTCGTCTGCCGGAGTTGCCGGAGGAGAGTTTTGCTCAGTTACTGtctcgtggtggctt                                                                                                                                                                  <  2:222177/250‑11 (MQ=255)                                                                   aCCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTa                                                                                                                                                                <  1:375991/251‑1 (MQ=255)                                                                        cGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACcttc                                                                                                                                                           <  2:391222/251‑1 (MQ=255)                                                                          gAAAGATCGTGCATTTCACGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGCGGAGGGGTTGGCTGAGTTACTGCCTCTTGAGGGCt                                                                                                                                                                   >  2:136297/1‑241 (MQ=255)                                                                          gAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCt                                                                                                                                                                   <  1:136297/241‑1 (MQ=255)                                                                              gATCGTGCATTTCAAGGCGCACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCggg                                                                                                                                                     <  1:210697/251‑1 (MQ=255)                                                                                    gCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTg                                                                                                                                               >  1:424470/1‑251 (MQ=255)                                                                                    gCATTCCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTAGGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGTCTGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCGGGAGCCACATGCTTGGTTGACGGACGTCCGGACGCGGCTGCCGGATTGTCCGGATGAGAGGTTGGCTGAGTTACTGCCTCCTGAGGGCGTTAACTGCTTCgtttaagt                                                                                                                                                <  2:410915/250‑8 (MQ=255)                                                                                         tCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGATCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATAcc                                                                                                                                          >  2:420900/1‑251 (MQ=255)                                                                                          cAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCt                                                                                                                                         <  1:162399/251‑1 (MQ=255)                                                                                            aGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTTAGGGCTTTACCTTCTTCGGGTAAGTGATTCCTGc                                                                                                                                       <  1:141292/251‑1 (MQ=255)                                                                                              ccactcacaACAATATCTTCGAGCTCGCTATCCTTCCTGTCGATATGGGACGAAAGGCCTTGTTGT‑TCGCTGGTTCACTGGTTGCCGGGCACCGCGCTGCCCAGATAATGAGCCTTCTGGGAACCGCAGTCTTGAGCCACTTGCTTGTTTTACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGATAGGTTGGCTGAGTTACTGCCTCTTGAGTGCGTTACCTTCTTCTGGTAAGTGATACCttcct                                                                                                                                     <  1:125837/246‑5 (MQ=255)                                                                                               cacacacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑Gt                                                                                                                                   <  1:13318/251‑1 (MQ=255)                                                                                                   cacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAgg                                                                                                                               >  1:168057/1‑251 (MQ=255)                                                                                                          atatCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCGTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTtgt                                                                                                                        <  1:314149/251‑1 (MQ=255)                                                                                                                 cGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAACGGCCTGGCTGT‑TGGCTGGTTAACTGCTGGTCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGCGTGGCCGGAGGCGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTCAGTGATACC‑GCCAGTCAGGAGTTTGTGCAGCgg                                                                                                                 >  2:151288/1‑251 (MQ=255)                                                                                                                  gAGCGTGCTAGCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTGCTTCGGGTAAGTGATACCGGCC‑TTCAGGTGTTTGTGCACCggg                                                                                                                <  2:28087/251‑1 (MQ=255)                                                                                                                         cTATCCGTCCGGTCGTTATGGGACGCAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGCGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAAc                                                                                                         >  2:154125/1‑251 (MQ=255)                                                                                                                           aTCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCt                                                                                                       >  1:416745/1‑251 (MQ=255)                                                                                                                                  cgtgcGTTATGGGGCGAAAGGCCTGGT‑GTGTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATTAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTTACGCGTCTGCCGGAGTGGCCGGAGGAGATGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATTCCTGCC‑GTCAGGGGTTTGTGCACCGGGCCATAACCTGCAGTTg                                                                                                <  2:255/247‑1 (MQ=255)                                                                                                                                  cgGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTg                                                                                                <  1:74522/251‑1 (MQ=255)                                                                                                                                      cGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGa                                                                                             >  2:185925/1‑250 (MQ=255)                                                                                                                                             ggACGAAAGGCCTGGTTGT‑GCGCGGGTTCACGGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATGGAACg                                                                                     <  2:229148/251‑1 (MQ=255)                                                                                                                                                cGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCt                                                                                  <  2:143747/251‑1 (MQ=255)                                                                                                                                                cGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCt                                                                                  >  1:355449/1‑251 (MQ=255)                                                                                                                                                    aGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGaca                                                                              <  1:225096/251‑1 (MQ=255)                                                                                                                                                                 tCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGa                                                                   <  1:330350/250‑1 (MQ=255)                                                                                                                                                                  cGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGc                                                                 <  1:350017/251‑1 (MQ=255)                                                                                                                                                                       gttaaCTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACTCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTagag                                                            <  1:77573/247‑1 (MQ=255)                                                                                                                                                                               gtggCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTGCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAac                                                    <  1:287339/251‑1 (MQ=255)                                                                                                                                                                                     ggtaaCCGCGCGGCACAGATCATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGGTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTTAGTTACTGCCTCTTGAGGGCTTTACCTTCTGCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGATCAGTAAACAcgcga                                              <  2:310413/248‑2 (MQ=255)                                                                                                                                                                                       gAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGc                                            >  2:357909/1‑251 (MQ=255)                                                                                                                                                                                           cgcgcgGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGAc                                        <  1:154125/251‑1 (MQ=255)                                                                                                                                                                                           cgcgcgGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGGGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGAc                                        >  2:231367/1‑251 (MQ=255)                                                                                                                                                                                               cgGCACACATAATGAGCCTTCTGCGAACCGCAGTCTGGAGCCACATAATTGGTTGACGGACGCCCTGACGCGTAAGCCGGAGTGGCCGGAGGCGAGGCTGGCTGAGTTTCTGCTTCTTGAGGGCTTTACCTTCTCCGGGAAAGCGATACCTGCA‑GTCCGCAGTTTGCGCACCGGGCCATAACCTGCAGTTGGGAATTGACCGCCTACCAACAGGAAATGAGCTAGCACAGTACACACGCGGGcggccattg                                    >  2:428062/1‑243 (MQ=255)                                                                                                                                                                                                       aTAAGGAGCCTTCTGGGAACCGAAGTCTGGAGCCACAGGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTGACTGCCTCTTGAGGGATTTCCCTTCTTCGGGTAAGTGATACGTGCC‑GTCAGGGGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGCGCTGACAATGGGTCTGaa                            <  2:424470/251‑1 (MQ=255)                                                                                                                                                                                                          aTGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTACTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAatat                         >  1:210344/1‑251 (MQ=255)                                                                                                                                                                                                               ccTTCTGGGACCCGCAGTTTGGAGCCACATGCTGGGTTTACGGACGTCCTTACGCGGCTGCCGGAGTGGCCGGAGGAGAGGTTTGCTGAGTTACTGCCTCTTTAGGGCTTTACCTTCTTCTGGTAAGTGATACCTGCC‑GTCAGGTGTTTTTGCTCCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACTGGAAATGAGCTAGAGCACTAAACACGCGGGGTTACAATGGGTCAGAATATGAGtg                    <  2:101781/251‑2 (MQ=255)                                                                                                                                                                                                                         aCCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGt          <  2:416745/251‑1 (MQ=255)                                                                                                                                                                                                                            gCAGTCTGGAGCCACATGCTGGGGTGACGGCCGTCCGGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGATGGCTTTACCTTCGTCGGGTAAGTGATACCTGTC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTGAGCCCGGTGTTAt       <  2:186587/251‑1 (MQ=255)                                                                                                                                                                                                                             cAGTCGGGAGCCACAGGCTTGGTGGACGGACGGCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGTCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTAtt      <  1:45201/251‑1 (MQ=255)                                                                                                                                                                                                                                 cTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTg  >  1:375186/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          TTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGT‑TCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCC‑GTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTG  >  NC_000913/4501259‑4501732

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