Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,076,627	G→C	S161T (AGT→ACT)	yihY →	BrkB family putative transporter, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,076,627	0	G	C	100.0%	277.4 / NA	77	S161T (AGT→ACT)	yihY	BrkB family putative transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (39/38);  total (39/38)

GATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCCGCGTAC  >  NC_000913/4076493‑4076761                                                                                                                                       |                                                                                                                                       gATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTcc                                                                                                                                    >  1:961815/1‑139 (MQ=255)      aCTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTAttt                                                                                                                               >  1:533888/1‑139 (MQ=255)       cTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTg                                                                                                                              >  2:101380/1‑139 (MQ=255)        tCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  2:616870/139‑1 (MQ=255)        tCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  2:1110881/139‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                            >  1:934767/1‑139 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             >  2:144434/1‑138 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  1:144434/138‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCACTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                            >  2:294102/1‑139 (MQ=255)            tCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCtctc                                                                                                                         >  2:1135210/1‑139 (MQ=255)              gATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTccc                                                                                                                       <  2:357502/139‑1 (MQ=255)                  gcgcGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   <  2:961815/139‑1 (MQ=255)                  gcgcGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   >  2:264409/1‑139 (MQ=255)                   cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcgc                                                                                                                  >  2:996752/1‑139 (MQ=255)                   cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcgc                                                                                                                  >  2:382273/1‑139 (MQ=255)                    gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCt                                                                                                                 >  2:405798/1‑139 (MQ=255)                     cgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTg                                                                                                                >  2:8619/1‑139 (MQ=255)                          aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcga                                                                                                           <  1:101380/139‑1 (MQ=255)                            tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                                         <  2:628886/139‑1 (MQ=255)                            tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                                         <  2:457435/139‑1 (MQ=255)                               cATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                                      <  2:934767/139‑1 (MQ=255)                               cATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                                      <  2:21991/139‑1 (MQ=255)                                aTCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtc                                                                                                     >  2:759822/1‑139 (MQ=255)                                 tCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                              >  2:1060150/1‑129 (MQ=255)                                 tCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                              <  1:1060150/129‑1 (MQ=255)                                 tCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtct                                                                                                    <  2:533888/139‑1 (MQ=255)                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCaa                                                                                                 >  1:679131/1‑139 (MQ=255)                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCaa                                                                                                 >  2:431571/1‑139 (MQ=255)                                      cgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                               <  1:996752/139‑1 (MQ=255)                                        cAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACt                                                                                             <  1:8619/139‑1 (MQ=255)                                            aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             >  2:315690/1‑103 (MQ=255)                                            aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  1:315690/103‑1 (MQ=255)                                             aaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCAt                                                                                        <  2:905565/139‑1 (MQ=255)                                                     cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACa                                                                                   >  1:266982/1‑136 (MQ=255)                                                     cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACa                                                                                   <  2:266982/136‑1 (MQ=255)                                                     cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                                >  2:231983/1‑139 (MQ=255)                                                     cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                                >  1:638294/1‑139 (MQ=255)                                                       aCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  1:915199/92‑1 (MQ=255)                                                       aCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGAGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCAACAACGtg                                                                              >  2:812282/1‑139 (MQ=255)                                                       aCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTCTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             >  2:915199/1‑92 (MQ=255)                                                        ccAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgc                                                                             <  1:1016171/139‑1 (MQ=255)                                                                      gTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCg                                                               <  1:1135210/139‑1 (MQ=255)                                                                                     gatgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACtt                                                                                                                                      >  2:1150511/1‑53 (MQ=255)                                                                                       tgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAt                                              <  1:657598/139‑1 (MQ=255)                                                                                        gatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             <  1:106690/59‑1 (MQ=255)                                                                                        gatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                             >  2:106690/1‑59 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtctc                                           >  2:119760/1‑139 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtctc                                           >  1:1033344/1‑139 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtctc                                           >  1:181264/1‑139 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtctc                                           >  2:972057/1‑139 (MQ=255)                                                                                                  cTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                                   <  1:382273/139‑1 (MQ=255)                                                                                                  cTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                                   <  1:972057/139‑1 (MQ=255)                                                                                                   tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                          >  1:1061135/1‑67 (MQ=255)                                                                                                   tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                          <  2:1061135/67‑1 (MQ=255)                                                                                                         cGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTCCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTa                            <  1:812282/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGCCAGTCTGTCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa                       <  2:1033344/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa                       <  1:119760/139‑1 (MQ=255)                                                                                                                ggCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCAtt                     >  2:564926/1‑139 (MQ=255)                                                                                                                ggCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCAtt                     >  2:423199/1‑139 (MQ=255)                                                                                                                 gCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGtt                                                             >  2:494745/1‑98 (MQ=255)                                                                                                                 gCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGtt                                                             <  1:494745/98‑1 (MQ=255)                                                                                                                 gCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTg                    <  1:264409/139‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   <  1:3603/37‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   >  2:3603/1‑37 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   >  1:1269499/1‑37 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                                   <  2:1269499/37‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                                   >  2:496131/1‑117 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                                   <  1:496131/117‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAc             <  1:405798/139‑1 (MQ=255)                                                                                                                         cAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAcc            <  2:441712/139‑1 (MQ=255)                                                                                                                         cAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAcc            <  1:164467/139‑1 (MQ=255)                                                                                                                           gTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                               >  1:152583/1‑54 (MQ=255)                                                                                                                           gTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                               <  2:152583/54‑1 (MQ=255)                                                                                                                            tCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATc         >  2:792228/1‑139 (MQ=255)                                                                                                                               ggCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                                         <  1:436908/40‑1 (MQ=255)                                                                                                                               ggCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                                         >  2:436908/1‑40 (MQ=255)                                                                                                                                   aTCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCCGCGTAc  >  2:268643/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       GATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCCGCGTAC  >  NC_000913/4076493‑4076761

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