Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,690,933	G→A	P559S (CCT→TCT)	bcsB ←	cellulose synthase periplasmic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,690,933	0	G	A	100.0%	229.9 / NA	65	P559S (CCT→TCT)	bcsB	cellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/36);  total (29/36)

ATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGGCGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGATGGTGGAG  >  NC_000913/3690797‑3691063                                                                                                                                         |                                                                                                                                   gtgatCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGacgc                                                                                                                                 >  1:529921/2‑140 (MQ=255)    atCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCttt                                                                                                                              <  2:69667/140‑1 (MQ=255)      cATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCg                                                                                                                            <  2:529921/140‑1 (MQ=255)      cATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCg                                                                                                                            >  2:191357/1‑140 (MQ=255)          atgTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCAt                                                                                                                        <  1:655127/140‑1 (MQ=255)            gTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCa                                                                                                                      >  1:99455/1‑140 (MQ=255)            gTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGCCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCa                                                                                                                      >  2:21386/1‑140 (MQ=255)             tCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCAc                                                                                                                     >  2:206278/1‑140 (MQ=255)               ggCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACgg                                                                                                                   >  1:208666/1‑140 (MQ=255)               ggCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACgg                                                                                                                   >  2:492890/1‑140 (MQ=255)               ggCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACgg                                                                                                                   >  1:401650/1‑140 (MQ=255)                gCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGt                                                                                                                  >  2:34112/1‑140 (MQ=255)                gCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGt                                                                                                                  >  1:488582/1‑140 (MQ=255)                    cTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATg                                                                                                              >  2:109395/1‑140 (MQ=255)                      ttGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGt                                                                                                            >  2:804774/1‑140 (MQ=255)                       tGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGtt                                                                                                           <  1:30541/140‑1 (MQ=255)                           cTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCg                                                                                                       <  2:187579/140‑1 (MQ=255)                              aaTGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACa                                                                                                    <  2:99455/140‑1 (MQ=255)                              aaTGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACa                                                                                                    >  1:619154/1‑140 (MQ=255)                                   tgctgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATcc                                                                                               >  2:548843/1‑140 (MQ=255)                                   tgctgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATcc                                                                                               >  1:507134/1‑140 (MQ=255)                                    gctgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCg                                                                                              <  1:592853/140‑1 (MQ=255)                                    gctgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCg                                                                                              <  2:619154/140‑1 (MQ=255)                                     ctgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGc                                                                                             >  2:559011/1‑140 (MQ=255)                                      tgcCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGcc                                                                                            >  1:116149/1‑140 (MQ=255)                                        ccATCCACGGTCACCGCCAGATTAATCTCCGGGAAGCCCCTCTGTGCGCCGATTAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTATACGCCTTCCGCATCCCGGTGATGGTTTGCGACAGATCCGCCAt                                                                                          >  2:548644/1‑140 (MQ=255)                                        ccATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCAt                                                                                          >  2:474539/1‑140 (MQ=255)                                        ccATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCAt                                                                                          >  2:353486/1‑140 (MQ=255)                                        ccATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCAt                                                                                          <  1:917144/140‑1 (MQ=255)                                        ccATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCAt                                                                                          >  2:219380/1‑140 (MQ=255)                                        ccATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGCTTGCGACAGATCCGCCAt                                                                                          >  1:488355/1‑140 (MQ=255)                                               cGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTg                                                                                   <  2:202193/140‑1 (MQ=255)                                                     ccGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATggg                                                                             <  1:353486/140‑1 (MQ=255)                                                         cAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGc                                                                         >  2:599875/1‑140 (MQ=255)                                                            aTTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCcg                                                                      >  1:906478/1‑140 (MQ=255)                                                             ttAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCcgc                                                                     <  2:358505/140‑1 (MQ=255)                                                                aTCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGtt                                                                  >  2:852742/1‑140 (MQ=255)                                                                    ccGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCa                                                              <  2:116149/140‑1 (MQ=255)                                                                      gggAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCaaa                                                            >  2:653000/1‑140 (MQ=255)                                                                      gggAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCaaa                                                            <  1:852742/140‑1 (MQ=255)                                                                       ggAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAg                                                           <  1:206278/140‑1 (MQ=255)                                                                        gAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAgg                                                          <  1:804774/140‑1 (MQ=255)                                                                                  tgtgCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGCTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATc                                                <  2:488355/140‑1 (MQ=255)                                                                                       gcCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCa                                           <  1:548843/140‑1 (MQ=255)                                                                                       gcCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCa                                           <  1:559011/140‑1 (MQ=255)                                                                                          gATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCg                                        <  1:824671/140‑1 (MQ=255)                                                                                             aaaaCCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGa                                     <  1:21386/140‑1 (MQ=255)                                                                                              aaaCCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGAt                                    <  1:653000/140‑1 (MQ=255)                                                                                                aCCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGa                                  <  2:401650/140‑1 (MQ=255)                                                                                                  cAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAg                                <  1:548644/140‑1 (MQ=255)                                                                                                     cAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTgg                             >  1:745989/1‑140 (MQ=255)                                                                                                     cAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTgg                             <  1:34112/140‑1 (MQ=255)                                                                                                      aGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGt                            <  2:906478/140‑1 (MQ=255)                                                                                                      aGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGt                            <  1:109395/140‑1 (MQ=255)                                                                                                      aGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGt                            <  2:488582/140‑1 (MQ=255)                                                                                                      aGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGt                            <  1:492890/140‑1 (MQ=255)                                                                                                             aGCAACGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACtt                     <  2:785234/140‑1 (MQ=255)                                                                                                                 aCGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCgag                 <  1:599875/140‑1 (MQ=255)                                                                                                                  cGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCgaga                >  1:212896/1‑140 (MQ=255)                                                                                                                  cGTTTCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCgaga                >  1:681805/1‑140 (MQ=255)                                                                                                                     ttCCATCTGTGCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGt             <  2:507134/140‑1 (MQ=255)                                                                                                                            tgtgCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGAtggt      <  2:212896/140‑1 (MQ=255)                                                                                                                             gtgCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGAtggtg     <  2:681805/140‑1 (MQ=255)                                                                                                                               gCTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGATGGTGGa   <  1:474539/140‑1 (MQ=255)                                                                                                                                cTTCGTTAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGATGGTGGAg  >  1:688162/1‑140 (MQ=255)                                                                                                                                     ttAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGAtggtg     >  2:706930/1‑132 (MQ=255)                                                                                                                                     ttAGACGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGAtggtg     <  1:706930/132‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   ATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTTCCATCTGTGCTTCGTTAGGCGCTTTCGGCATCACGGTGATGGTTTGCGACAGATCCGCCATCCGGCTGAATGGGAAGCCCGCGTTAGCAAAGGCGCGTAGATCCGGCATCGGGATGAAGTGGTAATACTTCGAGAAGTCGATGGTGGAG  >  NC_000913/3690797‑3691063

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