Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,146,009	A→G	D915G (GAT→GGT)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,146,009	0	A	G	100.0%	231.4 / NA	64	D915G (GAT→GGT)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/28);  total (36/28)

TCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGT  >  NC_000913/2145872‑2146144                                                                                                                                          |                                                                                                                                        tCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTc                                                                                                                                       >  1:270524/1‑140 (MQ=255) tCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTc                                                                                                                                       >  1:336384/1‑140 (MQ=255)             gCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGcc                                                                                                                           <  2:336384/140‑1 (MQ=255)             gCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGcc                                                                                                                           <  2:536070/140‑1 (MQ=255)             gCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGcc                                                                                                                           >  2:102729/1‑140 (MQ=255)               tCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTGTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCAt                                                                                                                         >  2:955285/1‑140 (MQ=255)                  cccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGcc                                                                                                                           >  2:895105/1‑135 (MQ=255)                  cccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGcc                                                                                                                           <  1:895105/135‑1 (MQ=255)                      cGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTg                                                                                                                  >  1:1045066/1‑140 (MQ=255)                      cGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTg                                                                                                                  >  1:949283/1‑140 (MQ=255)                           tCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGAccgc                                                                                                             >  1:693062/1‑140 (MQ=255)                           tCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGAccgc                                                                                                             >  1:783041/1‑140 (MQ=255)                           tCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGAccgc                                                                                                             >  1:877790/1‑140 (MQ=255)                                      cGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCa                                                                                                  <  1:672230/140‑1 (MQ=255)                                         aTCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGa                                                                                               <  2:111606/140‑1 (MQ=255)                                                   cgcgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                                                     >  2:525823/1‑140 (MQ=255)                                                   cgcgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                                                     >  2:684413/1‑140 (MQ=255)                                                     cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                                                                   >  2:357923/1‑140 (MQ=255)                                                     cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                                                                   >  1:617766/1‑140 (MQ=255)                                                     cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                                                                   >  1:625831/1‑140 (MQ=255)                                                     cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                                                                   >  1:700072/1‑140 (MQ=255)                                                     cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                                                                   <  2:936891/140‑1 (MQ=255)                                                        aaTTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGc                                                                                <  1:684413/140‑1 (MQ=255)                                                             gtgGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCaa                                                                           <  2:11556/140‑1 (MQ=255)                                                                gCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAg                                                                        <  2:783041/140‑1 (MQ=255)                                                                gCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAg                                                                        <  1:310585/140‑1 (MQ=255)                                                                gCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAg                                                                        >  2:611655/1‑140 (MQ=255)                                                                  tGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCg                                                                      <  2:877790/140‑1 (MQ=255)                                                                     tttCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                                                                   <  1:955285/140‑1 (MQ=255)                                                                         aCTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCa                                                               >  2:706502/1‑140 (MQ=255)                                                                           ttAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGt                                                             <  2:753934/140‑1 (MQ=255)                                                                            tAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                                                            <  2:700072/140‑1 (MQ=255)                                                                              aGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTaaa                                                          >  2:806098/1‑140 (MQ=255)                                                                               gCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAg                                                         >  2:231430/1‑140 (MQ=255)                                                                                         tGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt                                               <  2:949283/140‑1 (MQ=255)                                                                                         tGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt                                               <  1:193572/140‑1 (MQ=255)                                                                                         tGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt                                               <  2:1045066/140‑1 (MQ=255)                                                                                          gCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtc                                              >  1:235879/1‑140 (MQ=255)                                                                                          gCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtc                                              >  2:353935/1‑140 (MQ=255)                                                                                          gCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtc                                              >  2:849255/1‑140 (MQ=255)                                                                                           cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCa                                                      <  2:88379/131‑1 (MQ=255)                                                                                           cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCa                                                      >  1:88379/1‑131 (MQ=255)                                                                                               ggCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctccct                                         >  2:895503/1‑140 (MQ=255)                                                                                                gCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctccctc                                        >  1:745436/1‑140 (MQ=255)                                                                                                          tgatgaACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCg                              >  1:595161/1‑140 (MQ=255)                                                                                                              gaACAAACATTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGc                          >  2:437033/1‑140 (MQ=255)                                                                                                                 caaacaTTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCgg                       >  2:82907/1‑140 (MQ=255)                                                                                                                  aaacaTTTCGTCGTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGt                      <  1:231430/140‑1 (MQ=255)                                                                                                                       tttCGTCGTAGCTTCAGTGGTCCGGCGCTTGGCCATCCTCTTGACCGCCGAGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCCCCCCCTTTCCGTTGCCGGTTTGAg                 <  2:689128/140‑1 (MQ=255)                                                                                                                         tcgtcgTAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTa               >  1:421757/1‑140 (MQ=255)                                                                                                                          cgtcgtAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAg              <  2:274953/140‑1 (MQ=255)                                                                                                                          cgtcgtAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAg              >  1:253557/1‑140 (MQ=255)                                                                                                                          cgtcgtAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAg              >  2:757412/1‑140 (MQ=255)                                                                                                                            tcgtAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAgcg            >  2:211134/1‑140 (MQ=255)                                                                                                                               tAGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCa         <  1:895503/140‑1 (MQ=255)                                                                                                                                aGCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCAc        >  2:819945/1‑140 (MQ=255)                                                                                                                                 gCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACg       <  1:737333/140‑1 (MQ=255)                                                                                                                                 gCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACg       <  1:819945/140‑1 (MQ=255)                                                                                                                                 gCTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACg       <  2:235879/140‑1 (MQ=255)                                                                                                                                  cTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGc      >  2:575836/1‑140 (MQ=255)                                                                                                                                  cTTCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGc      <  2:266744/140‑1 (MQ=255)                                                                                                                                   ttCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGTCAGTTTGAGTAGCGCCACGCt     <  2:745436/140‑1 (MQ=255)                                                                                                                                   ttCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCt     <  1:437033/140‑1 (MQ=255)                                                                                                                                    tCAGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTg    <  2:253557/140‑1 (MQ=255)                                                                                                                                      aGCGGTCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGt  >  1:468672/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        TCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGT  >  NC_000913/2145872‑2146144

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