Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,059,964	C→T	*317* (TAG→TAA)	cbl ←	ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,059,964	0	C	T	100.0%	223.3 / NA	69	*317* (TAG→TAA)	cbl	ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/35);  total (34/35)

TAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTC  >  NC_000913/2059829‑2060093                                                                                                                                         |                                                                                                                                  tAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                <  2:1480866/139‑1 (MQ=255) tAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                <  2:626182/139‑1 (MQ=255)      tgtgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTg                                                                                                                           <  2:701211/139‑1 (MQ=255)          ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATaa                                                                                                                       >  1:1104367/1‑139 (MQ=255)          ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATaa                                                                                                                       >  1:1144251/1‑139 (MQ=255)          ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATaa                                                                                                                       >  2:77547/1‑139 (MQ=255)          ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATaa                                                                                                                       >  2:639627/1‑139 (MQ=255)               acaacGATTCCTCTGTTGTTTAGTCGGGAGAAAGGCGGGCTGTTAAACCGTATTTTTCTGGTTTGCGGACCGGCCGAGGGGGCCAAATTCTGGAAAAACCCCCTTTT‑TAGCGGGGTTTTTGTTttattttttttttata                                                                                                                  >  1:467198/1‑127 (MQ=255)               acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAAt                                                                                                                  <  1:77547/139‑1 (MQ=255)               acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAAt                                                                                                                  <  1:538335/139‑1 (MQ=255)                  acGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               <  2:1104367/139‑1 (MQ=255)                  acGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               <  1:243339/139‑1 (MQ=255)                        ccTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTc                                                                                                         >  1:75237/1‑139 (MQ=255)                            tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      <  1:1381664/122‑1 (MQ=255)                            tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      >  2:1381664/1‑122 (MQ=255)                                gTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCttt                                                                                                 <  2:1144251/139‑1 (MQ=255)                                 ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCTCTGGTTCGACTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCGAGCTTTTATGGCGGGAGTTTTGTTATATATGATAATCAATTTCCTCTTCACTGCTATc                                                                                                >  1:1304538/1‑139 (MQ=21)                                 ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                                >  1:1370767/1‑139 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCa                                                                                                                    <  2:432459/117‑1 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCa                                                                                                                    >  1:432459/1‑117 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCAc                                                                                                       >  1:1498240/1‑130 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCAc                                                                                                       <  2:1498240/130‑1 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCa                                                                                              >  1:1119264/1‑139 (MQ=255)                                   cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCa                                                                                              >  1:1180702/1‑139 (MQ=255)                                      tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                           >  2:1264591/1‑139 (MQ=255)                                        ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      >  1:1392357/1‑110 (MQ=255)                                        ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      <  2:1392357/110‑1 (MQ=255)                                        ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                         <  1:1264591/139‑1 (MQ=255)                                        ggTAGAACGACGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               <  2:829102/117‑1 (MQ=255)                                        ggTAGAACGACGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               >  1:829102/1‑117 (MQ=255)                                                  cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGa                                                                               >  1:610695/1‑139 (MQ=255)                                                  cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGC‑TTCCATCACCTGCCGCTTGat                                                                              <  1:639627/139‑1 (MQ=255)                                                   ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                              <  2:1370767/139‑1 (MQ=255)                                                    gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGata                                                                             >  2:54985/1‑139 (MQ=255)                                                               cGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                  <  1:1217131/139‑1 (MQ=255)                                                                 tATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTgg                                                                  <  1:95902/137‑2 (MQ=255)                                                                 tATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                  >  2:95902/1‑137 (MQ=255)                                                                 tATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACa                                                                >  1:1094438/1‑139 (MQ=255)                                                                   tGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGatat                                                                            >  2:1596126/1‑125 (MQ=255)                                                                   tGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGatat                                                                            <  1:1596126/125‑1 (MQ=255)                                                                   tGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGt                                                              >  1:985509/1‑139 (MQ=255)                                                                    gTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                             <  2:1180702/139‑1 (MQ=255)                                                                    gTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                             <  1:1314916/139‑1 (MQ=255)                                                                      cACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                                           <  2:1094438/139‑1 (MQ=255)                                                                       aCTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      >  1:1384617/1‑79 (MQ=255)                                                                       aCTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      <  2:1384617/79‑1 (MQ=255)                                                                        cTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                              >  1:712187/1‑118 (MQ=255)                                                                        cTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                              <  2:712187/118‑1 (MQ=255)                                                                          ggTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGCCCTGCAt                                                       <  1:1477332/139‑1 (MQ=255)                                                                           gTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                >  1:601160/1‑65 (MQ=255)                                                                           gTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                <  2:601160/65‑1 (MQ=255)                                                                           gTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCAtt                                                      >  1:389361/1‑139 (MQ=255)                                                                                aGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                >  1:367923/1‑60 (MQ=255)                                                                                aGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                <  2:367923/60‑1 (MQ=255)                                                                                     aGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCa                                                                                                        >  2:769342/1‑79 (MQ=255)                                                                                     aGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCa                                                                                                        <  1:769342/79‑1 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTg                                                                                                                           <  2:75280/59‑1 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTg                                                                                                                           >  1:75280/1‑59 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                    <  2:651106/82‑1 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                    >  1:651106/1‑82 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCa                                                        <  1:826907/126‑1 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCa                                                        >  2:826907/1‑126 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                   <  2:188075/131‑1 (MQ=255)                                                                                      gTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                   >  1:188075/1‑131 (MQ=255)                                                                                           a‑ggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                             >  1:1325694/1‑116 (MQ=255)                                                                                           a‑ggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                             <  2:1325694/116‑1 (MQ=255)                                                                                                                          ataGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCt        <  2:75237/139‑1 (MQ=255)                                                                                                                               gggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                   <  2:1172745/107‑1 (MQ=255)                                                                                                                               gggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                   >  1:1172745/1‑107 (MQ=255)                                                                                                                                ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  2:421722/1‑139 (MQ=255)                                                                                                                                ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  2:1192464/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                  TAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGA‑GGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTC  >  NC_000913/2059829‑2060093

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