BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	933,473	A→G	T343A (ACG→GCG)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	933,473	0	A	G	100.0%	360.9 / NA	114	T343A (ACG→GCG)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (60/54); total (60/54)

Only 100 of 124 total aligned reads displayed.

GGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTT  >  NC_000913/933383‑933562

ggTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGcg                                                                                         >  2:2809068/1‑93 (MQ=255)

ggTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGGTGAACCTGTGACTCAGGcg                                                                                         >  2:3013172/1‑93 (MQ=255)

 gTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGgcgc                                                                                        >  1:121010/1‑93 (MQ=255)

    cgcCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGcgccgc                                                                                     <  1:4220776/93‑1 (MQ=255)

    cgcCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGcgccgc                                                                                     >  2:1528978/1‑93 (MQ=255)

     gcgaaTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCgccgcc                                                                                    >  2:3525031/4‑93 (MQ=255)

      ccAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCt                                                                                   >  1:16337/1‑93 (MQ=255)

      ccAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCt                                                                                   <  2:3295630/93‑1 (MQ=255)

       cAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTg                                                                                  >  2:1359333/1‑93 (MQ=255)

       cAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTg                                                                                  >  2:3817961/1‑93 (MQ=255)

        aaTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGt                                                                                 >  2:504850/1‑93 (MQ=255)

        aaTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGGTGAACCTGTGACTCAGGCGCCGCCTGt                                                                                 >  2:3948034/1‑93 (MQ=255)

         aTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGtt                                                                                <  1:2333556/93‑1 (MQ=255)

         aTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGtt                                                                                >  2:100725/1‑93 (MQ=255)

         aTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGACTGtt                                                                                >  2:3572889/1‑93 (MQ=255)

           tACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGc                                                                              <  1:3556068/93‑1 (MQ=255)

            aCCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGcc                                                                             >  2:823308/1‑93 (MQ=255)

             ccGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCct                                                                            <  2:2847658/93‑1 (MQ=255)

              cGACCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctc                                                                           <  2:1656288/93‑1 (MQ=255)

              cGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctc                                                                           >  2:1174331/1‑93 (MQ=255)

              cGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctc                                                                           >  1:1570036/1‑93 (MQ=255)

               gAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctct                                                                          <  1:1273353/93‑1 (MQ=255)

               gAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctct                                                                          <  1:4641079/93‑1 (MQ=255)

               gAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCctct                                                                          >  1:963161/1‑93 (MQ=255)

                caCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTg                                                                         <  1:840893/92‑1 (MQ=255)

                aaCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTg                                                                         <  1:1733399/93‑1 (MQ=255)

                aaCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTg                                                                         >  1:1682614/1‑93 (MQ=255)

                aaCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTg                                                                         >  1:3408940/1‑93 (MQ=255)

                 aCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGt                                                                        >  2:103546/1‑93 (MQ=255)

                 aCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGACTGTTGCCTCTGt                                                                        >  2:800552/1‑93 (MQ=255)

                  ccTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGtt                                                                       >  2:1456159/1‑93 (MQ=255)

                  ccTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGtt                                                                       <  2:156263/93‑1 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  1:2323737/1‑93 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  1:2717874/1‑93 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  1:4397220/1‑93 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  1:3861039/1‑93 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      <  1:3698802/93‑1 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  1:1363027/1‑93 (MQ=255)

                   cTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTAAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  2:4596498/1‑93 (MQ=255)

                   cTGGCGCTGTAGAAGCTGCTGCTACGACGGCGACACAAAGCTTCGGTGCGCCGCGTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTg                                                                      >  2:4199173/1‑93 (MQ=255)

                    tGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGa                                                                     >  1:2942400/1‑93 (MQ=255)

                    tGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGa                                                                     >  2:1395470/1‑93 (MQ=255)

                     gTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGAt                                                                    <  1:4584149/93‑1 (MQ=255)

                      tCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATg                                                                   >  1:119892/1‑93 (MQ=255)

                       cGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGt                                                                  <  1:3013381/93‑1 (MQ=255)

                        gCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGtt                                                                 >  1:3012481/1‑93 (MQ=255)

                           gTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCa                                                              >  2:1552349/1‑93 (MQ=255)

                           gTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCa                                                              >  2:1359953/1‑93 (MQ=255)

                            tAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCAc                                                             >  2:4511199/1‑93 (MQ=255)

                              gcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCAccc                                                           >  2:1778729/1‑92 (MQ=255)

                              gcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCt                                                           <  1:2804801/93‑1 (MQ=255)

                              gcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCt                                                           <  1:3538776/93‑1 (MQ=255)

                                agcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTgc                                                         >  2:1445064/1‑93 (MQ=255)

                                agcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTgc                                                         <  2:2325940/93‑1 (MQ=255)

                                     ctgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAAc                                                    <  1:3588322/93‑1 (MQ=255)

                                     ctgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAAc                                                    <  2:1785650/93‑1 (MQ=255)

                                     ctgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAAc                                                    <  1:2760431/93‑1 (MQ=255)

                                      tgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGCTGTTCCACCTGCGCAAcc                                                   >  2:4213773/1‑93 (MQ=255)

                                      tgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAAcc                                                   <  2:2716444/93‑1 (MQ=255)

                                      tgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAAcc                                                   >  2:2507382/1‑93 (MQ=255)

                                       gctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCt                                                  >  2:3807577/1‑93 (MQ=255)

                                        ctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTa                                                 >  2:296043/1‑93 (MQ=255)

                                         tACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTAc                                                >  2:2527394/1‑93 (MQ=255)

                                          aCCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACa                                               <  1:3828607/93‑1 (MQ=255)

                                           ccACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAg                                              <  2:1304627/93‑1 (MQ=255)

                                           ccACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAg                                              >  2:2541547/1‑93 (MQ=255)

                                             aCGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTa                                            <  1:4477772/93‑1 (MQ=255)

                                             aCGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTa                                            >  1:2087723/1‑93 (MQ=255)

                                               ggCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGc                                          <  1:482564/93‑1 (MQ=255)

                                               ggCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGc                                          >  2:3100650/1‑93 (MQ=255)

                                                gCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGcc                                         <  2:538346/93‑1 (MQ=255)

                                                 cGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCt                                        >  1:2691582/1‑93 (MQ=255)

                                                      caAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAg                                   <  1:3585248/93‑1 (MQ=255)

                                                      caAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAg                                   <  1:3996576/93‑1 (MQ=255)

                                                         aGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCt                                >  2:2860878/1‑93 (MQ=255)

                                                         aGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCt                                >  2:3018518/1‑93 (MQ=255)

                                                          gCTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGGAGCCTg                               >  2:2121302/1‑93 (MQ=255)

                                                           cTGGGCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGt                              >  1:3692602/1‑93 (MQ=255)

                                                             gggCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTAc                            >  1:1080084/1‑93 (MQ=255)

                                                             gggCTGCGCCGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTAc                            <  2:1828795/93‑1 (MQ=255)

                                                                      cGGTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCcac                   <  2:3110710/93‑1 (MQ=255)

                                                                       ggTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCcaca                  <  1:254886/93‑1 (MQ=255)

                                                                       ggTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCcaca                  >  2:1633638/1‑93 (MQ=255)

                                                                        gTTGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACaa                 >  2:2126787/1‑93 (MQ=255)

                                                                         ttGAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACaaa                <  1:2831771/93‑1 (MQ=255)

                                                                           gAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACg              <  2:1165345/93‑1 (MQ=255)

                                                                           gAACCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACg              >  2:794860/1‑93 (MQ=255)

                                                                            aaCCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACgg             <  2:1119940/93‑1 (MQ=255)

                                                                            aaCCTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACgg             >  2:4647666/1‑93 (MQ=255)

                                                                              ccTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGga           <  2:3988321/93‑1 (MQ=255)

                                                                              ccTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGga           >  2:3502891/1‑93 (MQ=255)

                                                                              ccTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGga           <  2:81049/93‑1 (MQ=255)

                                                                               cTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGgag          >  1:2204482/1‑93 (MQ=255)

                                                                               cTGTGACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGgag          >  1:2893656/1‑93 (MQ=255)

                                                                                tgtgACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGgaga         <  1:538239/93‑1 (MQ=255)

                                                                                 gtgACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGagag        <  1:3423204/93‑1 (MQ=255)

                                                                                   gACTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGcc      <  2:1314813/93‑1 (MQ=255)

                                                                                     cTCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCgg    >  2:1543763/1‑93 (MQ=255)

                                                                                      tCAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt   >  2:4205254/1‑93 (MQ=255)

                                                                                       cAGGCGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGtt  <  2:3256192/93‑1 (MQ=255)

GGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTT  >  NC_000913/933383‑933562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`