Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	876,549	T→C	S405S (AGT→AGC)	dgcI →	putative diguanylate cyclase DgcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	876,549	0	T	C	100.0%	218.8 / NA	62	S405S (AGT→AGC)	dgcI	putative diguanylate cyclase DgcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/37);  total (25/37)

GTAGATAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCAT  >  NC_000913/876412‑876667                                                                                                                                          |                                                                                                                       gTAGATAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCAt                                                                                                                      <  2:1049214/140‑1 (MQ=255) gTAGATAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCAt                                                                                                                      <  1:894858/140‑1 (MQ=255)     aTAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCAt                                                                                                                      >  1:173212/1‑136 (MQ=255)     aTAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGgt                                                                                                                  <  2:1060523/140‑1 (MQ=255)          ggTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGc                                                                                                             <  2:253702/140‑1 (MQ=255)                 tCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCa                                                                                                      >  2:246979/1‑140 (MQ=255)                  cTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCaa                                                                                                     <  2:855981/140‑1 (MQ=255)                    aCCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAaat                                                                                                   >  2:39973/1‑140 (MQ=255)                           cggcggTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTc                                                                                            >  1:762524/1‑140 (MQ=255)                           cggcggTGACGAATTTGCGGGGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTc                                                                                            >  2:897049/1‑140 (MQ=255)                            ggcggTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTcc                                                                                           <  1:39973/140‑1 (MQ=255)                             gcggTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCt                                                                                          <  2:369210/140‑1 (MQ=255)                                   aCGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAaatgaa                                                                                                >  1:1133829/1‑128 (MQ=255)                                   aCGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAaatgaa                                                                                                <  2:1133829/128‑1 (MQ=255)                                   aCGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTc                                                                                    <  1:897049/140‑1 (MQ=255)                                           gCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                            >  1:384709/1‑140 (MQ=255)                                                 aTTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAAt                                                                      >  1:422273/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATTGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATGGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  1:738684/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGaaa                                                                        <  2:148813/128‑1 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGaaa                                                                        >  1:148813/1‑128 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  1:559999/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  1:650604/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  2:802430/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  2:737895/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  2:686011/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  1:356045/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  2:423035/1‑140 (MQ=255)                                                           gTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACACTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                            >  1:288104/1‑140 (MQ=255)                                                            tGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATg                                                           >  2:357302/1‑140 (MQ=255)                                                              gAACGCTGGGAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAg                                                         <  2:470613/140‑1 (MQ=255)                                                              gAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGTAAATGCTGGCGGATGAg                                                         >  1:733012/1‑140 (MQ=255)                                                              gAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAg                                                         <  2:212632/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  2:384709/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  2:337062/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  1:569336/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  2:559999/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  2:650604/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  2:922845/140‑1 (MQ=255)                                                               aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGc                                                        <  1:357302/140‑1 (MQ=255)                                                                aCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCg                                                       >  2:1074105/1‑140 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGTGGATGAGCgg                                                      <  1:906910/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  2:96578/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  1:623153/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  1:1074105/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  1:52332/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  1:93141/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  2:356045/140‑1 (MQ=255)                                                                 cGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACACTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCgg                                                      <  2:288104/140‑1 (MQ=255)                                                                    tGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTTAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTg                                                   >  2:663826/1‑140 (MQ=255)                                                                     ggAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGt                                                  >  1:672092/1‑140 (MQ=255)                                                                       aaGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGtat                                                <  1:967584/140‑1 (MQ=255)                                                                                 aTCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGa                                      >  2:638570/1‑140 (MQ=255)                                                                                    gATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaa                                   >  1:182424/1‑140 (MQ=255)                                                                                          aaaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCa                             <  2:425107/140‑1 (MQ=255)                                                                                             aGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAggg                          <  2:182424/140‑1 (MQ=255)                                                                                              gAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGc                         <  1:802430/140‑1 (MQ=255)                                                                                               aaaTTGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGca                        <  2:738684/140‑1 (MQ=255)                                                                                               aaaTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGTAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGca                        <  2:733012/140‑1 (MQ=255)                                                                                                aaTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGcac                       >  1:160795/1‑140 (MQ=255)                                                                                                aaTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGcac                       >  1:963036/1‑140 (MQ=255)                                                                                                        aTTTCCAGGTGGGAAGCATTAATTTTAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc               <  1:663826/140‑1 (MQ=255)                                                                                                        aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc               <  1:686011/140‑1 (MQ=255)                                                                                                        aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc               <  1:423035/140‑1 (MQ=255)                                                                                                                     aaGCATTAATTTAAGTACCAGCATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCat  <  1:638570/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                       GTAGATAAAGGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCAT  >  NC_000913/876412‑876667

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