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	92.7%	178.4 / NA	68	L81L (CTG→TTG)	yjgZ	uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/35);  minor base G (1/4);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTT  >  NC_000913/4501259‑4501744                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     gtcgTTACCCCCCGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTgacggac                                                                                                                                                                                                                                              <  2:327324/247‑1 (MQ=255)       ccctcctGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGa                                                                                                                                                                                                                                       <  2:507902/250‑1 (MQ=255)        ccttctGGACCGTTTCTTGGCCTTCGAGTCGGAGAGTGATGTAGCCGGGGCTGTCTTTCACCGACTGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGTAACCTCGCGGCATAGATAATGAGCCTTCTGGTAACCGCAGTCTGGAGCCTCATGCTTGGTTGACTGACGTCCTGAc                                                                                                                                                                                                                                      <  2:364593/247‑1 (MQ=255)                aCCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGc                                                                                                                                                                                                                              <  1:541825/251‑1 (MQ=255)                          gTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGcc                                                                                                                                                                                                                    >  1:332271/1‑251 (MQ=255)                                ccAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggagga                                                                                                                                                                                                              <  1:90515/251‑1 (MQ=255)                                     cGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGt                                                                                                                                                                                                         <  1:201817/251‑1 (MQ=255)                                       gAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTg                                                                                                                                                                                                       <  1:365693/251‑1 (MQ=255)                                       gAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGGGGCCgggggcggggttg                                                                                                                                                                                                       >  2:406570/1‑240 (MQ=255)                                       gAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGGTGACGGACGTCCTGACGCGTCGGCCGGAGTGGCCGGAGGAGAGGTTg                                                                                                                                                                                                       <  1:190494/251‑1 (MQ=255)                                        ggCGAGAGTGAGCCGGGTCTGACTTTCACCGACGGAAAGCGCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGGTATGGGACGAAAGGCCTGGTTGGTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTgg                                                                                                                                                                                                      <  2:79130/250‑1 (MQ=255)                                            agtggtAGCCGGGGGGGACTTTCCCCGACAGAAAGTTCGTGCACTTCATGGCTCGTACATCAATATCTTCGATCGTGCTATCCTTCCGGTCGTTTTGGGACGAAATGCTCGGTTGGTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATACTGATCCTTCTGTGAACCGCAGTCTGGAGCCACAGGCTTGGTTGAGTGACTTCCTGACGCGGCTGCCGGAGTGGCCGGAGGAGAGTTTg                                                                                                                                                                                                       <  2:62912/241‑1 (MQ=255)                                                      ggggCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCct                                                                                                                                                                                        <  2:332271/251‑1 (MQ=255)                                                           tgACTTTCACCGGCGGAAAGAGCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGACGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAACCTTCTGGGAACCGCAGGCTGGAGCCACATGCTGGGTTGAAGGACGACCTCACGCGTCTGCAGGAGTGGCCGGGGGTGAGGGTTGCTGAGATGCTGCCTCTTgg                                                                                                                                                                                   >  2:378622/2‑250 (MQ=255)                                                                  cACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCttt                                                                                                                                                                            >  1:249756/1‑251 (MQ=255)                                                                  cACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCtt                                                                                                                                                                             <  1:406570/250‑1 (MQ=255)                                                                    ccGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTAc                                                                                                                                                                          >  1:344138/1‑251 (MQ=255)                                                                            aaGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCg                                                                                                                                                                  <  1:372846/251‑1 (MQ=255)                                                                              gatggTGCATTGCAAGGCCCACTCAACAATAGCTGCTAGCGTGCTAGCCGGTCGGTCTGTATGGTTCGAAAGGCCTGGTTGGTCGCTGGTTCCCTGGGTGCCGGGTACCGCGCGTCTCAGATAATGAGACTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGGTGACGGACGTCCTGACGCGTCTGTCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCCTGAGGGCTTTTCCTTCTTAggg                                                                                                                                                                <  2:185197/247‑1 (MQ=255)                                                                                  gTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCTCAGTCTGGAGCCACATGCTTGGTTGCCGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGGGAGGTTGGCTGAGTTACTGCCTCTTGTGGGCTTTACCTTCTTCGGGTagg                                                                                                                                                            >  2:142775/1‑249 (MQ=255)                                                                                      aTTTCAAGGCACACACAACAATATCTGCGAGCGTGCTCTCCGTACGGTCGTTATGGGACGAAAGGCCTGGTTGTTAGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGGGGTTGGCTGGGTTACTGCCTCTTGGGGGCTTTACCTTCTTCGGGTaaagttc                                                                                                                                                        >  2:290277/1‑246 (MQ=255)                                                                                      aTTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGCGGTTGGCTGAGTTACTGCctct                                                                                                                                                                                      >  2:243131/1‑221 (MQ=255)                                                                                      aTTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCctct                                                                                                                                                                                      <  1:243131/221‑1 (MQ=255)                                                                                      aTTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTaa                                                                                                                                                             <  1:290277/246‑1 (MQ=255)                                                                                            aGGCACTCACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGc                                                                                                                                                  <  2:166832/251‑1 (MQ=255)                                                                                            aGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGc                                                                                                                                                  >  1:448963/1‑251 (MQ=255)                                                                                              gCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGATGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCg                                                                                                                                                <  1:4447/251‑1 (MQ=255)                                                                                               cacacacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGTCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGt                                                                                                                                               >  1:10288/1‑251 (MQ=255)                                                                                                 cacacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCa                                                                                                                                             <  1:239087/251‑1 (MQ=255)                                                                                                      aacaaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGgtgt                                                                                                                                        >  1:501685/1‑251 (MQ=255)                                                                                                       acaaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGtt                                                                                                                                       >  1:22617/1‑251 (MQ=255)                                                                                                                gtgagggtgttATCCGTCCGGTCGTTATGGGGCGAAAGGCTGGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCATAGATAATGAGCCTTCTGGTAACTGCAGTCTGGAGCCACATGCTTGGTTGCCGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACGGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATTCCTGCCGTCAGGTGTTTGTGCAcc                                                                                                                               <  2:447797/240‑1 (MQ=255)                                                                                                                  gAGCGTGCTATCAGTCCGGTCGTTATGGGACGAAAGGACTGGTTGTTCGCTGGTTCACAGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTGTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCATGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTATGTGATACCTGACGTCAGGTGTTTGTGCACCggg                                                                                                                            >  2:252045/1‑251 (MQ=255)                                                                                                                   ggCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATATCTGCCGTCAGGTGTTTGTGCACCGGGc                                                                                                                           <  1:148214/250‑1 (MQ=255)                                                                                                                    gCGTGTTATCCGTTCGTTAGCTATGGGACGAAAGGTCTAGTTTTTCGCGGTTTCTCTGGTGGACGGGAACCGCTCGGCACTGATAATGGGCCTTCTGGGAGCCGCAGTCTTGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGTAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGTGCTTTACCTGCTTCGGGTAAGGGATACCTGCCGTCAGGTGTTtgttcttcgggc                                                                                                                           <  2:333554/250‑10 (MQ=255)                                                                                                                      gTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCAt                                                                                                                        >  1:7310/1‑251 (MQ=255)                                                                                                                        gCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACCGGGGGCCGGGAGCCGCGCGGCACAGATAATGAGCCTTCTGGGAAACGCAGTCTGGAGCCACATGCTTGGGTGCCGGACGTCCTGACGCGTCTGCCGCAGTGGCCGGAGGAGAGGTTGGCTGAGTTTCTGCCTCTTGAGGTCTTTACCTTCTTCGGGTAAGTGACTCCGGCCGTCAGGTGTTTGTGCTCCGGGCCATa                                                                                                                       <  1:199927/250‑1 (MQ=255)                                                                                                                            tccgtccgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTg                                                                                                                                                                                    <  2:209058/185‑1 (MQ=255)                                                                                                                            tccgtccgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTg                                                                                                                                                                                    >  1:209058/1‑185 (MQ=255)                                                                                                                               gtcaggtagTTATGAGACGCAAGACCTGGTTGTTCGCTGGTTCACTGTAGGCCTGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCGGTCTGGAGCCACATACGTGGGTGACGGACGTCCTGACGCGGCTGCCGGCGGGGCCGGCGGAGAGGGTGGCTGATTTACTGCCTCTTTAGGGCTTTACCTTCTTCGGGGACGAGAGACCTGACTTCAGGTGTGTATGCGCCGGGCCCTAACATGc                                                                                                                 >  2:64/9‑249 (MQ=255)                                                                                                                                  cgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTg                                                                                                            >  1:525848/1‑251 (MQ=255)                                                                                                                                   ggTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACCGGTGGCCGGGAACCGCGCGGCACCGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCAAATGCTTGGTTGACGGCCGACCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGATGAGTTTATGCCTCTTGAGGGCTTTACCTTCTTCGGGTGAGTGATAACTGGCGTCAGGAGTTTTTGCAAAGAGACATAACCTGCAGTTgg                                                                                                           >  2:101805/1‑251 (MQ=255)                                                                                                                                       gTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAAt                                                                                                       <  1:142775/251‑1 (MQ=255)                                                                                                                                         tATGGGATGAAAGGGCTGGTTTTTCGCTGGTTCCCTGGGGGCCGGGACCCGCGCGGCACAGTTAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTTGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTGGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCTCCGGGCCATAACCTGCAGTTGGGAAtt                                                                                                      <  1:530402/250‑1 (MQ=255)                                                                                                                                          aTGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTg                                                                                                     <  1:218641/250‑1 (MQ=255)                                                                                                                                             ggACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGAGCGGCACAGATAATGAGCATTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGTTGAGTTACTGCCTCTTGAGGGCTTTACATTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGGGCACCGGGCCATAACCAGCAGTTGGGAATTGAACg                                                                                                 <  1:159495/251‑1 (MQ=255)                                                                                                                                              gACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGc                                                                                                >  1:512019/1‑251 (MQ=255)                                                                                                                                                cgaatggccagtttgttCGCTGGTTCGCTGGTGGCCGGGAACCGCGCGGCACAGATAAGGATCCTTCTGGGAACCGCAGTCTGGAGCCACATGCGTGGGTGACGGACGTCCTGACGCGGCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATCCCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTTGGAATTGAACGCCt                                                                                              <  2:177666/239‑1 (MQ=255)                                                                                                                                                 gAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTg                                                                                             >  1:259815/1‑251 (MQ=255)                                                                                                                                                  accgacgcccttacggtcTGGTTCACAGTTGGCCGGGAACCGCGCGGGACAGATAAGGGGCCATCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGGCTGCCGGGGTGGCCGGAGGAGAGGTTGGCTGAGTTTCTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGa                                                                                            <  2:533844/234‑1 (MQ=255)                                                                                                                                                    aGGCCTGGTTGTTCGCTGGTTCACTGGTGGGCGGGAACCGCGCGGCGCAGATACTGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGTCGCGTCTTCCGGAGTGGCCGGAGGAGGTGTTGGCTGAGTTACTGCCTCTTGAGTGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCGGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGaca                                                                                          <  1:293059/251‑1 (MQ=255)                                                                                                                                                                     gtttCACTGGTGGCCGGGAACCGCGCGGCAAAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCGTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTaga                                                                         <  2:109502/249‑1 (MQ=255)                                                                                                                                                                     ggTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACCTGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATACCCTGCAGTTGGGAATTGAACGCCTGacgatcggaagagcccccgt                                                                         >  1:533844/1‑233 (MQ=255)                                                                                                                                                                                      gAACCGCGCGGCACAGATAATGAGCCTTCTGGGAAACGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGATACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACCACAGGAAATGAGCTAGAGCAGTAACCACGCGGGc                                                        >  1:526544/1‑251 (MQ=255)                                                                                                                                                                                        acctagcgGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCCCATGCTGGGTTGACGGACGTCCTGACGCGTCTGCCGGATTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAAGTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTg                                                      <  2:179072/246‑1 (MQ=255)                                                                                                                                                                                        aCCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGG‑CTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTg                                                      >  1:288658/1‑250 (MQ=255)                                                                                                                                                                                               ggCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCCCATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAAt                                                 <  2:279600/249‑1 (MQ=255)                                                                                                                                                                                               ggCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAAt                                                 >  1:279600/1‑249 (MQ=255)                                                                                                                                                                                                     gATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGa                                         <  1:334381/251‑1 (MQ=255)                                                                                                                                                                                                         aTGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAatat                                     <  1:542609/251‑1 (MQ=255)                                                                                                                                                                                                          tGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATg                                    >  1:572165/1‑251 (MQ=255)                                                                                                                                                                                                                   tGGGATCCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCGGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCGCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAg                              <  2:122129/248‑1 (MQ=255)                                                                                                                                                                                                                                 tGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTTCCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTgg             >  1:170416/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGTAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGTGTTATTCCTgg             <  1:113825/250‑1 (MQ=255)                                                                                                                                                                                                                                    aGCCACATGCTGGGTTGACGGACGTCCTGACGCGTCTGCAGGAGTGACCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTATTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAAAGCCTGACAACAGGCAATGAACTAGAGCAGTAAACACGCGGGCTGACAGTGGGTCAGAATATGAGATAGCCCGGCGGTATTCCTggggg          >  2:141959/1‑250 (MQ=255)                                                                                                                                                                                                                                      ccACATGCTGGGTTGACGGACGTCCTGACGCGGCTGCCGGATGGGCCGGAGGAGAGGTTGGCTGTGTTACTGCCTCTTGAGTGCTTTACCTGCTTCGGGTAAGTGATACCTGCCGGCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGTGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCa        <  2:448963/251‑1 (MQ=255)                                                                                                                                                                                                                                           tGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATCTGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACt   >  1:91347/1‑251 (MQ=255)                                                                                                                                                                                                                                             cTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGACGTTATTCCTGGGGTCAGTACtt  >  2:465387/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     TTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTT  >  NC_000913/4501259‑4501744

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