Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,873,384	A→G	K448K (AAA→AAG)	dgcJ →	putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,873,384	0	A	G	100.0%	238.5 / NA	66	K448K (AAA→AAG)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/34);  total (32/34)

CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACA  >  NC_000913/1873248‑1873520                                                                                                                                         |                                                                                                                                         cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                       <  1:642854/140‑1 (MQ=255) cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                       <  2:720824/140‑1 (MQ=255) cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                       <  1:422212/140‑1 (MQ=255)        aGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAAc                                                                                                                                >  2:997837/1‑140 (MQ=255)            cAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACatat                                                                                                                            >  2:215711/1‑140 (MQ=255)            cAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACatat                                                                                                                            >  2:608382/1‑140 (MQ=255)               gCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgc                                                                                                                         <  1:221006/140‑1 (MQ=255)               gCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgc                                                                                                                         >  2:550856/1‑140 (MQ=255)                 gATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCgcgc                                                                                                                       <  1:27551/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     <  1:550856/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     <  2:714954/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     >  2:929157/1‑140 (MQ=255)                    tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccgc                                                                                                                    >  2:883007/1‑140 (MQ=255)                    tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccgc                                                                                                                    <  1:608382/140‑1 (MQ=255)                       aCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGa                                                                                                                 <  1:929157/140‑1 (MQ=255)                           tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATTAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAaga                                                                                                             >  2:245856/1‑140 (MQ=255)                           tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAaga                                                                                                             >  2:233640/1‑140 (MQ=255)                             gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATTAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAgaga                                                                                                           <  1:245856/140‑1 (MQ=255)                                  gTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGTACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGc                                                                                                      <  1:390238/140‑1 (MQ=255)                                  gTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGc                                                                                                      <  1:678970/140‑1 (MQ=255)                                           aTTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc                                                                                             >  1:951643/1‑140 (MQ=255)                                              aTGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGt                                                                                          >  2:239872/1‑140 (MQ=255)                                                   aTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGGCGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACGTATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTa                                                                                     >  2:741278/1‑140 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                    >  1:800277/1‑140 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                    >  1:758589/1‑140 (MQ=255)                                                     ccGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAcg                                                                                   <  1:883007/140‑1 (MQ=255)                                                     ccGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAcg                                                                                   >  1:762493/1‑140 (MQ=255)                                                         cTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGAt                                                                               >  2:478070/1‑140 (MQ=255)                                                          tCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATg                                                                              >  2:481017/1‑140 (MQ=255)                                                                  gATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGaaaa                                                                      <  1:545784/140‑1 (MQ=255)                                                                    tGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACg                                                                    <  1:239872/140‑1 (MQ=255)                                                                    tGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACg                                                                    <  1:595321/140‑1 (MQ=255)                                                                      aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  <  1:336499/140‑1 (MQ=255)                                                                      aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  <  2:800277/140‑1 (MQ=255)                                                                       aTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATa                                                                 <  1:62215/140‑1 (MQ=255)                                                                         tCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACg                                                               <  2:758589/140‑1 (MQ=255)                                                                           tGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  <  1:597553/135‑1 (MQ=255)                                                                           tGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  >  2:597553/1‑135 (MQ=255)                                                                                catTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatg                                                        >  2:355547/1‑140 (MQ=255)                                                                                 atTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatga                                                       <  1:816454/140‑1 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGCAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACat                                                         >  1:329597/1‑137 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                      >  1:1076124/1‑140 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                      <  1:233640/140‑1 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACat                                                         <  2:329597/137‑1 (MQ=255)                                                                                       gTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGta                                                 <  1:534479/140‑1 (MQ=255)                                                                                       gTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGta                                                 >  2:10456/1‑140 (MQ=255)                                                                                            ttCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAg                                            >  2:676192/1‑140 (MQ=255)                                                                                            ttCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAg                                            >  1:732359/1‑140 (MQ=255)                                                                                             tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGc                                           >  1:390890/1‑140 (MQ=255)                                                                                               gACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCtt                                         >  1:89206/1‑140 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                    <  2:23485/140‑1 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                    <  1:997837/140‑1 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                    <  1:481017/140‑1 (MQ=255)                                                                                                      cAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGa                                  <  1:478070/140‑1 (MQ=255)                                                                                                        aaTTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGc                                <  2:367434/140‑1 (MQ=255)                                                                                                           tGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGc                                >  1:760731/1‑137 (MQ=255)                                                                                                           tGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGc                                <  2:760731/137‑1 (MQ=255)                                                                                                            gcagcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGttt                            <  2:390890/140‑1 (MQ=255)                                                                                                                cacaACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTatat                        >  1:256027/1‑140 (MQ=255)                                                                                                                     cTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa                   >  1:879888/1‑140 (MQ=255)                                                                                                                        ccTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATaa                >  1:507204/1‑140 (MQ=255)                                                                                                                            aaCGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg            <  2:89206/140‑1 (MQ=255)                                                                                                                                 aTCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa       >  1:582642/1‑140 (MQ=255)                                                                                                                                 aTCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTATATATGTCAATAAGCAGAACaa       >  2:999407/1‑140 (MQ=255)                                                                                                                                  tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa      >  2:1003352/1‑140 (MQ=255)                                                                                                                                    gAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaaaa    <  2:419522/140‑1 (MQ=255)                                                                                                                                      aaagCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACa  <  2:418279/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACA  >  NC_000913/1873248‑1873520

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