Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,407,078	T→A	T240S (ACT→TCT)	dgcM ←	diguanylate cyclase DgcM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407,078	0	T	A	100.0%	163.0 / NA	52	T240S (ACT→TCT)	dgcM	diguanylate cyclase DgcM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (26/26);  total (26/26)

GCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGC  >  NC_000913/1406833‑1407318                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 gcAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAatat                                                                                                                                                                                                                                             >  1:11949/1‑251 (MQ=255)   aaCTTTCGAGGGGGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCa                                                                                                                                                                                                                                           <  2:183064/251‑1 (MQ=255)           aGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATACTGcacact                                                                                                                                                                                                                                   >  2:377441/1‑249 (MQ=255)                     cGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGAACAATACATAACATg                                                                                                                                                                                                                         >  1:192352/1‑251 (MQ=255)                         ggCGTATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTGTAAAACGATCGGTGTCGATGAGCAACAAGCTTTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACGGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGAAATATCATGCACAATACATCACATGAGCt                                                                                                                                                                                                                     <  1:433109/251‑1 (MQ=255)                           cgcATAACACCTCATCACCTTTAGAATGCCAATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTg                                                                                                                                                                                                                   >  1:156874/1‑251 (MQ=255)                                 acacCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGCTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGcc                                                                                                                                                                                                              >  2:242778/1‑250 (MQ=255)                                          caCCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTc                                                                                                                                                                                                     >  1:122229/1‑250 (MQ=255)                                             cTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCCTTTGACCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCCTCGCGTCATGATAAGCAGCATGTTCCCGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCactc                                                                                                                                                                                                 >  2:139581/1‑248 (MQ=255)                                                              ataTCGTGAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCGGTAATCCTGAGCGATGGGGAGATGCGGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTg                                                                                                                                                                                <  2:129420/249‑1 (MQ=255)                                                                     tAATGTTTTTAAAACGACCGTGGTCGAGGAGCAACAAGCTGTAACCCTGAGCGATGGCGAGATGCTGCATGTGGCCTGGTTCCGTAATGTGATAAAACGGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATg                                                                                                                                                                          <  2:192352/250‑1 (MQ=255)                                                                           ggttAAAACGGTCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATCCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTc                                                                                                                                                                   <  2:403735/249‑1 (MQ=255)                                                                               aaaaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAAc                                                                                                                                                               >  1:453944/1‑251 (MQ=255)                                                                               aaaaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAAc                                                                                                                                                               >  1:426547/1‑251 (MQ=255)                                                                                    gATCGGTGTCGATGATCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTAATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATcc                                                                                                                                                          <  2:274548/251‑1 (MQ=255)                                                                                      tCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGTCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGc                                                                                                                                                        >  1:295925/1‑251 (MQ=255)                                                                                            gcGATGACCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                                                                                                                                                  <  2:86571/250‑1 (MQ=255)                                                                                             cGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGACCCATCCGCCAGttt                                                                                                                                                  >  2:328516/1‑250 (MQ=255)                                                                                              gATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTAt                                                                                                                                                >  2:163271/1‑251 (MQ=255)                                                                                                tGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACAATCCGCCAGTTTATg                                                                                                                                               >  1:516449/1‑250 (MQ=255)                                                                                                    caacaaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACa                                                                                                                                           >  1:364081/1‑250 (MQ=255)                                                                                                    caacaaGCTGTAATCCTGAGCGATGGCGAGATGCTACATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACa                                                                                                                                           >  1:2379/1‑250 (MQ=255)                                                                                                     aacagGCTGTAATCCTGTGCGATGGAGAGAGGCTGCATTGGCCCTGGTTCCGTAATGTGATAAAACGGTCGCCTATTCAGTACTCCGGTCACCGCGTCATGGTGAGCATAATGTCCCAGTTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACaaa                                                                                                                                         <  2:156874/251‑1 (MQ=255)                                                                                                       caaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCATGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGATGCTCCTCCAGCCGTTTTTGCTCAGAAATATCCTGCACAATACATAAAATGCGATTGTCGCCATAAATTTCAATCGGTCCAGCATAGGTATGAAGATGACGCGTCGAAACATCCGCCATTTTATGAACaaac                                                                                                                                        >  2:358293/1‑249 (MQ=255)                                                                                                          gCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTc                                                                                                                                     >  1:215897/1‑250 (MQ=255)                                                                                                                agccggAGCGGTGGCGAGATGCTGCATTTGGCCTGGTGCCGGAATGTAATAAAACGGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGtt                                                                                                                              <  1:265629/246‑1 (MQ=255)                                                                                                                           gggCGAGATGGTGCATTTGGCCTGGTTCCGAAATGTGATAAAACGGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCg                                                                                                                   <  1:183069/250‑1 (MQ=255)                                                                                                                            ggCGAGATGCTGCATTCGGCCTGGTTCCGTAATGTGATAAAACTGCCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCg                                                                                                                   <  1:97510/250‑1 (MQ=255)                                                                                                                                atatgctgtatttgtctgggTTCCGTAATGTGATAAAACTGTCTCCGATTCATTAATCCGTTCATCTCGTCATGGTGAGCAGAATATTCCAGCAGCTCCTCCATCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATTAGCTTGTCGCCATAAATTTCAATCGGTCCTGCATAGGTCTGCACATGACGAGTCGAACCAGCCGCCAGTTTATGAACAAAATTCAAAGTTTTATGACCACCGTGTaa                                                                                                              <  2:488527/233‑1 (MQ=255)                                                                                                                                            tgtccTGTTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAAGCTGGTCATCTCTGCATGGGTAGCAGCATGTTCCTACGGCTCCTCCAGCCGTTTTTGCTCAGAAATATCAGGCTCAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGTTCTGCACATGACGATACGAACCATCCGTCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCCAGAttt                                                                                                   <  2:503213/247‑1 (MQ=255)                                                                                                                                             ggCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTc                                                                                                  >  2:471589/1‑250 (MQ=255)                                                                                                                                                 tGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCCCAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgc                                                                                              >  1:34234/1‑250 (MQ=255)                                                                                                                                                         tAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGc                                                                                     >  1:204634/1‑251 (MQ=255)                                                                                                                                                               gATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATgac                                                                                <  1:50262/250‑1 (MQ=255)                                                                                                                                                                   aaaCTGCCGCCGATTCAGTAATCCGGTGATCGCGTCATGGTGAGCAGCATGTCCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATgacgcgac                                                                           <  2:364081/251‑1 (MQ=255)                                                                                                                                                                          cGCCGATTCAGTAATCCGGGCATTGCGTCATGGTGAGCAGCATGTACCAGCTGCTCCTCCAGCCGTTGTTTCTCAGAAATATCATGCACAATACATAACAGGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGGCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTGTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAg                                                                    <  1:335961/251‑1 (MQ=255)                                                                                                                                                                               aTTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCatat                                                               >  2:188840/1‑251 (MQ=255)                                                                                                                                                                                 tCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATTTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCAtatt                                                              <  1:377441/250‑1 (MQ=255)                                                                                                                                                                                 tCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCAtattt                                                             >  1:345719/1‑251 (MQ=255)                                                                                                                                                                                                          gTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGcaca                                    <  1:401280/251‑1 (MQ=255)                                                                                                                                                                                                            gAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATAAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATc                                  <  1:210790/251‑1 (MQ=255)                                                                                                                                                                                                                agcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGttt                              >  1:241330/1‑251 (MQ=255)                                                                                                                                                                                                                 gcatatcccatctgtccctccAGCCGTTTGTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGTTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTc                             <  2:11949/235‑1 (MQ=255)                                                                                                                                                                                                                   aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCAt                           <  1:95415/251‑1 (MQ=255)                                                                                                                                                                                                                      ttCCAGCTGCTCCTCAAGCCGTTTTTGCTCAGAAATATCATGAACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGCCGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGt                        >  1:418038/1‑251 (MQ=255)                                                                                                                                                                                                                           gctgctCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTACATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTtata                    >  1:81716/1‑250 (MQ=255)                                                                                                                                                                                                                              gctCCTCCAGCCGTTTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATAGTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCa                <  1:30329/251‑1 (MQ=255)                                                                                                                                                                                                                                     cAGCCGTGTTTGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGGCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAAt         <  1:251356/251‑1 (MQ=255)                                                                                                                                                                                                                                          gTTTTTGCTCTGAAATTTCATGCAAAATAAATAACATGAGCGTGTCGCCATAAAGTTCAATCGGTCCGGCATAGGTCGGGACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAATGTTTATGCCCTCCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCTCGAGCATATTTATTTCCCAGGTATGTTCCTGGCACATCGTTTCATGGTTATATCCATAGAAATTGAgc    <  2:322126/251‑1 (MQ=255)                                                                                                                                                                                                                                             tttGCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAgcgc  <  1:328516/250‑1 (MQ=255)                                                                                                                                                                                                                                                gCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCt                                          <  2:189552/207‑1 (MQ=255)                                                                                                                                                                                                                                                gCTCAGAAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCt                                          >  1:189552/1‑207 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGC  >  NC_000913/1406833‑1407318

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