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%	234.1 / NA	65	K448K (AAA→AAG)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/31);  total (34/31)

TATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACA  >  NC_000913/1873251‑1873520                                                                                                                                      |                                                                                                                                         tattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTg                                                                                                                                    >  2:157562/1‑140 (MQ=255)        tCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACata                                                                                                                             >  1:338150/1‑140 (MQ=255)        tCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACata                                                                                                                             >  1:318155/1‑140 (MQ=255)            gCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgc                                                                                                                         >  1:89306/1‑140 (MQ=255)             cGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgcg                                                                                                                        <  1:408360/140‑1 (MQ=255)                ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     <  2:89306/140‑1 (MQ=255)                ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     <  1:157562/140‑1 (MQ=255)                ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     <  2:102922/140‑1 (MQ=255)                ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                     >  1:837438/1‑140 (MQ=255)                        tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAaga                                                                                                             >  2:702155/1‑140 (MQ=255)                         cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAagag                                                                                                            >  1:563608/1‑140 (MQ=255)                          gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAgaga                                                                                                           <  2:750959/140‑1 (MQ=255)                          gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAgaga                                                                                                           >  2:238531/1‑140 (MQ=255)                               gTAAAAGTGATTATGCCTTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGc                                                                                                      >  2:329724/1‑140 (MQ=255)                                 aaaaGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt                                                                                                    <  1:329724/140‑1 (MQ=255)                                 aaaaGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCGGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGTTCGGCtt                                                                                                    <  2:691377/140‑1 (MQ=255)                                  aaaGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc                                                                                                   >  1:274809/1‑140 (MQ=255)                                       gATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTc                                                                                              <  1:702155/140‑1 (MQ=255)                                       gATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTc                                                                                              <  1:111023/140‑1 (MQ=255)                                        aTTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc                                                                                             >  2:514725/1‑140 (MQ=255)                                        aTTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc                                                                                             <  1:374494/140‑1 (MQ=255)                                        aTTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc                                                                                             >  2:138111/1‑140 (MQ=255)                                         ttATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCg                                                                                            >  1:912323/1‑140 (MQ=255)                                                aTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTa                                                                                     <  1:138111/140‑1 (MQ=255)                                                 tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                    >  2:45824/1‑140 (MQ=255)                                                 tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                    >  2:579708/1‑140 (MQ=255)                                                       tCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATg                                                                              >  1:148105/1‑140 (MQ=255)                                                            ggCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATgaaaga                                                                         >  1:62979/1‑140 (MQ=255)                                                             gCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATgaaagaa                                                                        <  1:45824/140‑1 (MQ=255)                                                               gATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGaaaa                                                                      <  2:614746/140‑1 (MQ=255)                                                                aTGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAAc                                                                     <  1:181221/140‑1 (MQ=255)                                                                 tGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACg                                                                    <  2:148105/140‑1 (MQ=255)                                                                 tGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACg                                                                    <  2:912323/140‑1 (MQ=255)                                                                 tGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACg                                                                    <  1:238531/140‑1 (MQ=255)                                                                   aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  <  2:361832/140‑1 (MQ=255)                                                                   aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                                  <  1:630037/140‑1 (MQ=255)                                                                     ttCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGaaaa                                                                      <  1:558066/134‑1 (MQ=255)                                                                     ttCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGaaaa                                                                      >  2:558066/1‑134 (MQ=255)                                                                     ttCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATAc                                                                >  1:415700/1‑140 (MQ=255)                                                                              atTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatga                                                       <  2:270585/140‑1 (MQ=255)                                                                               ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                      <  2:62979/140‑1 (MQ=255)                                                                               ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                      >  1:899989/1‑140 (MQ=255)                                                                               ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                      <  2:563608/140‑1 (MQ=255)                                                                                  ttGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCg                                                   >  1:346131/1‑140 (MQ=255)                                                                                    gTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGta                                                 >  2:533695/1‑140 (MQ=255)                                                                                         ttCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAg                                            >  2:48712/1‑140 (MQ=255)                                                                                         ttCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAg                                            >  1:747451/1‑140 (MQ=255)                                                                                          tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATaagct                                           >  1:226450/1‑137 (MQ=255)                                                                                          tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGc                                           >  1:237173/1‑140 (MQ=255)                                                                                          tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGc                                           >  2:860529/1‑140 (MQ=255)                                                                                          tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCACAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGc                                           >  2:763885/1‑140 (MQ=255)                                                                                                 cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                    <  1:533695/140‑1 (MQ=255)                                                                                                 cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                    <  2:747451/140‑1 (MQ=255)                                                                                                   cAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGa                                  <  1:603421/140‑1 (MQ=255)                                                                                                   cAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCATATAAAGCTTCCGATGa                                  >  2:650633/1‑140 (MQ=255)                                                                                                    aaaTTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAg                                 <  2:522744/140‑1 (MQ=255)                                                                                                      aTTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCACAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCg                               <  1:763885/140‑1 (MQ=255)                                                                                                       ttGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGt                              <  2:318155/140‑1 (MQ=255)                                                                                                           agcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtat                          <  2:338150/140‑1 (MQ=255)                                                                                                            gcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtata                         >  2:180930/1‑140 (MQ=255)                                                                                                                aaCTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                     >  1:719260/1‑140 (MQ=255)                                                                                                                 aCTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa                    <  2:632264/140‑1 (MQ=255)                                                                                                                         aaCGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAATACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtata                         <  1:371826/127‑1 (MQ=255)                                                                                                                         aaCGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtata                         >  2:371826/1‑127 (MQ=255)                                                                                                                               tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa      >  1:166571/1‑140 (MQ=255)                                                                                                                               tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa      >  1:12543/1‑140 (MQ=255)                                                                                                                               tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCATATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa      >  1:392598/1‑140 (MQ=255)                                                                                                                                   aaagCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACa  <  2:274809/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         TATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACA  >  NC_000913/1873251‑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: ‑