Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	803,662	C→A	L54I (CTT→ATT)	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,662	0	C	A	98.4%	245.0 / NA	64	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/30);  total (34/30)
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.98e-01

GCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑ACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGT  >  NC_000913/803425‑803907                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       gCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACGATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGt                                                                                                                                                                                                                                            <  1:216716/250‑1 (MQ=255)       ggAATTATTTATAGTATATCCATTTATTATCTTTCTGCGCGCTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTAGCTGAAAAAGGCTTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCGCCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCGCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATAGATAA‑AATGAAAATTAAGTTTGAt                                                                                                                                                                                                                                       <  2:458575/248‑1 (MQ=255)       cGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTc                                                                                                                                                                                                                                     >  1:246835/1‑251 (MQ=255)       cGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACAGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGAtt                                                                                                                                                                                                                                      <  2:90900/250‑1 (MQ=255)             attTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGc                                                                                                                                                                                                                               <  2:286591/251‑1 (MQ=255)             attTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATATGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTatc                                                                                                                                                                                                                               >  2:328950/1‑249 (MQ=255)              tttAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGcc                                                                                                                                                                                                                              >  2:41230/1‑251 (MQ=255)                 agagTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCctc                                                                                                                                                                                                                            <  2:30578/250‑1 (MQ=255)                    gTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAggtctttcgctggggatttcct                                                                                                                                                                                                                        >  1:143473/1‑230 (MQ=255)                        atCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                                                                                                                                                                                                                     <  2:233812/250‑1 (MQ=255)                        atCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                                                                                                                                                                                                                     >  1:119714/1‑250 (MQ=255)                          ccATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGatan                                                                                                                                                                                                                   <  2:246835/250‑2 (MQ=255)                                       ttCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                                                                                                                                                                                                      >  1:430994/1‑250 (MQ=255)                                             gcACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc                                                                                                                                                                                               >  2:436623/1‑251 (MQ=255)                                                          cgcAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGc                                                                                                                                                                                  >  2:508951/1‑251 (MQ=255)                                                          cgcAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGc                                                                                                                                                                                  >  2:355382/1‑251 (MQ=255)                                                                  tCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTCAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGt                                                                                                                                                                           >  1:521511/1‑250 (MQ=255)                                                                     ggAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGt                                                                                                                                                                           <  1:508951/247‑1 (MQ=255)                                                                     ggAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTAt                                                                                                                                                                         <  1:477901/249‑1 (MQ=255)                                                                               atgatCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTccc                                                                                                                                                             <  2:140071/251‑1 (MQ=255)                                                                                    cAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGc                                                                                                                                                        >  2:375975/1‑251 (MQ=255)                                                                                     aaGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGcc                                                                                                                                                       <  1:387731/251‑1 (MQ=255)                                                                                              gAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTAtgtg                                                                                                                                                <  1:424775/249‑1 (MQ=255)                                                                                                aaaaGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTg                                                                                                                                             >  1:182880/1‑250 (MQ=255)                                                                                                         gtTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGcc                                                                                                                                   >  1:446504/1‑251 (MQ=255)                                                                                                          tttCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGcc                                                                                                                                   >  1:540069/1‑250 (MQ=255)                                                                                                                 ccAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactca                                                                                                                           <  2:182880/251‑1 (MQ=255)                                                                                                                   aGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAGAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactcact                                                                                                                         >  1:279322/1‑251 (MQ=255)                                                                                                                      aataaCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTg                                                                                                                      >  1:50259/1‑251 (MQ=255)                                                                                                                        aaaCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTgcg                                                                                                                    <  2:50259/250‑1 (MQ=255)                                                                                                                          aCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGcc                                                                                                                  <  2:119714/251‑1 (MQ=255)                                                                                                                            gAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGt                                                                                                                <  1:554862/251‑1 (MQ=255)                                                                                                                                         aagaaCATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAtt                                                                                                   >  1:506678/1‑251 (MQ=255)                                                                                                                                             aCATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACg                                                                                               >  2:46975/1‑251 (MQ=255)                                                                                                                                                  tCACCGGCGCGATTAAAAGAGAAGAAGCCAACAGAGGCACTATTGCCTGGGCTATTCTCTCTTCGCATAATAAGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGCGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGTGTCACTATTAACGGtgtt                                                                                          >  2:26657/1‑249 (MQ=255)                                                                                                                                                          gAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGt                                                                                  >  1:259558/1‑251 (MQ=255)                                                                                                                                                          gAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGATAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGACACTATTAACGGTGATGACCATGt                                                                                  >  1:376899/1‑251 (MQ=255)                                                                                                                                                            aaTTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTc                                                                                                               <  2:530900/220‑1 (MQ=255)                                                                                                                                                            aaTTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTc                                                                                                               >  1:530900/1‑220 (MQ=255)                                                                                                                                                                aaaaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTgg                                                                            >  2:325548/1‑251 (MQ=255)                                                                                                                                                                aaaaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTg                                                                             <  1:434704/250‑1 (MQ=255)                                                                                                                                                                  aaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTa                                                                                                                                                    >  2:382444/1‑177 (MQ=255)                                                                                                                                                                  aaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGtt                                                                          <  2:506678/251‑1 (MQ=255)                                                                                                                                                                  aaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTCCCCGCTGCCGTa                                                                                                                                                    <  1:382444/177‑1 (MQ=255)                                                                                                                                                                   aaaGAAGAAGCCAAGAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGttt                                                                         <  2:279322/251‑1 (MQ=255)                                                                                                                                                                    aagaagaagCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTa                                                                        <  1:233307/251‑1 (MQ=255)                                                                                                                                                                         gaagCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATcggc                                                                   <  1:263762/251‑1 (MQ=255)                                                                                                                                                                          aagCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATcggcg                                                                  <  1:436623/251‑1 (MQ=255)                                                                                                                                                                           agCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATcggcgg                                                                 >  1:362476/1‑251 (MQ=255)                                                                                                                                                                                        aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATgg                                                    >  1:18799/1‑251 (MQ=255)                                                                                                                                                                                           aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCgg                                                 <  2:54735/251‑1 (MQ=255)                                                                                                                                                                                             tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑t                                               <  1:72329/250‑1 (MQ=255)                                                                                                                                                                                             tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑t                                               <  2:111023/250‑1 (MQ=255)                                                                                                                                                                                                  ggTCTATTCTCTCTTCGCATAATACATCCGGAAATATGGATAACACT‑AAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACATTTGTTGGTATCGTACAGACCGCTAAAGCGTCCGGTATGGAACGCTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGCCGCCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGTCGGCCCCGCGttcttgcggccttctt                                          >  2:116819/1‑237 (MQ=255)                                                                                                                                                                                                      tATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GAtttttttgg                                     >  1:92109/1‑247 (MQ=255)                                                                                                                                                                                                        ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑AATTACAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCt                                   >  1:150083/1‑251 (MQ=255)                                                                                                                                                                                                                tCGCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCa                              <  1:41230/248‑1 (MQ=255)                                                                                                                                                                                                                  gCATAATACGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATAttn                         <  2:259558/251‑2 (MQ=255)                                                                                                                                                                                                                    ataataCGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGGCTCCGCATAtttcg                       >  1:116925/1‑248 (MQ=255)                                                                                                                                                                                                                      aataCGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGGCTCCGCCTATTGCggg                     >  1:108469/1‑250 (MQ=255)                                                                                                                                                                                                                         aCGTCCGGAAATATGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTATTTTTTTGCTCCCCCATATTgggttct                   >  2:57822/1‑245 (MQ=255)                                                                                                                                                                                                                                    tatGGATAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGGCTCCGCagttttccgttcttcctacctat        >  2:83040/1‑228 (MQ=255)                                                                                                                                                                                                                                         aTAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCGCATTTTTCGGGCATCCCTCAATATTTGCgg  >  1:15074/1‑250 (MQ=255)                                                                                                                                                                                                                                         aTAA‑AATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGt  <  2:92109/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       GCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAA‑ACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGT  >  NC_000913/803425‑803907

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