Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,628,895	(C)5→4	coding (44/2244 nt)	pdeF ←	cyclic di‑GMP phosphodiesterase PdeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,891	0	C	.	100.0%	277.1 / NA	64	coding (48/2244 nt)	pdeF	cyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (32/32);  total (32/32)

CCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAGCCCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACT  >  NC_000913/2628752‑2629025                                                                                                                                            |                                                                                                                                       cccAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAGc                                                                                                                                        >  2:148595/1‑140 (MQ=255)  ccAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAGcc                                                                                                                                       <  1:144084/140‑1 (MQ=255)         aaaaaCATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCAt                                                                                                                               <  2:892359/140‑1 (MQ=255)                tcatcatGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCAc                                                                                                                        <  2:493929/140‑1 (MQ=255)                 catcatGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACg                                                                                                                       <  1:426487/140‑1 (MQ=255)                 catcatGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACg                                                                                                                       <  2:756536/140‑1 (MQ=255)                             catcaGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTtatata                                                                                                           >  2:200433/1‑140 (MQ=255)                              atcaGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGCCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATaa                                                                                                          >  1:344231/1‑140 (MQ=255)                                caGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGt                                                                                                        >  1:920865/1‑140 (MQ=255)                                     ccAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCAt                                                                                                   >  1:13441/1‑140 (MQ=255)                                     ccAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCAt                                                                                                   >  1:797607/1‑140 (MQ=255)                                     ccAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCAt                                                                                                   >  1:519813/1‑140 (MQ=255)                                            cAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGttt                                                                                            <  1:375497/140‑1 (MQ=255)                                             aGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTc                                                                                           >  2:377653/1‑140 (MQ=255)                                             aGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTc                                                                                           >  1:1067029/1‑140 (MQ=255)                                              gATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCa                                                                                          >  1:431433/1‑140 (MQ=255)                                               aTAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCAt                                                                                         <  2:819954/140‑1 (MQ=255)                                                    agagAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATg                                                                                    >  1:472668/1‑140 (MQ=255)                                                      agaAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCt                                                                                  >  1:1089614/1‑140 (MQ=255)                                                        aaaaaCCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAg                                                                                       <  2:888031/133‑1 (MQ=255)                                                        aaaaaCCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAg                                                                                       >  1:888031/1‑133 (MQ=255)                                                           aaCCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTcagca                                                                             <  1:200433/140‑1 (MQ=255)                                                             ccTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGa                                                                           <  1:103203/140‑1 (MQ=255)                                                                             atGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGa                                                           <  2:797607/140‑1 (MQ=255)                                                                                  cgAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAaggag                                                      <  2:1067029/140‑1 (MQ=255)                                                                                  cgAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAaggag                                                      <  2:13441/140‑1 (MQ=255)                                                                                  cgAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAaggag                                                      <  1:377653/140‑1 (MQ=255)                                                                                   gAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAaggagg                                                     <  1:148595/140‑1 (MQ=255)                                                                                    aaaGTATTAATGTGGGCGAAAATGGCCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGgt                                                    <  2:344231/140‑1 (MQ=255)                                                                                      aGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGgtgt                                                  >  2:629903/1‑140 (MQ=255)                                                                                               gtgGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACa                                         >  2:958191/1‑140 (MQ=255)                                                                                                tgGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAg                                        >  1:514675/1‑140 (MQ=255)                                                                                                 gggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAag                                                         <  1:907990/122‑1 (MQ=255)                                                                                                 gggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAag                                                         >  2:907990/1‑122 (MQ=255)                                                                                                 gggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAgg                                       <  2:920865/140‑1 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGAc                                          <  2:997552/136‑1 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGAc                                          >  1:997552/1‑136 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGc                                      >  1:492677/1‑140 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGc                                      >  1:993145/1‑140 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGc                                      >  2:210170/1‑140 (MQ=255)                                                                                                  ggCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGc                                      >  1:750537/1‑140 (MQ=255)                                                                                                    cGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCAt                                    >  1:335615/1‑140 (MQ=255)                                                                                                        aaTGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGAATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCa                                >  2:817370/1‑140 (MQ=255)                                                                                                          tGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCAcc                              <  1:59502/140‑1 (MQ=255)                                                                                                             gCAAGATTAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGa                           >  1:691154/1‑140 (MQ=255)                                                                                                                   ttAAAGAAGGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGa                     >  2:416129/1‑140 (MQ=255)                                                                                                                        gaagGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAatgat                <  1:416129/140‑1 (MQ=255)                                                                                                                        gaagGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAatgat                <  1:629903/140‑1 (MQ=255)                                                                                                                          agGCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAAtgatga              <  2:1089614/140‑1 (MQ=255)                                                                                                                           ggCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATgatgat             >  2:969347/1‑140 (MQ=255)                                                                                                                           ggCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATgatgat             >  1:311100/1‑140 (MQ=255)                                                                                                                           ggCAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATgatgat             >  2:1032434/1‑140 (MQ=255)                                                                                                                             cAGGAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGatta           >  1:671335/1‑140 (MQ=255)                                                                                                                               ggAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGattatt         <  2:993145/140‑1 (MQ=255)                                                                                                                               ggAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGattatt         <  2:519813/140‑1 (MQ=255)                                                                                                                               ggAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGattatt         <  1:542229/140‑1 (MQ=255)                                                                                                                               ggAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGattatt         <  2:472668/140‑1 (MQ=255)                                                                                                                               ggAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGattatt         <  2:335615/140‑1 (MQ=255)                                                                                                                                gAACAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtc        >  1:352625/1‑140 (MQ=255)                                                                                                                                 aaCAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtca       <  2:691154/140‑1 (MQ=255)                                                                                                                                 aaCAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtca       <  2:625585/140‑1 (MQ=255)                                                                                                                                 aaCAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtca       <  2:311100/140‑1 (MQ=255)                                                                                                                                 aaCAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtca       <  1:210170/140‑1 (MQ=255)                                                                                                                                  aCAGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtcat      <  1:969347/140‑1 (MQ=255)                                                                                                                                    aGCGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATtcatca    <  1:958191/140‑1 (MQ=255)                                                                                                                                      cGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACt  >  1:444022/1‑140 (MQ=255)                                                                                                                                      cGGAAG‑CCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACt  >  1:114573/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                       CCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGGAACAGCGGAAGCCCCCACCATTTATCACGTATTTTTATATAAGTTGCATTCAGTTTCATAGATGCTCAGCAGAATCCCCCACATCCTGAAGGAGGTGTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACT  >  NC_000913/2628752‑2629025

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: ‑