BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,028,548	G→A	D58D (GAC→GAT)	hspQ ←	heat shock protein, hemimethylated DNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,028,548	0	G	A	100.0%	1006.7 / NA	243	D58D (GAC→GAT)	hspQ	heat shock protein, hemimethylated DNA‑binding protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (130/113); total (130/113)

Only 100 of 246 total aligned reads displayed.

GCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGA  >  NC_000913/1028450‑1028645

gcCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTatcatc                                                                                                 >  1:1507314/1‑101 (MQ=255)

  cAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGACGTTATCAtcct                                                                                               >  1:1507316/1‑101 (MQ=255)

     cTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCa                                                                                            >  2:1507319/1‑101 (MQ=255)

      tCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCAt                                                                                           >  1:1507322/1‑101 (MQ=255)

      tCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCAt                                                                                           >  2:1507324/1‑101 (MQ=255)

         tCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCAc                                                                                        >  1:1507328/1‑101 (MQ=255)

         tCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCAc                                                                                        >  2:1507327/1‑101 (MQ=255)

          ccATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACt                                                                                       <  2:1506836/101‑1 (MQ=255)

           cATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTa                                                                                      <  1:1506860/101‑1 (MQ=255)

           cATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTa                                                                                      >  1:1507330/1‑101 (MQ=255)

           cATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTa                                                                                      >  1:1507333/1‑101 (MQ=255)

                atggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGa                                                                                 >  1:1507336/1‑101 (MQ=255)

                atggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGa                                                                                 <  1:1506837/101‑1 (MQ=255)

                  ggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATa                                                                               >  2:1507339/1‑101 (MQ=255)

                  ggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATa                                                                               <  1:1507093/101‑1 (MQ=255)

                   gCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATAc                                                                              >  1:1507341/1‑101 (MQ=255)

                   gCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATAc                                                                              <  1:1506878/101‑1 (MQ=255)

                   gCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATAc                                                                              <  1:1506782/101‑1 (MQ=255)

                      gTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCAc                                                                           >  2:1507344/1‑101 (MQ=255)

                      gTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCAc                                                                           >  1:1507345/1‑101 (MQ=255)

                        tCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACgg                                                                         >  2:1507348/1‑101 (MQ=255)

                         cAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGa                                                                        >  2:1507350/1‑101 (MQ=255)

                         cAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGa                                                                        <  2:1506839/101‑1 (MQ=255)

                          aGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAg                                                                       >  2:1507353/1‑101 (MQ=255)

                          aGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAg                                                                       <  2:1507169/101‑1 (MQ=255)

                           ggATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAgc                                                                      <  1:1507053/101‑1 (MQ=255)

                              tGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAgcggc                                                                   <  2:1506867/101‑1 (MQ=255)

                               gTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAgcggcg                                                                  <  2:1506932/101‑1 (MQ=255)

                               gTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAgcggcg                                                                  >  2:1507357/1‑101 (MQ=255)

                                 tCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGcgcg                                                                >  1:1507359/1‑101 (MQ=255)

                                   aTCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGaa                                                              <  1:1507192/101‑1 (MQ=255)

                                   aTCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGaa                                                              >  1:1507361/1‑101 (MQ=255)

                                    tCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAg                                                             >  2:1507364/1‑101 (MQ=255)

                                      cTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCt                                                           <  1:1506995/101‑1 (MQ=255)

                                      cTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCt                                                           <  1:1507065/101‑1 (MQ=255)

                                      cTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCt                                                           >  2:1507367/1‑101 (MQ=255)

                                       tGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTc                                                          <  2:1507108/101‑1 (MQ=255)

                                       tGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTc                                                          >  2:1507373/1‑101 (MQ=255)

                                       tGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTc                                                          >  1:1507372/1‑101 (MQ=255)

                                       tGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTc                                                          >  1:1507369/1‑101 (MQ=255)

                                        gCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCg                                                         <  1:1507004/101‑1 (MQ=255)

                                         cAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGt                                                        <  2:1506966/101‑1 (MQ=255)

                                          aaTTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTc                                                       <  2:1506951/101‑1 (MQ=255)

                                            ttCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCTTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTc                                                       <  2:1506750/99‑1 (MQ=255)

                                             tCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCAtt                                                    <  1:1506831/101‑1 (MQ=255)

                                             tCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCAtt                                                    <  2:1506993/101‑1 (MQ=255)

                                             tCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCAtt                                                    >  2:1507377/1‑101 (MQ=255)

                                              cGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTc                                                   <  1:1506937/101‑1 (MQ=255)

                                              cGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTc                                                   >  1:1507378/1‑101 (MQ=255)

                                                cTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCAc                                                 >  1:1507381/1‑101 (MQ=255)

                                                cTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCAc                                                 >  2:1507382/1‑101 (MQ=255)

                                                 tACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCAcc                                                <  1:1507032/101‑1 (MQ=255)

                                                   cTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGc                                              >  1:1507385/1‑101 (MQ=255)

                                                      aGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCCCGGAGCGGCGCGAAGCTCGTCATTCACCGCCAg                                           <  2:1506928/101‑1 (MQ=255)

                                                      aGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAg                                           >  1:1507389/1‑101 (MQ=255)

                                                        cTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCt                                         >  2:1507390/1‑101 (MQ=255)

                                                         tGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtc                                        <  2:1506640/101‑1 (MQ=255)

                                                          gggCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtca                                       >  1:1507393/1‑101 (MQ=255)

                                                           ggCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtcat                                      >  2:1507395/1‑101 (MQ=255)

                                                             cTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtcatca                                    <  1:1506902/101‑1 (MQ=255)

                                                             cTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtcatca                                    <  2:1506821/101‑1 (MQ=255)

                                                             cTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtcatca                                    >  2:1507398/1‑101 (MQ=255)

                                                                cAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGc                                 >  2:1507402/1‑101 (MQ=255)

                                                                cAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGc                                 >  1:1507401/1‑101 (MQ=255)

                                                                cAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGc                                 <  2:1507076/101‑1 (MQ=255)

                                                                cAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCAGCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGc                                 <  2:1506936/101‑1 (MQ=255)

                                                                 aGCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCg                                >  1:1507405/1‑101 (MQ=255)

                                                                  gCCAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAg                                   <  2:1506859/97‑1 (MQ=255)

                                                                    cAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACg                             >  1:1507410/1‑101 (MQ=255)

                                                                    cAGGTAGGTATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACg                             >  2:1507408/1‑101 (MQ=255)

                                                                     aggtaggtATGAACCGGTAGGCCGTTATCATCCTCCATCACTa                                                                                      <  2:1506721/43‑1 (MQ=255)

                                                                     aggtaggtATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACgg                            <  2:1506905/101‑1 (MQ=255)

                                                                     aggtaggtATGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACgg                            >  1:1507412/1‑101 (MQ=255)

                                                                       gtaggtaTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACg                             >  2:1507416/1‑98 (MQ=255)

                                                                         aggtaTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCa                        >  1:1507420/1‑101 (MQ=255)

                                                                         aggtaTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCa                        >  1:1507418/1‑101 (MQ=255)

                                                                          ggtaTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCag                       <  2:1506800/101‑1 (MQ=255)

                                                                           gtaTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCaga                      <  1:1506962/101‑1 (MQ=255)

                                                                            taTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGc                                              >  1:1507424/1‑76 (MQ=255)

                                                                            taTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                     <  1:1507001/101‑1 (MQ=255)

                                                                            taTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                     <  1:1507162/101‑1 (MQ=255)

                                                                             aTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagaga                    <  1:1507125/101‑1 (MQ=255)

                                                                             aTGAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagaga                    >  1:1507426/1‑101 (MQ=255)

                                                                               gAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGc                  <  1:1506799/101‑1 (MQ=255)

                                                                               gAACCGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGc                  >  2:1507428/1‑101 (MQ=255)

                                                                                  ccGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGaa               <  2:1507003/101‑1 (MQ=255)

                                                                                  ccGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGaa               >  1:1507430/1‑101 (MQ=255)

                                                                                   cGGTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAAt              <  2:1506970/101‑1 (MQ=255)

                                                                                    ggTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATa             <  2:1507198/101‑1 (MQ=255)

                                                                                     gTAGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATaa            >  2:1507433/1‑101 (MQ=255)

                                                                                       aGGCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAAc          <  1:1506917/101‑1 (MQ=255)

                                                                                        ggCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAAcc         <  1:1506650/101‑1 (MQ=255)

                                                                                         gCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCg        >  2:1507436/1‑101 (MQ=255)

                                                                                         gCCGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCg        <  2:1507057/101‑1 (MQ=255)

                                                                                          ccGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCgg       >  1:1507438/1‑101 (MQ=255)

                                                                                           cGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCggg      >  2:1507439/1‑101 (MQ=255)

                                                                                           cGTTATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCggg      <  2:1507025/101‑1 (MQ=255)

                                                                                             ttATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTc    >  1:1507442/1‑101 (MQ=255)

                                                                                             ttATCATCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTc    <  1:1507149/101‑1 (MQ=255)

                                                                                               atcatcCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGa  >  1:1507446/1‑101 (MQ=255)

GCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGA  >  NC_000913/1028450‑1028645

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