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	100.0%	165.7 / NA	43	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/21);  total (22/21)

TCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAAT  >  NC_000913/803430‑803899                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              aacatcatcaacggtggtgagcATTTATTATCTTTCTGCGCACTTCTCGGTGCGGAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTTGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                                                                                                                                                                                                             <  2:358847/230‑1 (MQ=255)   gAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                                                                                                                                                                                                           <  1:528513/251‑1 (MQ=255)     attattTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCCCGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                                                                                                                                                                                                         >  1:227194/1‑251 (MQ=255)                   atCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGCCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                                                                                                                                                                                                            >  1:552303/1‑250 (MQ=255)                     ccATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGatat                                                                                                                                                                                                          <  2:545477/250‑1 (MQ=255)                                  ttCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt                                                                                                                                                                                                <  1:231734/247‑1 (MQ=255)                                    cTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAtt                                                                                                                                                                                          >  1:251639/1‑251 (MQ=255)                                        gcACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc                                                                                                                                                                                      >  1:333847/1‑251 (MQ=255)                                          aCTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc                                                                                                                                                                                      >  2:482616/1‑249 (MQ=255)                                                aggTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGcg                                                                                                                                                                              <  2:544176/250‑2 (MQ=255)                                                 ggTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACATACCGCt                                                                                                                                                                              >  2:516833/1‑250 (MQ=255)                                                   tGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTaa                                                                                                                                                                            >  2:539000/1‑250 (MQ=255)                                                         gATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTcc                                                                                                                                                                     <  2:251639/251‑1 (MQ=255)                                                                  aGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGa                                                                                                                                                             >  1:128424/1‑250 (MQ=255)                                                                                  gTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAGTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgttg                                                                                                                                                                                               >  1:330179/1‑200 (MQ=255)                                                                                  gTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAGTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgttg                                                                                                                                                                                               <  2:330179/200‑1 (MQ=255)                                                                                  gTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCACGCTGCCGt                                                                                                                                            >  2:208933/1‑251 (MQ=255)                                                                                    tATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGAAGAAGAACATTTCACAGGAGAAATTAAAAAAGAAGAAGCCAAAAAAGGCAATATTGCCTGGTCTATTCTCTCTTCGCAGAAAACGTCCGGAAATATGGATAAAATTAAAATAAAGTTTGATTCATTAGACTATCAAGATATTACCTTTGCTGGTATTGTAGAGAACGCTAAAGCGTCCGGTATGGAACGTTTACCGCTgacga                                                                                                                                            >  1:110932/1‑245 (MQ=255)                                                                                                gCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAAc                                                                                                                              <  2:559028/251‑1 (MQ=255)                                                                                                    gtTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGcc                                                                                                                          <  2:137208/251‑1 (MQ=255)                                                                                                    gtTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGcc                                                                                                                          >  1:545983/1‑251 (MQ=255)                                                                                                              aGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactcact                                                                                                                <  2:430063/251‑1 (MQ=255)                                                                                                               gTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactcactc                                                                                                               >  1:153836/1‑251 (MQ=255)                                                                                                                 aataaCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCActct                                                                                                              <  2:55106/250‑1 (MQ=255)                                                                                                                         aataatTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTcg                                                                                                     <  1:539000/251‑1 (MQ=255)                                                                                                                            aatTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTcggcg                                                                                                  >  2:78207/1‑251 (MQ=255)                                                                                                                                       aaCATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc                                                                                                                                                                                      >  1:571198/1‑156 (MQ=255)                                                                                                                                       aaCATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc                                                                                                                                                                                      <  2:571198/156‑1 (MQ=255)                                                                                                                                       aaCATTTCACCGGCAAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCTCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAGTAAc                                                                                       <  2:471458/251‑1 (MQ=255)                                                                                                                                                           aaaaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTgg                                                                   >  1:393452/1‑251 (MQ=255)                                                                                                                                                            aaaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGt                                                                  <  1:482616/251‑1 (MQ=255)                                                                                                                                                            aaaaaGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGt                                                                  >  2:300013/1‑251 (MQ=255)                                                                                                                                                                              aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAAATCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTTGTATTGTACAGACCGCTAAAGCGTCCGGGATGGAACGTTTCCCGCTGCCGGAGGTGCTTACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCTGCGGCCCAGCGt                                                <  2:227194/251‑1 (MQ=255)                                                                                                                                                                                  cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATg                                            >  1:187589/1‑251 (MQ=255)                                                                                                                                                                                         gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTCCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTAt                                     <  2:552303/251‑1 (MQ=255)                                                                                                                                                                                                  aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgtg                              <  1:516833/249‑1 (MQ=255)                                                                                                                                                                                                            cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCatat                   <  1:78207/250‑1 (MQ=255)                                                                                                                                                                                                               ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTAtttttttggcccctcacatttcc               >  2:254765/1‑233 (MQ=255)                                                                                                                                                                                                                    aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGAtttttttggcttcccctatttgggtcct          >  1:221221/1‑228 (MQ=255)                                                                                                                                                                                                                       tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTGTCTCCCCATCTTGTGGtcctccc       >  2:496716/1‑250 (MQ=255)                                                                                                                                                                                                                         cGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccatc     <  1:256855/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGAGTTTTTGTGCCttcccccattgcgggttttcatca    >  1:104132/1‑227 (MQ=255)                                                                                                                                                                                                                            aaaTATGGTAAAAATTAAAATTAAGTTTGAATCATTAGCCTCTCTCGATGTTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATGAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAat  <  2:544806/251‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              TCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAAT  >  NC_000913/803430‑803899

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