Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	45,854	C→T	H16H (CAC→CAT)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	45,854	0	C	T	100.0%	215.3 / NA	63	H16H (CAC→CAT)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/35);  total (28/35)

AGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGT  >  NC_000913/45614‑46099                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCAttt                                                                                                                                                                                                                                             <  1:169257/251‑1 (MQ=255)     tGACGGAAGTGTGCGCTTCGATTACGCCGGATGCCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCCTCGCCGCATTTTgtt                                                                                                                                                                                                                                         >  2:6569/1‑249 (MQ=255)       aCGGCAGTGTGCGCTTCGATTACGCCGGATGGCGGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCtgt                                                                                                                                                                                                                                       <  2:433378/251‑1 (MQ=255)             gtgnGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGaa                                                                                                                                                                                                                                 <  2:605245/251‑1 (MQ=255)              tgtgCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGa                                                                                                                                                                                                                                  >  1:395463/1‑249 (MQ=255)                        gATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCg                                                                                                                                                                                                                      <  1:382707/251‑1 (MQ=255)                         aTTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTACAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGGCCGt                                                                                                                                                                                                                     >  2:234313/1‑251 (MQ=255)                             cGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTGGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCt                                                                                                                                                                                                                 <  1:402614/251‑1 (MQ=255)                               ccGGATTTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGTCTCCTTTGGTGTGGAGGTCCGTTACGGCTGATGTTGGTTTGAGACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAGTAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGGGGGGAAGCGGTGGTCCGTTTCTgg                                                                                                                                                                                                               <  2:116127/251‑1 (MQ=255)                                  gATGTCTGGCGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCTAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATg                                                                                                                                                                                                            <  1:119556/251‑1 (MQ=255)                                     gTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATGCCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGGCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGtt                                                                                                                                                                                                         <  1:6569/251‑1 (MQ=255)                                       ccggAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGAGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTAAGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGAAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGAAGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGGTTCTGGATGGNTAt                                                                                                                                                                                                       >  2:590491/3‑251 (MQ=255)                                          gAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATGCATCGCCGCATTTTGCTGGGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTa                                                                                                                                                                                                    <  2:261372/251‑1 (MQ=255)                                             tGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTg                                                                                                                                                                                                 >  1:389050/1‑251 (MQ=255)                                                ggggccTGTCGCAGTGTGGGGCTGGGGTCAGCGCTGGAACAGTGGGAATGCCCGTGCGGCACATTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCGAAAAACAGGAATAACCATGCAACCGGCCAAAAGCTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATGTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTa                                                                                                                                                                                              <  2:460563/247‑1 (MQ=255)                                                       gTCGCATTCTGGGTCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGTAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAATCTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATtg                                                                                                                                                                                       <  2:429824/251‑1 (MQ=255)                                                       gTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGAt                                                                                                                                                                                         <  1:507962/249‑1 (MQ=255)                                                       gTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATtg                                                                                                                                                                                       >  1:83150/1‑251 (MQ=255)                                                           cATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATtggcgt                                                                                                                                                                                   >  2:506797/1‑251 (MQ=255)                                                            aTTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACGTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATtggcgtg                                                                                                                                                                                  <  2:264408/251‑1 (MQ=255)                                                                 ggggctggggTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCt                                                                                                                                                                             <  1:31178/251‑1 (MQ=255)                                                                             ggCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGAc                                                                                                                                                                 <  1:543620/251‑1 (MQ=255)                                                                                                   ccGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCt                                                                                                                                           <  1:560964/251‑1 (MQ=255)                                                                                                         ggCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTgg                                                                                                                                     <  1:418637/251‑1 (MQ=255)                                                                                                                       ggAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTgg                                                                                                                        <  1:344824/250‑1 (MQ=255)                                                                                                                               gTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGa                                                                                                                >  2:67680/1‑250 (MQ=255)                                                                                                                                 tACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAAc                                                                                                              <  1:506797/250‑1 (MQ=255)                                                                                                                                   cGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGACACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCt                                                                                                           >  2:423893/1‑251 (MQ=255)                                                                                                                                          tGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCggg                                                                                                    >  2:520464/1‑251 (MQ=255)                                                                                                                                             tGGTTTGATACGTAACGCCGCACTGACTATCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACGTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCt                                                                                                  <  2:310243/250‑1 (MQ=255)                                                                                                                                                 ttGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTc                                                                                              <  2:98068/250‑1 (MQ=255)                                                                                                                                                  tGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCg                                                                                             <  1:621195/250‑1 (MQ=255)                                                                                                                                                   gATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGtt                                                                                           <  1:615397/251‑1 (MQ=255)                                                                                                                                                   gATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGtt                                                                                           <  1:581251/251‑1 (MQ=255)                                                                                                                                                      aCGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTNggg                                                                                        >  2:587290/1‑250 (MQ=255)                                                                                                                                                          aaCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACAt                                                                                    <  2:186109/251‑1 (MQ=255)                                                                                                                                                               cgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTg                                                                               >  1:552224/1‑251 (MQ=255)                                                                                                                                                                    tGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTgg                                                                          >  1:33693/1‑251 (MQ=255)                                                                                                                                                                          tcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGcac                                                                                      >  2:449498/1‑233 (MQ=255)                                                                                                                                                                          tcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGcac                                                                                      <  1:449498/233‑1 (MQ=255)                                                                                                                                                                          tcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTatat                                                                    >  2:316683/1‑251 (MQ=255)                                                                                                                                                                          tcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTatat                                                                    >  1:383115/1‑251 (MQ=255)                                                                                                                                                                           cATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTtat                                                                      <  1:246299/248‑1 (MQ=255)                                                                                                                                                                                  aaaaaCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGAt                                                             >  2:470084/1‑250 (MQ=255)                                                                                                                                                                                  aaaaaCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATa                                                            <  1:458516/251‑1 (MQ=255)                                                                                                                                                                                   aaanCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGTGAACGCTCGCCGGGCGGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATaa                                                           <  2:608090/251‑1 (MQ=255)                                                                                                                                                                                    aaaCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATaaa                                                          <  1:550868/251‑1 (MQ=255)                                                                                                                                                                                       cagtaaTAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATTATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGTGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTc                                                       <  1:42534/247‑1 (MQ=255)                                                                                                                                                                                          gAATAACCATGCACACGTCCAGAAAATTTGACGCTCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGCTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTTAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGa                                                    <  1:330618/251‑1 (MQ=255)                                                                                                                                                                                          gAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCGCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGa                                                    <  2:553066/251‑1 (MQ=255)                                                                                                                                                                                            aTAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACg                                                  >  1:66899/1‑251 (MQ=255)                                                                                                                                                                                               aCCATGCAACCGTCAAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGATCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGAGAGAGTCGGgcggcg                                               >  2:204951/1‑251 (MQ=255)                                                                                                                                                                                                ccATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGgcgc                                              >  1:71671/1‑251 (MQ=255)                                                                                                                                                                                                     cAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAAt                                         >  1:637162/1‑251 (MQ=255)                                                                                                                                                                                                         cGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGTCTTGCTGGTCGCGGGAACGCCCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAACGTCGGACGGCGCAAAATGTTc                                     >  2:509434/1‑251 (MQ=255)                                                                                                                                                                                                         cGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                                     >  1:172991/1‑251 (MQ=255)                                                                                                                                                                                                         cGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCAAATCGCTGTTTGGTTATATTTCCGATAAAGGCGGACGGCGCAAAATNTTc                                     >  2:594165/1‑251 (MQ=255)                                                                                                                                                                                                              aGAAACTTTGACGATATCAAATTCTCCACTATTCATCGCCGCATTTTGATGTGGGGAAGCGGTGGACAGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGATGAAAATGGACGCTGACTGGATTGGCTTGCAGGGCGCGGGAAGGCCAGCCGGGCTGGTCGTTGGCACATCGCTGTTGGGTTATATTTACGATAAAGCCGGACGGCGCAAAATGttaatta                                 >  1:253236/1‑245 (MQ=255)                                                                                                                                                                                                                aaaCTTTGACGATCTCAAATTCTCCTCTATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTg                              >  2:565383/1‑251 (MQ=255)                                                                                                                                                                                                                                           tATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATc                                                                                   >  2:101570/1‑171 (MQ=255)                                                                                                                                                                                                                                           tATTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATc                                                                                   <  1:101570/171‑1 (MQ=255)                                                                                                                                                                                                                                            aTTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATc     <  1:423893/248‑1 (MQ=255)                                                                                                                                                                                                                                            aTTCATCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGt  >  2:4274/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      AGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGGCTGGGGTCGGCGCTGGAACAGTGGGAATACCCGCGCGGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGT  >  NC_000913/45614‑46099

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