Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,439	C→T	A738T (GCA→ACA)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,439	0	C	T	100.0%	263.2 / NA	69	A738T (GCA→ACA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/34);  total (35/34)

GTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTG  >  NC_000913/186295‑186586                                                                                                                                                 |                                                                                                                                                    gTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTg                                                                                                                                               >  1:112579/1‑151 (MQ=255)   ggCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGa                                                                                                                                              >  1:316402/1‑150 (MQ=255)             tACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAAtt                                                                                                                                   >  2:437436/1‑151 (MQ=255)              aCTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAAtt                                                                                                                                   <  2:289519/150‑1 (MQ=255)              aCTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTg                                                                                                                                  >  1:576672/1‑151 (MQ=255)               cattGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgc                                                                                                                                 <  2:154437/149‑1 (MQ=255)                   gCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcggcg                                                                                                                             <  2:451081/151‑1 (MQ=255)                   gCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcggcg                                                                                                                             <  2:45406/151‑1 (MQ=255)                   gCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcggcg                                                                                                                             <  1:993722/151‑1 (MQ=255)                    cGTCAGTACTTGCTCCAGACCAAACCGAATAACCGCATGACGATCTGCGGACAGCGGGTTTCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGGGAAAATTTGTGTGTCTTGAACACTTAAATTGcggcgg                                                                                                                            <  2:531664/151‑1 (MQ=255)                         gTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTa                                                                                                                        >  1:124498/1‑150 (MQ=255)                           aCTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAAtt                                                                                                                     >  1:230762/1‑151 (MQ=255)                              tGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTc                                                                                                                    <  1:2546/149‑1 (MQ=255)                              tGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTc                                                                                                                    <  1:590359/149‑1 (MQ=255)                                 tCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcac                                                                                                               >  1:253302/1‑151 (MQ=255)                                  ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcaca                                                                                                              <  1:49372/151‑1 (MQ=255)                                  ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGCTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGCCTAATTCGGcac                                                                                                               >  2:693935/1‑150 (MQ=255)                                     gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAc                                                                                                           <  1:1022067/151‑1 (MQ=255)                                      aCCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAg                                                                                                             >  1:334873/1‑148 (MQ=255)                                      aCCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcg                                                                                                          <  1:437462/151‑1 (MQ=255)                                      aCCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcg                                                                                                          >  1:661817/1‑151 (MQ=255)                                        cAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcggc                                                                                                        <  1:731089/151‑1 (MQ=255)                                         aaaCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATTTTCATACCGTTGCGAGTGGTGAAAATTTGTGTGTCGTGAAGACTTAAATTGCTGCCGTCCAATTCGCCACAGAcggc                                                                                                        <  2:123476/150‑1 (MQ=255)                                                  aaCCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAAc                                                                                                 >  1:805088/1‑148 (MQ=255)                                                       cATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAg                                                                                          <  2:576640/150‑1 (MQ=255)                                                           aCGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCgg                                                                                     >  1:991928/1‑151 (MQ=255)                                                            cGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCgg                                                                                     >  1:828882/1‑150 (MQ=255)                                                            cGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGt                                                                                    >  1:591613/1‑151 (MQ=255)                                                               tCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTcc                                                                                  >  1:796571/1‑150 (MQ=255)                                                               tCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTcc                                                                                  >  1:365234/1‑150 (MQ=255)                                                                 tGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCggg                                                                               >  1:653218/1‑151 (MQ=255)                                                                  gCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGt                                                                                    <  1:377526/145‑1 (MQ=255)                                                                  gCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTc                                                                                   <  2:620666/146‑1 (MQ=255)                                                                  gCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCgg                                                                                <  1:992058/149‑1 (MQ=255)                                                                  gCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCg                                                                                 <  2:860632/148‑1 (MQ=255)                                                                      aCAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCa                                                                          <  1:1069150/151‑1 (MQ=255)                                                                       cAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTcc                                                                           >  1:493326/1‑149 (MQ=255)                                                                        aGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTa                                                                                                                                       <  1:232986/88‑1 (MQ=255)                                                                        aGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTa                                                                                                                                       >  2:232961/1‑88 (MQ=255)                                                                           gggTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAAGGCGGTCCGGGCTCCAGATa                                                                      >  2:650475/1‑150 (MQ=255)                                                                                  cATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATCCCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGGGGCAAACAGATAAAGGCGGTCCGGGCTACAGATAAAAAtctc                                                              >  2:218454/1‑151 (MQ=255)                                                                                        gTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGGTAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGcc                                                        >  1:917970/1‑151 (MQ=255)                                                                                        gTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGcc                                                        >  2:124874/1‑151 (MQ=255)                                                                                         ttCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGcc                                                        >  2:926705/1‑150 (MQ=255)                                                                                           ccAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTc                                                      <  1:365315/150‑1 (MQ=255)                                                                                             aGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAc                                                   >  1:543036/1‑151 (MQ=255)                                                                                                aCGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAcgcg                                                <  1:434157/151‑1 (MQ=255)                                                                                                aCGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAcgcg                                                >  1:381567/1‑151 (MQ=255)                                                                                                 cGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGt                                               <  2:175475/151‑1 (MQ=255)                                                                                                   aTAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAg                                             <  2:492815/151‑1 (MQ=255)                                                                                                      aaGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCt                                          >  1:362771/1‑151 (MQ=255)                                                                                                                tCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTa                                  >  2:54887/1‑149 (MQ=255)                                                                                                                  gCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGGGGTGTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGACTTCACGCGTTGCCTGCGGGCTAAGCa                              <  2:633504/151‑1 (MQ=255)                                                                                                                  gCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCa                              <  1:166945/151‑1 (MQ=255)                                                                                                                   ccATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCCGGCCCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGc                               >  2:443323/1‑149 (MQ=255)                                                                                                                       aCCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATAc                          <  2:974803/150‑1 (MQ=255)                                                                                                                       aCCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATAc                          <  2:509779/150‑1 (MQ=255)                                                                                                                        ccGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCa                        <  1:794442/151‑1 (MQ=255)                                                                                                                             gcgAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCgg                    >  1:408705/1‑150 (MQ=255)                                                                                                                             gcgAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGt                   <  1:438635/151‑1 (MQ=255)                                                                                                                               gAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGtt                  <  2:83171/150‑1 (MQ=255)                                                                                                                                 gtggtgAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGc               <  1:829796/151‑1 (MQ=255)                                                                                                                                  tggtgAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGGAGCCTGCGGGCTAAGCAATACCAGCGGTTTg                >  2:192049/1‑149 (MQ=255)                                                                                                                                   ggtgAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCt              <  2:238663/150‑1 (MQ=255)                                                                                                                                     tgAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTa            >  1:415989/1‑150 (MQ=255)                                                                                                                                       aaaaTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAAt         <  2:424364/151‑1 (MQ=255)                                                                                                                                        aaaTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATc        >  2:568460/1‑151 (MQ=255)                                                                                                                                          aTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTTCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTaaa          >  2:480419/1‑147 (MQ=255)                                                                                                                                           tttGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGGGCCTCCACGCGTAGCCTGCGGGCTAAGCAATAACAGCGGTTTGCTTAAATCATg     >  2:853911/1‑151 (MQ=255)                                                                                                                                             tgtgtgTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGc    <  1:189109/150‑1 (MQ=255)                                                                                                                                              gtgtgtCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTg  <  1:61930/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTG  >  NC_000913/186295‑186586

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