Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,658,260	T→C	C64C (TGT→TGC)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,260	0	T	C	100.0%	216.4 / NA	59	C64C (TGT→TGC)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/39);  total (20/39)

CCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCG  >  NC_000913/1658124‑1658396                                                                                                                                         |                                                                                                                                         ccGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTCCGTCAAATGTGGGAGCCGCTgcgcc                                                                                                                                       >  2:1493357/1‑139 (MQ=255)  cgcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCAc                                                                                                                                      <  2:1480573/140‑1 (MQ=255)   gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACt                                                                                                                                     <  2:1263776/140‑1 (MQ=255)   gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACt                                                                                                                                     <  1:468683/140‑1 (MQ=255)       tAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGt                                                                                                                                 <  1:1160456/140‑1 (MQ=255)         ggTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCt                                                                                                                               <  2:655911/140‑1 (MQ=255)         ggTTCTCTGGCGCTGGCTGCAGGCGGCTTTTCTTTTCCGGTTACCCCGCGCAATTCAGCAGCAGCGGGACCACCCGCCCCCGAAAAAGGGGGCTGGGGGGCCCGTTCCGTCAACCTCGGTAGCCGCTGCGCGCTTCgcat                                                                                                                               >  2:869695/1‑137 (MQ=255)          gTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAACGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTa                                                                                                                              <  2:1213442/140‑1 (MQ=255)           ttCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTAc                                                                                                                             <  2:98615/140‑1 (MQ=255)              tctGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATg                                                                                                                          <  2:723701/140‑1 (MQ=255)                tGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGtt                                                                                                                        <  2:662456/140‑1 (MQ=255)                tGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGtt                                                                                                                        <  2:1252701/140‑1 (MQ=255)                  gcgcTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTaa                                                                                                                      <  2:1197509/140‑1 (MQ=255)                   cgcTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTaaa                                                                                                                     <  2:1018550/140‑1 (MQ=255)                        gctgCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATaa                                                                                                                <  2:1274328/140‑1 (MQ=255)                          tgCAGGCGGTTTTTTTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAATGGTCTGGGGGGCCCGGTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATg                                                                                                              >  2:701441/1‑140 (MQ=255)                            cAGGCGGTTTTTCTTTGCCGTTTACCCCGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGGAAAAGTGGTCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                            >  2:1472013/1‑140 (MQ=255)                                  gTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGAcc                                                                                                      >  2:1068976/1‑140 (MQ=255)                                         tttGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTgg                                                                                               <  1:135734/140‑1 (MQ=255)                                           tGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaa                                                                                             <  1:632769/140‑1 (MQ=255)                                           tGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaa                                                                                             <  1:393970/140‑1 (MQ=255)                                             ccGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAAc                                                                                           >  2:1089782/1‑140 (MQ=255)                                                    cccTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAAt                                                                                    <  2:426458/140‑1 (MQ=255)                                                       tGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATAcc                                                                                 >  1:1185075/1‑140 (MQ=255)                                                        gcgcAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCg                                                                                >  1:278721/1‑140 (MQ=255)                                                          gcAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGc                                                                              >  1:293681/1‑140 (MQ=255)                                                              tGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCg                                                                          >  2:105822/1‑140 (MQ=255)                                                                 cgccgccgcGGTACCACAGGCCCGCGAAAAAGTGGGCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATg                                                                       <  2:520180/133‑1 (MQ=255)                                                                  gcagcagcGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGa                                                                      <  2:278721/140‑1 (MQ=255)                                                                      cagcGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTAc                                                                  >  1:759338/1‑140 (MQ=255)                                                                          ggTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCa                                                              <  1:869695/140‑1 (MQ=255)                                                                              caacaGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCa                                                          <  1:1449404/140‑1 (MQ=255)                                                                                acaGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATc                                                        >  1:667920/1‑140 (MQ=255)                                                                                acaGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATc                                                        >  2:336414/1‑140 (MQ=255)                                                                                   ggCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAgg                                                     >  1:1508998/1‑140 (MQ=255)                                                                                    gCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGt                                                    <  1:1068976/140‑1 (MQ=255)                                                                                    gCCCGCGAAAAAGCGGTCTGGGGTGCCTGTTACGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGAGGAAACTGACAATACCGGCCGCGGTGAGAACGGCAACCATCAgga                                                    >  2:56701/1‑139 (MQ=255)                                                                                     cccGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTa                                                   <  2:960294/140‑1 (MQ=255)                                                                                      ccGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTAc                                                  >  1:599856/1‑140 (MQ=255)                                                                                      ccGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTAc                                                  >  1:43955/1‑140 (MQ=255)                                                                                         cgAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTAcgcg                                               >  2:957252/1‑140 (MQ=255)                                                                                              aaGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGt                                          >  2:66826/1‑140 (MQ=255)                                                                                                   gTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcg                                     <  2:1215412/140‑1 (MQ=255)                                                                                                    tCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcgc                                    <  1:344586/140‑1 (MQ=255)                                                                                                      tGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCgg                                  >  1:992907/1‑140 (MQ=255)                                                                                                       ggggTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGt                                 <  2:346223/140‑1 (MQ=255)                                                                                                         ggTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCg                               <  1:633431/140‑1 (MQ=255)                                                                                                         ggTGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCg                               <  2:759338/140‑1 (MQ=255)                                                                                                           tGCCTGTTCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCt                             <  1:726679/140‑1 (MQ=255)                                                                                                             ccTGTTCCGTCAACTGTGGTAGCCGCTGCGCGCTTCGTTTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTTACAATTCCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTcc                           <  2:997245/140‑1 (MQ=255)                                                                                                                 ttCCGTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATcc                       <  1:105822/140‑1 (MQ=255)                                                                                                                     gTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgccg                   <  2:599856/140‑1 (MQ=255)                                                                                                                     gTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgccg                   <  1:428492/140‑1 (MQ=255)                                                                                                                     gTCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgccg                   <  1:701441/140‑1 (MQ=255)                                                                                                                      tCAACTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCgg                  >  1:565592/1‑140 (MQ=255)                                                                                                                         aCTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt               >  1:1545719/1‑140 (MQ=255)                                                                                                                         aCTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt               <  2:942007/140‑1 (MQ=255)                                                                                                                         aCTGTGGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt               <  2:1185075/140‑1 (MQ=255)                                                                                                                           tgtgGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt             >  1:422563/1‑140 (MQ=255)                                                                                                                           tgtgGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt             <  1:1472013/140‑1 (MQ=255)                                                                                                                             tgGTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa           >  1:1316879/1‑140 (MQ=255)                                                                                                                             gggTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa           <  2:664146/139‑1 (MQ=255)                                                                                                                              ggTAGCCGCTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTaa          >  1:1215949/1‑140 (MQ=255)                                                                                                                                      cTGCGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCg  <  1:553017/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCG  >  NC_000913/1658124‑1658396

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