Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,557,426	T→C	L418L (TTA→TTG)	eutB ←	ethanolamine ammonia‑lyase, large subunit, heavy chain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,557,426	0	T	C	100.0%	239.3 / NA	61	L418L (TTA→TTG)	eutB	ethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/24);  total (37/24)

CAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCA  >  NC_000913/2557279‑2557573                                                                                                                                                    |                                                                                                                                                    cAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCaa                                                                                                                                                   <  1:1035614/150‑1 (MQ=255) cAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAAc                                                                                                                                                  >  1:701437/1‑151 (MQ=255)   aTCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTg                                                                                                                                                >  1:312045/1‑151 (MQ=255)     cTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGAc                                                                                                                                              >  2:865542/1‑151 (MQ=255)       gTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGc                                                                                                                                            >  1:779692/1‑151 (MQ=255)        tttttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGc                                                                                                                                            >  2:778266/1‑150 (MQ=255)        tttttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGc                                                                                                                                            >  1:15310/1‑150 (MQ=255)          tttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGcac                                                                                                                                          <  1:865590/150‑1 (MQ=255)          tttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGcac                                                                                                                                          >  1:811512/1‑150 (MQ=255)          tttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGcac                                                                                                                                          >  1:521263/1‑150 (MQ=255)           ttGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGcac                                                                                                                                          >  1:1024774/1‑149 (MQ=255)              aTCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTgg                                                                                                                                     >  1:755675/1‑151 (MQ=255)                ccATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCg                                                                                                                                   >  1:309011/1‑151 (MQ=255)                 cATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCgg                                                                                                                                  >  1:662792/1‑151 (MQ=255)                  aTGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGt                                                                                                                                 >  1:727610/1‑151 (MQ=255)                  aTGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGt                                                                                                                                 <  1:656461/151‑1 (MQ=255)                   tGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCgg                                                                                                                                  >  1:975601/1‑149 (MQ=255)                        tGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTg                                                                                                                           >  1:948162/1‑151 (MQ=255)                           tATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGa                                                                                                                         <  2:989984/150‑1 (MQ=255)                           tATCTCCGCGTCATCAGAAGAACAGAGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCATTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTg                                                                                                                           >  2:967210/1‑148 (MQ=255)                                  gcgTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGTCCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCggtg                                                                                                                 <  2:296935/151‑1 (MQ=255)                                    gTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCggtggt                                                                                                               <  2:143615/151‑1 (MQ=255)                                     tcatcaGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTc                                                                                                              >  1:187846/1‑151 (MQ=255)                                           gaagaaCAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGt                                                                                                          >  1:241299/1‑149 (MQ=255)                                            aagaaCAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCCTGGAATGCGGGGGGCTGGTAGt                                                                                                       >  2:285018/1‑151 (MQ=255)                                               aaCAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCCGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTg                                                                                                     >  1:914686/1‑150 (MQ=255)                                                  aGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCa                                                                                                 <  1:290556/151‑1 (MQ=255)                                                   gTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCat                                                                                                >  1:123493/1‑151 (MQ=255)                                                     gACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCatg                                                                                               <  2:193653/150‑1 (MQ=255)                                                     gACGGATCGCCCGCCCGTTTGGTCAGGCCACCGTTTGCCTTAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGGTGAGCAACTGACGCACAGTGGCGGTATCGAGGAATGCGGTGGTCTAGTAGTTGAGCatga                                                                                              <  2:638752/151‑1 (MQ=255)                                                          aTCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGtca                                                                                         <  1:469878/151‑1 (MQ=255)                                                          aTCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGtc                                                                                          >  2:764364/1‑150 (MQ=255)                                                             gcccgcccGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGtcatc                                                                                       >  2:821895/1‑150 (MQ=255)                                                                cgcccgTTTGGTCAAGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGTCTGGTAGTTGAGCATGAAGTCATCAccc                                                                                   >  2:639745/1‑151 (MQ=255)                                                                       ttGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGc                                                                             >  1:739813/1‑150 (MQ=255)                                                                             aGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATcccc                                                                       >  1:716806/1‑150 (MQ=255)                                                                              ggCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCGAACTCCGGTGACGGGCGCAGGTTGAGCTACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCa                                                                      <  2:997954/150‑1 (MQ=255)                                                                               gCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCatg                                                                    <  2:63088/151‑1 (MQ=255)                                                                                 gACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCatgat                                                                  <  1:705924/151‑1 (MQ=255)                                                                                  aCCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCatgatg                                                                 >  1:969424/1‑151 (MQ=255)                                                                                      tttGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGt                                                             >  1:162206/1‑151 (MQ=255)                                                                                         gCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGc                                                          >  1:564260/1‑151 (MQ=255)                                                                                         gCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGc                                                          <  2:934301/151‑1 (MQ=255)                                                                                         gCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGc                                                          <  1:170780/151‑1 (MQ=255)                                                                                          ccATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGc                                                          >  1:645220/1‑150 (MQ=255)                                                                                                 gCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCGACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGGAGTTGCAGCCTGCg                                                  >  2:897542/1‑151 (MQ=255)                                                                                                 gCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGc                                                   <  2:995/150‑1 (MQ=255)                                                                                                  cccATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCGACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCgg                                                 <  1:897589/151‑1 (MQ=255)                                                                                                       gCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCg                                            >  1:661477/1‑151 (MQ=255)                                                                                                        cTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGGAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCCGCCTGCGGTGGCg                                            >  2:969342/1‑150 (MQ=255)                                                                                                          ttCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGa                                           <  2:486190/149‑1 (MQ=255)                                                                                                              agccagcGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGa                                     <  2:389563/151‑1 (MQ=255)                                                                                                               gccagcGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGat                                    <  2:851030/151‑1 (MQ=255)                                                                                                                  agcGTTCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGatcat                                 >  2:113434/1‑151 (MQ=255)                                                                                                                       tCAAACTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGt                            >  1:219199/1‑151 (MQ=255)                                                                                                                         aaaCTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGttt                          >  1:1012995/1‑151 (MQ=255)                                                                                                                         aaaCTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGttt                          <  2:975550/151‑1 (MQ=255)                                                                                                                          aaCTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGtttt                         >  1:350230/1‑151 (MQ=255)                                                                                                                           aCTCCGGTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTc                        <  2:461769/151‑1 (MQ=255)                                                                                                                                 gTGACGGGCGCAGGTTGAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAg                  <  1:117098/151‑1 (MQ=255)                                                                                                                                                 gAGCAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtca  <  1:632157/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCA  >  NC_000913/2557279‑2557573

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