Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,151	G→A	P217S (CCT→TCT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,151	0	G	A	95.6%	196.8 / ‑3.8	69	P217S (CCT→TCT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (34/31);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGC  >  NC_000913/4306904‑4307394                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       cGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGc                                                                                                                                                                                                                                                     >  1:264410/1‑251 (MQ=255)     gATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGttt                                                                                                                                                                                                                                                 >  1:294462/1‑251 (MQ=255)      aTGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGCGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGtttt                                                                                                                                                                                                                                                >  1:458413/1‑251 (MQ=255)              atatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACTGGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGTAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTGTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATgg                                                                                                                                                                                                                                        <  2:538039/251‑1 (MQ=255)               tatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGCTCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATgg                                                                                                                                                                                                                                        >  1:231926/1‑250 (MQ=255)                atTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGa                                                                                                                                                                                                                                       <  1:612088/250‑1 (MQ=255)                   gctgcACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGa                                                                                                                                                                                                                                    >  1:540441/1‑250 (MQ=255)                         cccTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGaaa                                                                                                                                                                                                                                  >  1:384067/1‑246 (MQ=255)                         cccTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATc                                                                                                                                                                                                                             <  2:175691/251‑1 (MQ=255)                           cTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCCGGATCCCCGCCTCGGGATCGGACAGATTAATGCTGGGGGCAATGGCCCGTGCCCGGTTTTC‑AA‑GCCGACTCTGGAAGAAAGAGGCGGTTTCAGCCTGGCCGaaaaaaacgc                                                                                                                                                                                                                           >  2:586932/1‑245 (MQ=255)                            tGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGCGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGAGTTTTC‑AA‑GCCGAATCTGGACGAAAGAGGAGTTTTCCGCATGGACGAAAAGATcgcg                                                                                                                                                                                                                          >  1:611747/1‑251 (MQ=255)                             gAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑CA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGCCGAAAAggactcgc                                                                                                                                                                                                                         >  2:505445/1‑244 (MQ=255)                             gAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCgcgc                                                                                                                                                                                                                         >  1:565467/1‑251 (MQ=255)                                   cATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATggg                                                                                                                                                                                                                   <  1:599845/251‑1 (MQ=255)                                   cATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATggg                                                                                                                                                                                                                   >  1:450372/1‑251 (MQ=255)                                        aaGTCAGAAGAGGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAAtt                                                                                                                                                                                                              <  1:527521/251‑1 (MQ=255)                                            cAGAAGATGAGGGAGCAATAAAGCGCACGACCGGATGCGGCAGTGGACGGCGCAGGTACTTTGGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTGC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCg                                                                                                                                                                                                           <  2:334132/250‑1 (MQ=255)                                               aagaTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGaaa                                                                                                                                                                                                                                  <  1:403884/224‑1 (MQ=255)                                                agaTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCt                                                                                                                                                                                                      <  1:534287/251‑1 (MQ=255)                                                        cagcAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCg                                                                                                                                                                                              <  1:350224/249‑1 (MQ=255)                                                         agcAATAAAGCGCACGACCTGATGCGGTAGGTGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGt                                                                                                                                                                                             <  2:450372/250‑1 (MQ=255)                                                            aaTAAAGCGCAGGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGGCATGGCAACGAGAGTCTCGCGGGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTAcc                                                                                                                                                                                          <  2:264410/251‑1 (MQ=255)                                                               aaaGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGAAATTTCGGGGCTGTTGCTCGTACCAGc                                                                                                                                                                                       >  1:432624/1‑251 (MQ=255)                                                                   cgcACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCCATAGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCt                                                                                                                                                                                   >  1:29/1‑251 (MQ=255)                                                                      aCGACCTGATGCGGGAGTGGACGGCGTAGGTTCTTTGGGGTCATGTCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATTACGCCGCCGCCCAGGATCACCGCGTCGGGATCGATCAGATTAATGCTGGGGGCAAGGGCCCGTGCCGCGTTTTC‑AA‑GCATACTCGGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAGTTTCGGGGCTGTTGTTCGTTCCAGCGGCTTAg                                                                                                                                                                                <  2:48684/251‑1 (MQ=255)                                                                      aCGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGGGGGAAGGCGGGCATATCCATACCGCCACCGCCCAGGATCACGGCATCGGGATCGAAAAGATTAAAGCTGGGGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTATGGCCGAAAGGAGCGTTTTCCGCCAGGACGAAAAGATCGCGCAATGGGAAATTTCGGGGCTGTTGTTCgaccccgagtctt                                                                                                                                                                                  >  2:676895/1‑237 (MQ=255)                                                                             gATGCGGCAGTGGTCGGCGAAGGGACGTTGGGGGCACGGCAAAGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCTCGTATTC‑AA‑GCAGACTCTGGAAGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCAtt                                                                                                                                                                         <  2:29/251‑1 (MQ=255)                                                                                   gCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCagag                                                                                                                                                                   <  1:517175/251‑1 (MQ=255)                                                                                        ggACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAAtt                                                                                                                                                              >  1:648427/1‑251 (MQ=255)                                                                                         gACGGAGCAGGAATTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACAGCCCCGGATTAACGCATCGGGATCGAACAGATTACTGCTGGTGGCAATGGCCAGTGCCGCGTTTTCCGA‑G‑AGATGCTGGCCGGAAGAGGGGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGTGGATGTTGTGCATACTAGCGTGATAGCGCCCTTCCagagtgaatg                                                                                                                                                             >  2:303111/1‑245 (MQ=255)                                                                                          aCGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTgg                                                                                                                                                            >  1:335361/1‑251 (MQ=255)                                                                                                  ggTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAgg                                                                                                                                                    >  1:636353/1‑251 (MQ=255)                                                                                                         ttGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCAccc                                                                                                                                               >  2:99037/1‑249 (MQ=255)                                                                                                               cATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGTACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCAGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGc                                                                                                                                       >  1:125247/1‑251 (MQ=255)                                                                                                                       cGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                                                               >  1:452694/1‑251 (MQ=255)                                                                                                                         ttagTCTCGCGAGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGACCAGATTAATGCTGGTGGCACTGGCCCGTGCCGCGTTTTC‑AA‑GCCGACTCTGGACGACAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCCATTGGAAATTTCGGGGCTGTTGTTCGTACACGCGGCTTCGCCCCATTCCAGGGCTATTAGTTTCCAGGcccaccagcttccccccccccc                                                                                                                              >  2:14178/3‑229 (MQ=255)                                                                                                                            gTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCCGCGTCGTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGGTTTCCCCACgcgcggcg                                                                                                                          >  2:32603/1‑247 (MQ=255)                                                                                                                                   gTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGCTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCAAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTc                                                                                                                  >  1:248994/1‑251 (MQ=255)                                                                                                                                      ggAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCa‑ta                                                                                                              <  1:443644/251‑1 (MQ=255)                                                                                                                                      ggAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCGGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGGTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCa‑ta                                                                                                              <  2:261153/251‑1 (MQ=255)                                                                                                                                        aaGGCGGGCATTTCCATCACGCCTCCGCCCAGGATCACCGCATCGTGAACGAACAGATGAATGCTGGTGGCAATTGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATAGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACTCGTGCAGTGTTGG‑GTCagat                                                                                                              <  2:384067/250‑4 (MQ=255)                                                                                                                                                 atatCCAGCACGCCACCGCCCAGGAGCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTATC‑AA‑GCAGACGCTGTACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCGAGGATTGCCACACGCGCATTGTTGG‑GTCA‑TATCTCC‑CAgggg                                                                                                  <  2:66808/251‑1 (MQ=255)                                                                                                                                                      cATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGatat                                                                                             >  1:593537/1‑251 (MQ=255)                                                                                                                                                          aCGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGAcc                                                                                         <  1:666551/251‑1 (MQ=255)                                                                                                                                                                   cccAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCa                                                                                                                                                                      <  1:291728/168‑1 (MQ=255)                                                                                                                                                                   cccAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCAGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGCGGCGTTTTCCGCATGGACGAAAAGACCGCGCCATGGGTACTTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTcaa                                                                                                                                                                      >  2:291728/1‑166 (MQ=255)                                                                                                                                                                    ccAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCt                                                                               <  1:622148/251‑1 (MQ=255)                                                                                                                                                                           cACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAGGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                        <  1:205863/251‑1 (MQ=255)                                                                                                                                                                           cACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                        >  1:598491/1‑251 (MQ=255)                                                                                                                                                                                 aTCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑Gt                                                                                                                   >  1:173488/1‑204 (MQ=255)                                                                                                                                                                                 aTCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑Gt                                                                                                                   <  2:173488/204‑1 (MQ=255)                                                                                                                                                                                      gATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                             >  1:516773/1‑251 (MQ=255)                                                                                                                                                                                       aTCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                                                            >  1:411342/1‑251 (MQ=255)                                                                                                                                                                                         agAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACCCGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCAAGTTCGCCTGCCACACCGTGTGCACCCGTc                                                          <  1:95218/250‑1 (MQ=255)                                                                                                                                                                                            aCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTCCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTTTTAGCGCCATTCCAGAGCAATTGTTTTCCAGGCCCCCAGGATTGCCACACGCGCAGTGTTGG‑GTAA‑TTTCCCC‑CAGGGGAATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAc                                                       <  1:324868/251‑1 (MQ=255)                                                                                                                                                                                                aTTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAg                                                                                                                                                     >  1:113144/1‑156 (MQ=255)                                                                                                                                                                                                aTTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAg                                                                                                                                                     <  2:113144/156‑1 (MQ=255)                                                                                                                                                                                                 ttAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGAAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGACCCCAGGATTGCCACACGCGCAGGGTTGG‑GTCA‑TATCTCCACAGGTG‑ATATGACCCAGTTCGCCTGCCAACCCGTGTGccaccttccaccgagc                                                  >  2:504931/1‑236 (MQ=255)                                                                                                                                                                                                        tggtggCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGGAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACACAGGATTGCCCAACGCCCAGGGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCaactttc                                           >  2:155737/1‑245 (MQ=255)                                                                                                                                                                                                         gggggCAATGGCCCGGGCCGCGTTTAC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACAAAAAGATCGCGCAAGGGGTAATTCCGGGGCTGTTGTTCTTCCCAGCTTCTTAGCGCCTTTCCAGAGCAATTGGTTCCCAGGCCCCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCCCC‑CAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                           <  1:379490/247‑1 (MQ=255)                                                                                                                                                                                                            ggCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTTTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGTTTTCCAGGCCCCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CGGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                       <  1:99037/251‑1 (MQ=255)                                                                                                                                                                                                             gCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGAAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑ATCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTAATcc                                       >  1:374372/1‑250 (MQ=255)                                                                                                                                                                                                              cAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTAATCcac                                     >  1:116999/1‑251 (MQ=255)                                                                                                                                                                                                                       gTGCCGC‑TTTTC‑AAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTGCGGGGTTGTTTTTCGTCCCAGCGCCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                            <  1:505445/251‑1 (MQ=255)                                                                                                                                                                                                                         gCCGCGTTTTC‑AA‑GCAGTCTTTGGACGAAAGAGGCTTTTTCCGCTGGGACAAAAAGATCGCGCAATGGGTAATTCCGGGGCTGTTGTTCGTACAAGGTTCTTAGCGCCATCCCAGAGAATTTGTTTTCCAGGCCCCCGGGATTGCCACACGCGCAGTGTTGG‑TTCA‑TATCTCC‑CAGGGGGATATGCCCCAGTTCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                          <  1:604633/251‑1 (MQ=255)                                                                                                                                                                                                                             cgTTTTC‑AA‑GCAGACTCTGGACAAAAGAGGCTTTTCCCGCATGGACGAAAAGATCGCGAAAGGGGTAATTTCGGGGCTGTTGTCCGTCCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTTTTGG‑TTCA‑TATCTCC‑CAGGGGGATATGACCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                      <  1:32603/251‑1 (MQ=255)                                                                                                                                                                                                                               ttttC‑AT‑GCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCCATTGGTTTCAAGGCACCCAGGATTGCCACACGCGCAGTGATGG‑ACCA‑TATCTAC‑CAGGGGGATATGACCCAGTTCGCCGGACACAACGTGTGCACCCGACCACGGGGCACCGTTCATCCAACCCGCGAACCCCATccc                    >  2:92318/1‑251 (MQ=255)                                                                                                                                                                                                                                                gggACGAAAGAGGTGTTTCCCGAATGGACGAAAAGACCGCGCAATGTGTATTTTCGGGGTTGTTGTTCGTCCCAGCGTCTTACCCCCATCCAAGAGCATGTGGTTGCCAGGCCCCCAGGATGCCCACCCGCGCAG‑GTTGGTGTCA‑TATCTCC‑CAGGGGGATATACCCCGGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg      <  1:92318/249‑1 (MQ=255)                                                                                                                                                                                                                                                   aCGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCAATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑G‑CACTATCTCC‑CAGGGGGATATGTCCCAGTTCGCCAGCCACACCGTGTGCACCCGTCCACGGCGCCCCGTTCATCCACACTGCGCACCCCATCCCCGCACCGAAATCGgccgc  >  2:670091/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       CGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTC‑AA‑GCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGG‑GTCA‑TATCTCC‑CAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGC  >  NC_000913/4306904‑4307394

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