Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,048,541	A→T	6.2%	T1202S (ACC→TCC)	yeeJ →	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,048,541	0	A	T	6.2%	227.5 / 4.8	64	T1202S (ACC→TCC)	yeeJ	putative adhesin
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base T (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

AGCGCATGTCACGCTGAAAGGTAAAAAAGCGGGCACGCATACGGTTACCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACTTCACCAGCAACGCAGCGACAG  >  NC_000913/2048258‑2048826                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              aGCGCATGTCACGCTGAAAGGTAAAAAAGCGGGCACGCATACGGTTACCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAAc                                                                                                                                                                                                                                                                                          <  1:359487/289‑1 (MQ=255)             gCTGAAAGGTAAAAAAGCGGGCACGCATACGGTTACCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCAt                                                                                                                                                                                                                                                                               >  2:658376/1‑288 (MQ=255)                  aaGGTAAAAAAGCAGGCACGCAAACGGTTACCGCAACGCTGGCTAATAACCATCCCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATcgcg                                                                                                                                                                                                                                                                           <  2:169683/287‑1 (MQ=255)                               gggCACGCATACGGTTACCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATTCCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCa                                                                                                                                                                                                                                                            >  1:542210/1‑289 (MQ=255)                                               aCCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATcgcg                                                                                                                                                                                                                                                                           >  2:359057/1‑258 (MQ=255)                                               aCCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATcgcg                                                                                                                                                                                                                                                                           <  1:359057/258‑1 (MQ=255)                                               aCCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTg                                                                                                                                                                                                                                            >  1:565633/1‑289 (MQ=255)                                               aCCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTg                                                                                                                                                                                                                                            >  1:623976/1‑289 (MQ=255)                                               aCCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTg                                                                                                                                                                                                                                            >  1:491793/1‑289 (MQ=255)                                                 cGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTgag                                                                                                                                                                                                                                          <  2:244905/289‑1 (MQ=255)                                                 cGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTgag                                                                                                                                                                                                                                          >  1:256765/1‑289 (MQ=255)                                                    aaCGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGagaa                                                                                                                                                                                                                                        <  1:691714/288‑1 (MQ=255)                                                           ggTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGCGTCTGGCATCGATCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGACGACGGTTa                                                                                                                                                                                                                                >  2:190143/1‑289 (MQ=255)                                                           ggTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTa                                                                                                                                                                                                                                >  2:671128/1‑289 (MQ=255)                                                                  aCAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATc                                                                                                                                                                                                                         >  1:443444/1‑289 (MQ=255)                                                                   cAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCa                                                                                                                                                                                                                        <  2:654270/289‑1 (MQ=255)                                                                   cAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCa                                                                                                                                                                                                                        <  2:565633/289‑1 (MQ=255)                                                                   cAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCa                                                                                                                                                                                                                        <  1:658354/289‑1 (MQ=255)                                                                     aTACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACTCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCCCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTACTACCAACGAGTCTGCCATCGCGCACGCCACTCGCGCAGGCGTTTCCTATGGCAAGAATAGCGTTACTGCTTccat                                                                                                                                                                                                                      >  2:116773/1‑286 (MQ=255)                                                                     aTACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATcgcg                                                                                                                                                                                                                                                                           <  2:470747/236‑1 (MQ=255)                                                                     aTACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATcgcg                                                                                                                                                                                                                                                                           >  1:470747/1‑236 (MQ=255)                                                                           gTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCtaa                                                                                                                                                                                                                >  1:270238/1‑289 (MQ=255)                                                                                        ggTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCATCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGc                                                                                                                                                                                                   <  1:89530/289‑1 (MQ=255)                                                                                               tttGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGacaaca                                                                                                                                                                                            <  1:658376/289‑1 (MQ=255)                                                                                                        gACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGc                                                                                                                                                                                   >  1:421052/1‑289 (MQ=255)                                                                                                                              tCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAgcggc                                                                                                                                                             >  1:48604/1‑289 (MQ=255)                                                                                                                                    aGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCaa                                                                                                                                                        <  2:329905/288‑1 (MQ=255)                                                                                                                                        atCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCa                                                                                                                                                                                                                                                                                >  2:332735/1‑164 (MQ=255)                                                                                                                                        atCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCa                                                                                                                                                                                                                                                                                <  1:332735/164‑1 (MQ=255)                                                                                                                                        atCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAAtt                                                                                                                                                   >  1:672934/1‑289 (MQ=255)                                                                                                                                        atCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAAtt                                                                                                                                                   >  1:765381/1‑289 (MQ=255)                                                                                                                                                  gagaTCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTg                                                                                                                                          >  1:152425/1‑288 (MQ=255)                                                                                                                                                     aTCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGgcgc                                                                                                                                      <  2:672934/289‑1 (MQ=255)                                                                                                                                                      tCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGgcgc                                                                                                                                      <  2:443444/288‑1 (MQ=255)                                                                                                                                                          aGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTg                                                                                                                                                                                                                     >  2:11651/1‑205 (MQ=255)                                                                                                                                                          aGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTg                                                                                                                                                                                                                     <  1:11651/205‑1 (MQ=255)                                                                                                                                                               aTGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAg                                                                                                                             >  2:230377/1‑288 (MQ=255)                                                                                                                                                                tGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGAc                                                                                                                           <  1:671128/289‑1 (MQ=255)                                                                                                                                                                tGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGAc                                                                                                                           >  1:647506/1‑289 (MQ=255)                                                                                                                                                                      cGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCg                                                                                                                                                                                                                                                        <  2:748963/158‑1 (MQ=255)                                                                                                                                                                      cGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCg                                                                                                                                                                                                                                                        >  1:748963/1‑158 (MQ=255)                                                                                                                                                                         tAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAAt                                                                                                                   <  2:152425/288‑1 (MQ=255)                                                                                                                                                                                aCGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTa                                                                                                           <  2:270238/289‑1 (MQ=255)                                                                                                                                                                                                  aaaGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGc                                                                                                                                                                                                   <  1:438216/183‑1 (MQ=255)                                                                                                                                                                                                  aaaGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGc                                                                                                                                                                                                   >  2:438216/1‑183 (MQ=255)                                                                                                                                                                                                         aGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATTCCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGgcgc                                                                                                                                      >  1:374898/1‑237 (MQ=255)                                                                                                                                                                                                         aGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATTCCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGgcgc                                                                                                                                      <  2:374898/237‑1 (MQ=255)                                                                                                                                                                                                          gTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGtc                                                                                 >  1:21941/1‑289 (MQ=255)                                                                                                                                                                                                          gTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGt                                                                                  >  2:365329/1‑288 (MQ=255)                                                                                                                                                                                                                cAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCAcc                                                                           >  1:500091/1‑289 (MQ=255)                                                                                                                                                                                                                  aTGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATTCCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGc                                                                         <  2:542210/289‑1 (MQ=255)                                                                                                                                                                                                                       gTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCcaca                                                                     <  2:500091/288‑1 (MQ=255)                                                                                                                                                                                                                             aagcTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTg                                                              <  2:256765/286‑1 (MQ=255)                                                                                                                                                                                                                              atCTTCCGTTAACATTCAGCTCGGCCTCTTAAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGa                                                             <  2:307523/289‑1 (MQ=255)                                                                                                                                                                                                                                    cGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAAc                                                                                                                                                                                                                                                                                          <  2:171427/62‑1 (MQ=255)                                                                                                                                                                                                                                    cGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAAc                                                                                                                                                                                                                                                                                          >  1:171427/1‑62 (MQ=255)                                                                                                                                                                                                                                         aaaTTCAGCTCAGCCTCTTCAGGACTCACCATTACCCCGGGAGTAAGTAATACCAACGAGGCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTGTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCt                                                  <  2:421052/287‑1 (MQ=255)                                                                                                                                                                                                                                                  tcagcCGCTTCAGGACTCACCCTGACCCCGGGGGTAAGTAATACCAACGAGTCTGGCATCGGGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTg                                          <  2:161406/288‑1 (MQ=255)                                                                                                                                                                                                                                                                  tCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGAGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACCGTCGTTGATACTACTGGCTTTCCGGTGAAAGGTGTGCCTGTGa                          >  1:559215/1‑288 (MQ=255)                                                                                                                                                                                                                                                                    aCCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGATAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACt                       >  1:722939/1‑289 (MQ=255)                                                                                                                                                                                                                                                                    aCCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACt                       >  2:136764/1‑289 (MQ=255)                                                                                                                                                                                                                                                                    aCCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACt                       >  2:167482/1‑289 (MQ=255)                                                                                                                                                                                                                                                                          aCCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACTTCACCa                 >  1:204642/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                          aTACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCTGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGGGAAATGTGTTACTGTGAACTTCCCCAGCACCGCcgcgccca  >  2:412758/1‑284 (MQ=255)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              AGCGCATGTCACGCTGAAAGGTAAAAAAGCGGGCACGCATACGGTTACCGCAACGCTGGGTAATAACAATACCAGTGATTCGCAGCCGGTAACATTTGTGGCGGACAAAGCCTCGGCTCAGGTTGTCCTGCAGATATCAAAAGATGAGATCACAGGTAATGGCGTCGATAGCGCAACGCTAACTGCAACGGTTAAAGATCAGTTCGACAATGAGGTGAATAATCTTCCGGTAACATTCAGCTCAGCCTCTTCAGGACTCACCCTGACCCCGGGAGTAAGTAATACCAACGAGTCTGGCATCGCGCAGGCCACTCTCGCAGGCGTTGCCTTTGGTGAGAAGACGGTTACTGCATCACTGGCTAATAATGGTGCCAGCGACAACAAAACTGTGCATTTTATTGGCGACACAGCGGCGGCAAAAATTATCGAGTTGGCGCCTGTCCCAGACAGCATAATCGCCGGTACCCCGCAGAACAGCTCCGGCAGCGTCATCACCGCCACAGTCGTTGATAATAATGGCTTTCCGGTGAAAGGTGTGACTGTGAACTTCACCAGCAACGCAGCGACAG  >  NC_000913/2048258‑2048826

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