Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,426,929	C→T	C165Y (TGC→TAC)	yjfZ ←	DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,426,929	0	C	T	100.0%	164.0 / NA	48	C165Y (TGC→TAC)	yjfZ	DUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/24);  total (24/24)

GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAA  >  NC_000913/4426689‑4427166                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              ggCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCCTTCGCTCGAAACCCAAACGCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTGCAAAACATGGGGTAAATCAGAATACGCGTTTAGCAATTCGGCAATGATGAAACCGAGGTCATTATAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGACCGAGTGTTGGTAGTTCAGATGACCAGTAATTGGGCAt                                                                                                                                                                                                                                     <  2:703839/251‑1 (MQ=255)     ttAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGCGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATACAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTgg                                                                                                                                                                                                                                 >  1:497367/1‑251 (MQ=255)                cGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGt                                                                                                                                                                                                                      <  2:56876/251‑1 (MQ=255)                  aGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTGGGGCAGGAAGGGTTa                                                                                                                                                                                                                    <  2:361650/251‑1 (MQ=255)                        aCTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGt                                                                                                                                                                                                              >  2:211200/1‑251 (MQ=255)                              gTCTCGCTCCTGCCGCACGCGCTCGAAACCCAAACCCGAGTAGATGGCGGCTGAACTTGACAAGGAACGCAGACCAATTTTCAAATCATGGGGGAAATCAGAATACGCGTTTAGCAAAGCAGCAATGATGAAGCCGAGGGCATTAGAAACGCCGGTAACACGGGGAAGTTGTTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAGTTAGGCATTGGGGTATGAAgggctagaacttttccga                                                                                                                                                                                                        <  2:658376/251‑16 (MQ=255)                               tctcGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAACTCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGCTGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTAt                                                                                                                                                                                                       >  1:47639/1‑251 (MQ=255)                               tctcGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAACGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTAt                                                                                                                                                                                                       >  1:632482/1‑251 (MQ=255)                                         gCGGCACTCGCGCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTggtt                                                                                                                                                                                              <  2:93515/250‑3 (MQ=255)                                           ggCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACACCGTATCCTATCGTTTGGCCGGt                                                                                                                                                                                           >  1:368705/1‑251 (MQ=255)                                             cACTCGCTCGACACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCTTTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCggaat                                                                                                                                                                                         >  2:250134/1‑248 (MQ=255)                                                        aacccaaacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTa                                                                                                                                                                              <  1:601035/251‑1 (MQ=255)                                                           ccaaacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTAcc                                                                                                                                                                            >  1:494960/1‑250 (MQ=255)                                                              aacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATaa                                                                                                                                                                        <  2:128407/251‑1 (MQ=255)                                                                  cGAGTAGCTTGAGGCTGCAATTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGAAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACAATAAcc                                                                                                                                                                      >  2:419985/1‑249 (MQ=255)                                                                  cGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTGGGCCGGTATAGTCTGGAGTACCATAACCAg                                                                                                                                                                    <  2:466901/251‑1 (MQ=255)                                                                   gAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAgg                                                                                                                                                                   <  2:251196/251‑1 (MQ=255)                                                                     gTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTc                                                                                                                                                                 <  1:369292/251‑1 (MQ=255)                                                                        gATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTccc                                                                                                                                                               >  1:696906/1‑250 (MQ=255)                                                                                               aaCGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtggatggtgtat                                                                                                                                       <  2:696906/251‑1 (MQ=255)                                                                                                            tttCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGtt                                                                                                                          >  1:620511/1‑251 (MQ=255)                                                                                                             tacAAAGCATGGGGTAAATCAGAATACACGTTTAGCAATGCAGCAATGATGAATACGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTAGCTTTTGGCCTGTATAGTCTGGAGTACCATAAGCAGGTCCCAAAACCACTGCGTGTATTGTGTATTGACTGTGGGGtt                                                                                                                          <  2:326207/248‑1 (MQ=255)                                                                                                                       ggggTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTGCTGTTGGGTTCGAAGGATTCg                                                                                                               <  2:129306/251‑1 (MQ=255)                                                                                                                              aTCAGAATACGCGTTTAGCAAGGCAGCAATGAGGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCtt                                                                                                        <  2:599379/251‑1 (MQ=255)                                                                                                                                  gAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaa                                                                                                    <  2:632388/251‑1 (MQ=255)                                                                                                                                   aaTACGCGTCTAGCAATGCCGCAATGATGAATCCGTGGGCATTAGAAACGCCTGTCACAGTGCCAAGTTGGAAGATCTAAAGATCGATTGTTGGTAGTTCAGATGTCCAGTAATTAGGCATTTGGGCAGGAAGGGTTACACCGTCTCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCCCTTCGTGTATTGTGTATTTACTGTTGGGTTAGAAGGATTCGACATCTTGCaat                                                                                                   >  2:716172/1‑251 (MQ=255)                                                                                                                                     tACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTg                                                                                                                              <  1:116036/222‑1 (MQ=255)                                                                                                                                           tttAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATAttt                                                                                           >  1:487406/1‑251 (MQ=255)                                                                                                                                            ttAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTTCAACGTATCATATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGCAGGATTCGACATCTTGCAATAACTATTTc                                                                                          >  2:583226/1‑251 (MQ=255)                                                                                                                                                   tGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTAc                                                                                                                                                                                                                   >  2:513908/1‑123 (MQ=255)                                                                                                                                                   tGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTAc                                                                                                                                                                                                                   <  1:513908/123‑1 (MQ=255)                                                                                                                                                      agcaATGATGAATCCGCGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCgg                                                                                >  1:148970/1‑251 (MQ=255)                                                                                                                                                                   ccGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCAc                                                                    <  1:68094/250‑1 (MQ=255)                                                                                                                                                                      aGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCTGTTAACCTCGCACCcgc                                                                <  2:494960/251‑1 (MQ=255)                                                                                                                                                                          cATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTa                                                             >  2:18546/1‑250 (MQ=255)                                                                                                                                                                              aGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTGCTGCGGTGAACCTCGCACCCGCGTAAGCtt                                                        <  1:599706/251‑1 (MQ=255)                                                                                                                                                                               gaatcccctgaaaaagtgctaggTTGGTCTTTATCAAGATCGACTGTTGGTAGCTCAGATGAATAGTAATTAGGCATTTGGGCAGGAAGGGATACAACGTATCCTATCGCTTGTCCGACATAGTCTCGAATACCATAACCAGTTCCCAAAACCACTTCGTGAATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTCCGTTTACCCTCGCACCAGCGTAAGCttt                                                       >  2:103377/23‑251 (MQ=255)                                                                                                                                                                                                 cGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCGCCCGCGTAAGCTTTATTGTATTTCTGGGCGTg                                     <  2:233719/251‑1 (MQ=255)                                                                                                                                                                                                   aaGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtatt                                                                                                                                             <  2:65865/145‑1 (MQ=255)                                                                                                                                                                                                   aaGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtatt                                                                                                                                             >  1:65865/1‑145 (MQ=255)                                                                                                                                                                                                   aaGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGa                                    >  1:400951/1‑250 (MQ=255)                                                                                                                                                                                                             atatAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTa                          >  2:485120/1‑250 (MQ=255)                                                                                                                                                                                                                taAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGAATCCTATCGTTTGGCCGGTATAGTCTGCAGTACAATACCCAGGTCCCAAAACCACTTCGTGTCTTGTGTATTTACTGTTGGGTTCGCAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAAACTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATATAATGt                      >  1:692203/1‑251 (MQ=255)                                                                                                                                                                                                                  aaGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAg                    <  1:211200/251‑1 (MQ=255)                                                                                                                                                                                                                      tCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAACACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTc                >  2:407153/1‑251 (MQ=255)                                                                                                                                                                                                                                   gTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTtccatcca   <  2:497367/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTtccatcc    >  2:591276/1‑250 (MQ=255)                                                                                                                                                                                                                                     tCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCaa  >  2:97863/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAA  >  NC_000913/4426689‑4427166

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