Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,221	T→C	L193L (CTA→CTG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,221	0	T	C	88.1%	101.9 / NA	42	L193L (CTA→CTG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (16/21);  minor base A (5/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 4.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTT  >  NC_000913/4306931‑4307499                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                       tGAGCGCCATTAAAGTCAGAAGATTATGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAgcgc                                                                                                                                                                                                                                                                                   <  2:373626/296‑1 (MQ=255)   aGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTcc                                                                                                                                                                                                                                                                             <  2:312056/300‑1 (MQ=255)       ccATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCCATGGCCAATGCCGCGTTTTCAAGCAGACTCTGGAAGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAGTTTCGGGGCTGTTGTTCGTACCAGCGTCTAAGCGCCATTCAaga                                                                                                                                                                                                                                                                          >  2:220239/1‑299 (MQ=255)       ccATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCaga                                                                                                                                                                                                                                                                          >  2:471241/1‑299 (MQ=255)          ttAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCaa                                                                                                                                                                                                                                                                      <  2:689339/300‑1 (MQ=255)                                                       ggCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGt                                                                                                                                                                                                                          <  2:718432/299‑1 (MQ=255)                                                       ggCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGt                                                                                                                                                                                                                          <  1:94916/299‑1 (MQ=255)                                                       ggCAGGGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGt                                                                                                                                                                                                                          <  2:718535/299‑1 (MQ=255)                                                        gCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCTGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCa                                                                                                                                                                                                                        >  1:38495/1‑300 (MQ=255)                                                        gCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCa                                                                                                                                                                                                                        >  1:452759/1‑300 (MQ=255)                                                            tGGAAGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGt                                                                                                                                                                                                                          <  2:164520/294‑1 (MQ=255)                                                             ggACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATAtct                                                                                                                                                                                                                   <  1:183710/300‑1 (MQ=255)                                                             ggACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGAGTGGGAAGGCGGGCATATTCATCACGCAACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCAGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACACGCAGTGTTGGGTCatattt                                                                                                                                                                                                                   >  2:103565/1‑298 (MQ=255)                                                                    gcAGGTACTTTTGGTTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                            <  1:501416/300‑1 (MQ=255)                                                                    gcAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCTGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                            <  1:103525/300‑1 (MQ=255)                                                                      aGGTACTTTTGTTTCATGGCAACGAGAGTCTCGCTTTGGAAGGCGGGCATATTTATCACGCCACTGCCTAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAATGGTCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                                                                           <  2:214787/299‑1 (MQ=255)                                                                      aGGTACTTTTGGTTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                                                                           <  1:99864/299‑1 (MQ=255)                                                                          aCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                                                                                                                                                                      <  1:745470/300‑1 (MQ=255)                                                                           cTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                                                                                                                                                                      >  1:215905/1‑299 (MQ=255)                                                                              ttGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                                                                                  <  1:238970/300‑1 (MQ=255)                                                                                                         tGGGAAGGCGGGCATATTCATCACGCCACCGCCCAGGATCACCGCATCGGGAGCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAATGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAcc                                                                                                                                                                       <  1:283850/300‑1 (MQ=255)                                                                                                         tGGGAAGGCGGGCATATCCAGCACGCCACCGCCCACGATCACCACATCGGGATCGAACAGAGTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTGCAAGCAGACTCTGGACGGAAGGTGCGTTGTCCGTATGGACGAAAAGATCGCGCAATGGGTAATTTCGGTGCTGTTGTTCGTACCAGCATCTCAGCGCCAGTCCAGAGCAATTGGGTTCCAGGCACCCAGAATTGCCACACGCGCAGTGTTGGGt                                                                                                                                                                                                                          <  2:164363/249‑1 (MQ=255)                                                                                                                                 gCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                                                                                                <  2:215905/299‑1 (MQ=255)                                                                                                                                    aCCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGTGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTGCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                            <  1:656139/300‑1 (MQ=255)                                                                                                                                       gCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAGTTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCACAACGTGTGCCCCCGTCCACGGCGCACCGTTCATCCACACTGc                                                                                                                                         >  2:517488/1‑300 (MQ=255)                                                                                                                                       gCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCCGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAAGGCGAACCGTTCACCCACACt                                                                                                                                           >  2:517273/1‑298 (MQ=255)                                                                                                                                        cccAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGc                                                                                                                                         >  1:352026/1‑299 (MQ=255)                                                                                                                                                                    gATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCCGGGGGATATGACCCAGTTCGCATGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACACCATCCCCGTACCGAGATAGgcc                                                                                                             >  2:491562/1‑299 (MQ=255)                                                                                                                                                                    gATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTAAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCAACACTGCGAACCCCATCCCCGTACCGAGATAGgcc                                                                                                             >  2:436490/1‑299 (MQ=255)                                                                                                                                                                                 ggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                               <  1:371487/300‑1 (MQ=255)                                                                                                                                                                                            gTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGc                                                                                    >  1:389252/1‑300 (MQ=255)                                                                                                                                                                                            gTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCAAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCAGCCAGAACCAGTTGTTGCGTAAGGc                                                                                    >  1:722821/1‑300 (MQ=255)                                                                                                                                                                                            gTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTAAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCAAACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGc                                                                                    >  1:386891/1‑300 (MQ=255)                                                                                                                                                                                                          aGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAGAAGATCGCGCAATGGGTGATTTCGGGGCTGTTGTTCGTACCAGCGTCTAAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCAAACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCAAACCCCATCCCCGTACCGAGAt                                                                                                                  >  2:700840/1‑256 (MQ=255)                                                                                                                                                                                                                ctctGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccaa                                                                 >  2:679390/1‑297 (MQ=255)                                                                                                                                                                                                                      aCGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCTAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGt                                                          <  1:217743/300‑1 (MQ=255)                                                                                                                                                                                                                                   tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGt                                             <  1:654698/300‑1 (MQ=255)                                                                                                                                                                                                                                                                gggTAGTTTCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCGGGATTACCAAACGCGAAGTGTTGGGTCATATCTACCAGGGGGATATGACCCCGTTCGCCTGCCACACCATGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCACCCCATCCCCGTACCGATATAGGCCGCCAGACCCAGCTGTTGCGTAAGGCGGTTTTCTACTCCGTACCAGGAGCGCTGCAGGTTAACGTCGAGGGAAAACTCAACCGGACAATTCAGc                >  1:666133/1‑300 (MQ=255)                                                                                                                                                                                                                                                                      tttCGGGGCTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTAAGTCCAAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCCACCGGACAATTCAGCGTAttt          >  2:722839/1‑300 (MQ=255)                                                                                                                                                                                                                                                                              cTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCtt  >  1:565385/1‑300 (MQ=255)                                                                                                                                                                                                                                                                              cTGTTGTTCGTACCAGCGTCTCAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCAAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCtt  >  1:186682/1‑300 (MQ=255)                                                                                                                                                                                                                                                                              cTGTTGTTCGTACCAGCGTCTAAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCtt  >  1:565396/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                       TGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTT  >  NC_000913/4306931‑4307499

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