Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,572,542	G→A	P322S (CCG→TCG)	eutD ←	phosphate acetyltransferase EutD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,572,542	0	G	A	96.8%	197.4 / NA	64	P322S (CCG→TCG)	eutD	phosphate acetyltransferase EutD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/21);  minor base C (2/0);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCTGCAACATCACCACGAACCATGGCAGTACACAGGCCGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACTTCCGGCACGAAGGCGGCG  >  NC_000913/2572299‑2572783                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   gCTGCAACATCACCACGAACCATGGCCGTACACAGCCCGCCGCCAATCTGCTTAACACAAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑CGCt                                                                                                                                                                                                                                            >  2:257055/1‑251 (MQ=255)      aaCATCACCACGAACCATGGCAGTACACAGGCCGCCGCCAAGCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGccac                                                                                                                                                                                                                                       <  2:389217/251‑1 (MQ=255)        cATCACCACGAACCATGGCAGTACAAAGGCCGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGccacc                                                                                                                                                                                                                                      >  2:201718/1‑250 (MQ=255)             ccacGAACCATGGCAGTACACAGACCGCCCCCAATATGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGCTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGCACGGGACCGTTCATTCAACCAGTGTTTGTAACCTGCTTTCGCGGTTCACTTCTGTCTGACGCGACCC‑AGATGCCACCAGCGcc                                                                                                                                                                                                                                >  2:413472/1‑251 (MQ=255)                  aaCCATGGCAGTACACAGGCCGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTc                                                                                                                                                                                                                           >  1:280546/1‑251 (MQ=255)                                  aGGCCGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACg                                                                                                                                                                                                           >  2:101270/1‑251 (MQ=255)                                    gccgccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCt                                                                                                                                                                                                         <  2:112304/251‑1 (MQ=255)                                      cgccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTAc                                                                                                                                                                                                       >  1:72255/1‑251 (MQ=255)                                             aaTCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACg                                                                                                                                                                                                >  1:273238/1‑251 (MQ=255)                                                   cTTAACACCACCCAACTTCACCCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCACTCAGTGCAACCAGGCCCCGGGTTTAGCTCCTTCCTCATGCTTCCATTCTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCCGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTCTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGCTAATTTCCTGCACGCTACCACCACgagag                                                                                                                                                                                            >  2:371935/1‑249 (MQ=255)                                                          ccaaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCCGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGc                                                                                                                                                                                   >  2:194980/1‑251 (MQ=255)                                                                  ctcgACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCACCCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGGTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCATCGACAGCTCGATAATTTCCTGCACCATACAACAACGACAGAGATAGTGCATAGgctc                                                                                                                                                                           >  2:317732/4‑249 (MQ=255)                                                                     gACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGc                                                                                                                                                                        >  1:83222/1‑251 (MQ=255)                                                                                       cATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTCCAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTgg                                                                                                                                                      >  2:202670/1‑251 (MQ=255)                                                                                          cgcgTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTAGCAGTGGcc                                                                                                                                                    <  1:362699/250‑1 (MQ=255)                                                                                          cgcgTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTGTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGcc                                                                                                                                                    <  2:159027/250‑1 (MQ=255)                                                                                              tCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGAACATTCCCAATGCTTCCACTGTGCTTTCATATTTATCAGGGTCCAGAATGAGACCGTTCATTCAACCAGTGTTTGTAAACTGCATTCGCGGCTCACTTCTGTCTGACGCGCCAC‑AGCTGCCACCAGCGCCAGCTCGATACTTTCCTGCACGATACAACCACGAGAGCGATCGCGCATCGGCGCGGCAAGTCCATGTATCAGTGGCCCgccg                                                                                                                                               >  1:250055/1‑248 (MQ=255)                                                                                               cAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCCCCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACgg                                                                                                                                              >  2:181342/1‑251 (MQ=255)                                                                                               cAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACg                                                                                                                                               <  2:4988/250‑1 (MQ=255)                                                                                                aGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGCGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCCCCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGCCGGc                                                                                                                                             >  1:42413/1‑251 (MQ=255)                                                                                                       tCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGatat                                                                                                                                      >  2:120136/1‑251 (MQ=255)                                                                                                          aTCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATccg                                                                                                                                   >  1:42341/1‑251 (MQ=255)                                                                                                            cAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTATTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATccgcc                                                                                                                                 >  2:95938/1‑251 (MQ=255)                                                                                                             aGTGCAACCAGGCCCCGGGTTGCGATCATTCCTAATGCTTACATTGTGCTTTCCTATGTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCGTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCGGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCCCGATATccgccg                                                                                                                                <  2:261898/251‑1 (MQ=255)                                                                                                                 cAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGACCGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTc                                                                                                                            <  2:42341/251‑1 (MQ=255)                                                                                                                           ccGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGAt                                                                                                                  <  2:83222/251‑1 (MQ=255)                                                                                                                            cGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTATTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGAtt                                                                                                                 <  1:95938/251‑1 (MQ=255)                                                                                                                              ggTTTCGATACTCCCTAATGCTTCCATGGGGCTTTCCCTTTTATCATGTTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACGGCGTTCGCGGTTCACTGCTGTCTGACGCGACAC‑AGCTGCCGCCAGCGCAAGCTCGATAATTTCCTTCACGCTCCAACCACGAGATAGATCGTGCATCGGCGCGGCAATTCTTTGTATCAGTGGCCCGGCGGCACGATATCCGCCGAGTCGTTGTGCGAtttt                                                                                                               <  2:131903/251‑1 (MQ=255)                                                                                                                               gTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATAAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGAttttt                                                                                                              >  1:407496/1‑250 (MQ=255)                                                                                                                                            cTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGccacca                                                                                                                                                                                                                                     <  1:112770/119‑1 (MQ=255)                                                                                                                                            cTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGccacca                                                                                                                                                                                                                                     >  2:112770/1‑119 (MQ=255)                                                                                                                                                    tCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATTATTTCCTGCACGCTACAACCACGAGAGAGATCGTGc                                                                                                                                                                                   >  2:288155/1‑161 (MQ=255)                                                                                                                                                    tCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGc                                                                                                                                                                                   <  1:288155/161‑1 (MQ=255)                                                                                                                                                     ccATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTcc                                                                                         <  1:183468/250‑1 (MQ=255)                                                                                                                                                            gCTTTCCTCTTTATCAGGGTCCAGAACTGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACgg                                                                                 >  1:332758/1‑251 (MQ=255)                                                                                                                                                            gCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGCCgg                                                                                 >  2:12869/1‑251 (MQ=255)                                                                                                                                                             cTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACgg                                                                                 >  1:353482/1‑250 (MQ=255)                                                                                                                                                              tttCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGa                                                                                >  2:390146/1‑250 (MQ=255)                                                                                                                                                                  ctctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAAc                                                                           >  1:94165/1‑251 (MQ=255)                                                                                                                                                                  ctctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAAc                                                                           >  1:78070/1‑251 (MQ=255)                                                                                                                                                                   tctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAAcc                                                                          >  1:97944/1‑251 (MQ=255)                                                                                                                                                                    ctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCCACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAAcc                                                                          >  1:228789/1‑250 (MQ=255)                                                                                                                                                                      ttATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATc                                                                       >  1:112480/1‑251 (MQ=255)                                                                                                                                                                                   gAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGc                                                          >  1:333228/1‑251 (MQ=255)                                                                                                                                                                                     aCGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑ATCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAACAACCCTCACATTGGCCTTGccc                                                        >  1:187500/1‑251 (MQ=255)                                                                                                                                                                                      cGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCt                                                       <  2:192225/251‑1 (MQ=255)                                                                                                                                                                                      cGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCt                                                       <  2:72255/251‑1 (MQ=255)                                                                                                                                                                                        ggACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGt                                                     >  1:218711/1‑251 (MQ=255)                                                                                                                                                                                                cattcaACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGCCAC‑AGCTGCCACAAGCGCCAGGTCGATAATTTCCTGCACGCTACAACCACGAGAGCGATCGTGCATAGGCGCGGCAAGGCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGGGACGTTGTGCGATTTTGTAACCAATATTTGCGGCTTCCAGCGACGGAAAAACCATCGCCTTGACCTTGCCCTGTAACGGGCt                                             >  1:38620/1‑251 (MQ=255)                                                                                                                                                                                                    caACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGc                                          >  1:191131/1‑250 (MQ=255)                                                                                                                                                                                                        cagttttttgaaaCTGCTTTCTCGGTTCACTTCTGTCTGACGCCACCC‑AGCTGCCACCAGCGCCAGCTCCATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTTTATCAGGGGCCCGACGGCTCGATATCCGCCGAGGCGTTGTGCGATTTTGTAACCAATATTTCCGGCGTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGgcg                                     <  1:413397/241‑1 (MQ=255)                                                                                                                                                                                                          gtgtTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCt                                                                                                                                                                 <  2:345173/125‑1 (MQ=255)                                                                                                                                                                                                          gtgtTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCt                                                                                                                                                                 >  1:345173/1‑125 (MQ=255)                                                                                                                                                                                                            gtTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGCGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTTGCCTTGCCCTGTAGCGGGCTGGCAGgcgcgtt                                 >  2:34022/1‑248 (MQ=255)                                                                                                                                                                                                             tttGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACCC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCttt                                 <  1:127334/250‑1 (MQ=255)                                                                                                                                                                                                               tGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTg                              >  2:344687/1‑251 (MQ=255)                                                                                                                                                                                                                 tAAACTGCTTTCGCGGTTCACTTCGGCCTGACGCG‑CACCAGCTGCCACCAGAGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTgc                             <  2:323648/250‑1 (MQ=255)                                                                                                                                                                                                                   aaCTACTTTCGCGGTTCACTTCTGGCTGACGCGACCC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGcgcc                          <  2:95953/251‑1 (MQ=255)                                                                                                                                                                                                                      tGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCCCATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCgccgcc                       >  2:231097/1‑251 (MQ=255)                                                                                                                                                                                                                      tGCTTTCGCGGTTCACTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCgccgcc                       >  2:103420/1‑251 (MQ=255)                                                                                                                                                                                                                         tttCGCGGTTCACTTCTGTATGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGACCGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGGAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACt                    >  2:343574/1‑251 (MQ=255)                                                                                                                                                                                                                               ggTTCCCTTCTGTCTGACGCGACAC‑AGCTGCCCCCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCTTTGTATCAGTGGCCCGGCGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACTTCCGGc              <  1:375799/251‑1 (MQ=255)                                                                                                                                                                                                                                    aCTTCTGTCTGACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACTTCCGGCACGa          >  1:52731/1‑250 (MQ=255)                                                                                                                                                                                                                                            cggACGCGACAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACTTCCGGCACGAAGgcggcg  <  1:268545/248‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GCTGCAACATCACCACGAACCATGGCAGTACACAGGCCGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCAC‑AGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATCAGTGGCCCGACGGCACGATATCCGCCGAGTCGTTGTGCGATTTTGTAACCAATATTTCCGGCTTCCAGCGACGGAAAAACCATCACATTGGCCTTGCCCTGTAGCGGGCTGGCAGGCGCTTTTTGCGCCGCCACTTCCGGCACGAAGGCGGCG  >  NC_000913/2572299‑2572783

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