BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,408	A→G	intergenic (+117/‑28)	b4251 → / → b4252	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,408	0	A	G	100.0%	1063.6 / NA	302	intergenic (+117/‑28)	b4251/b4252	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (149/153); total (149/153)

Only 100 of 308 total aligned reads displayed.

AAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAT  >  NC_000913/4472318‑4472497

aaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCg                                                                                         <  2:188391/93‑1 (MQ=255)

 aGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgt                                                                                        >  1:2768816/1‑93 (MQ=255)

 aGCCCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgt                                                                                        <  1:4719280/93‑1 (MQ=255)

  gcgcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgtt                                                                                       >  1:4595125/1‑93 (MQ=255)

    gcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                     >  1:664767/1‑93 (MQ=255)

    gcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                     >  1:3064878/1‑93 (MQ=255)

     ccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgtt                                                                                    <  1:3205149/93‑1 (MQ=255)

     ccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgtt                                                                                    <  1:2745247/93‑1 (MQ=255)

      cTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGttt                                                                                   >  2:3902688/1‑93 (MQ=255)

       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  2:6475171/1‑93 (MQ=255)

       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  2:3459471/1‑93 (MQ=255)

       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  1:891760/1‑93 (MQ=255)

         gcGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaa                                                                                <  2:148401/92‑1 (MQ=255)

         ccGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaa                                                                                >  1:1157715/1‑93 (MQ=255)

          cGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaat                                                                               <  1:1782918/93‑1 (MQ=255)

          cGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaat                                                                               >  1:4386724/1‑93 (MQ=255)

            aTAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                             >  2:1804775/1‑93 (MQ=255)

            aTAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                             <  1:1558938/93‑1 (MQ=255)

            aTAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                             >  1:592193/1‑93 (MQ=255)

             tAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaattaa                                                                            <  2:6340211/93‑1 (MQ=255)

              aaGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAg                                                                           >  2:2803743/1‑93 (MQ=255)

              aaGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAg                                                                           >  1:2596349/1‑93 (MQ=255)

              aaGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAg                                                                           >  2:1152556/1‑93 (MQ=255)

               aGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                          >  2:2758893/1‑93 (MQ=255)

                gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         >  1:820938/1‑93 (MQ=255)

                gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         <  1:5803296/93‑1 (MQ=255)

                gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         >  1:4830237/1‑93 (MQ=255)

                attGTGATAAGACAATAAGTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         <  1:195775/92‑1 (MQ=255)

                  tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGt                                                                       <  2:3183094/93‑1 (MQ=255)

                  tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGt                                                                       <  1:4575897/93‑1 (MQ=255)

                   gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                                      <  2:3715162/93‑1 (MQ=255)

                    tgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                                     <  2:2216073/93‑1 (MQ=255)

                    tgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                                     <  1:4105121/93‑1 (MQ=255)

                    tgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                                     <  1:142540/93‑1 (MQ=255)

                      aTAAGACAATAATTCAAGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaa                                                                   >  1:4763175/1‑93 (MQ=255)

                       tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaa                                                                  >  1:5230377/1‑93 (MQ=255)

                        aaGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                 <  2:4260835/93‑1 (MQ=255)

                            cAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                             >  2:4959606/1‑93 (MQ=255)

                             aataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatga                                                            <  2:1366476/93‑1 (MQ=255)

                              ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                                           >  1:3425880/1‑93 (MQ=255)

                              ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                                           >  1:5064872/1‑93 (MQ=255)

                               taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatc                                                          <  1:2987308/93‑1 (MQ=255)

                                aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatca                                                         <  1:5030462/93‑1 (MQ=255)

                                aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatca                                                         <  2:3327415/93‑1 (MQ=255)

                                 atTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                                        <  2:5652558/93‑1 (MQ=255)

                                  ttCACGCGTTAAGTCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                       <  1:513402/93‑1 (MQ=255)

                                  ttCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                       <  1:977501/93‑1 (MQ=255)

                                   tCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCg                                                      <  2:2396569/93‑1 (MQ=255)

                                    cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCgg                                                     >  2:5270839/1‑93 (MQ=255)

                                      cgcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaa                                                   >  1:5873651/1‑93 (MQ=255)

                                       gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                  <  1:6055195/93‑1 (MQ=255)

                                        cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                                                 >  2:3423079/1‑93 (MQ=255)

                                         gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATGGGGAatt                                                >  1:2477668/1‑92 (MQ=255)

                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               <  1:1741134/93‑1 (MQ=255)

                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               >  1:6211952/1‑93 (MQ=255)

                                           tAAGGCTAGCGAAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATAtt                                              <  1:6315916/93‑1 (MQ=255)

                                            aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                             >  2:4218799/1‑93 (MQ=255)

                                            aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                             >  2:5496943/1‑93 (MQ=255)

                                             aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                            >  2:5825604/1‑93 (MQ=255)

                                             aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                            <  2:6165337/93‑1 (MQ=255)

                                              ggCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                                           >  1:6448588/1‑93 (MQ=255)

                                              ggCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                                           >  1:1187847/1‑93 (MQ=255)

                                               gCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                                          <  2:5692101/93‑1 (MQ=255)

                                                cTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                         <  2:609946/93‑1 (MQ=255)

                                                cTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                         <  2:5941624/93‑1 (MQ=255)

                                                   gCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                                      >  2:3684604/1‑93 (MQ=255)

                                                   gCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                                      >  1:5466096/1‑93 (MQ=255)

                                                      gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                                   >  2:263748/1‑93 (MQ=255)

                                                       aaTTTATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAg                                  <  2:4043366/93‑1 (MQ=255)

                                                       aaTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAg                                  <  2:1742587/93‑1 (MQ=255)

                                                         ttGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                >  1:1017865/1‑93 (MQ=255)

                                                         ttGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                <  2:2325504/93‑1 (MQ=255)

                                                         ttGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                >  2:1239772/1‑93 (MQ=255)

                                                           taTTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                              <  2:167710/92‑1 (MQ=255)

                                                           gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                              >  1:2192188/1‑93 (MQ=255)

                                                             ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTgg                            >  1:5379818/1‑93 (MQ=255)

                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  1:3935183/93‑1 (MQ=255)

                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  1:2906875/93‑1 (MQ=255)

                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  2:484971/93‑1 (MQ=255)

                                                               aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                          >  2:757145/1‑93 (MQ=255)

                                                               aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                          >  1:4234917/1‑93 (MQ=255)

                                                               aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                          <  1:146540/93‑1 (MQ=255)

                                                               aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGCGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                          >  1:2381623/1‑93 (MQ=255)

                                                                tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTa                         <  2:6379413/93‑1 (MQ=255)

                                                                tCTTTTCGTATAACGATAGAAATGAAGCGTTGATTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTa                         >  1:718654/1‑93 (MQ=255)

                                                                  ttttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                       >  2:5748771/1‑93 (MQ=255)

                                                                      cGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAg                   >  1:2237376/1‑93 (MQ=255)

                                                                      cGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAg                   <  2:5061787/93‑1 (MQ=255)

                                                                       gTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAgg                  >  1:1243876/1‑93 (MQ=255)

                                                                         ataACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAg                >  2:851092/1‑93 (MQ=255)

                                                                          taACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGt               >  2:1897667/1‑93 (MQ=255)

                                                                            aCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTa             >  1:3316222/1‑93 (MQ=255)

                                                                             cGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc            <  2:1972591/93‑1 (MQ=255)

                                                                               aTAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGc          >  1:3861485/1‑93 (MQ=255)

                                                                                 aGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCa        <  2:2078756/93‑1 (MQ=255)

                                                                                  gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg       >  2:3095199/1‑93 (MQ=255)

                                                                                  gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg       <  2:2946618/93‑1 (MQ=255)

                                                                                   aaaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAgg      <  1:5570026/93‑1 (MQ=255)

                                                                                   aaaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAgg      >  2:4639209/1‑93 (MQ=255)

                                                                                       gAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAt  >  2:2474568/1‑93 (MQ=255)

AAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAT  >  NC_000913/4472318‑4472497

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