BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,588	C→G	R86T (AGG→ACG)	gppA ←	guanosine pentaphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,588	0	C	G	100.0%	1013.5 / NA	284	R86T (AGG→ACG)	gppA	guanosine pentaphosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (155/129); total (155/129)

Only 100 of 290 total aligned reads displayed.

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  2:1521985/1‑93 (MQ=255)

 ccTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                                        >  1:858705/1‑93 (MQ=255)

  ctctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGc                                                                                       >  1:1501482/1‑93 (MQ=255)

   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      >  2:5858231/1‑93 (MQ=255)

    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     >  2:1100904/1‑93 (MQ=255)

    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     >  1:4356923/1‑93 (MQ=255)

    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     >  2:3521962/1‑93 (MQ=255)

     ttCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                                    <  2:4345841/93‑1 (MQ=255)

      tCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                                   <  2:5282714/93‑1 (MQ=255)

        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 <  2:3781460/93‑1 (MQ=255)

        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 <  2:1967303/93‑1 (MQ=255)

          cGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGa                                                                               >  1:4677718/1‑93 (MQ=255)

           gCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGAc                                                                              >  1:4655275/1‑93 (MQ=255)

           gCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGAc                                                                              <  2:4795644/93‑1 (MQ=255)

            cTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                                             <  1:4196238/93‑1 (MQ=255)

            cTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                                             <  1:3395339/93‑1 (MQ=255)

             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  2:969407/1‑93 (MQ=255)

             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  1:718554/1‑93 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:1743047/1‑93 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  1:1781761/1‑93 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:4426810/1‑93 (MQ=255)

                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         >  2:1463360/1‑93 (MQ=255)

                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         <  1:3468419/93‑1 (MQ=255)

                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         <  1:589200/93‑1 (MQ=255)

                 cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                                        >  1:3911542/1‑93 (MQ=255)

                 cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                                        >  2:4026004/1‑93 (MQ=255)

                  aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                                       >  2:1287474/1‑93 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      >  1:2876376/1‑93 (MQ=255)

                    cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                                     <  1:1163494/93‑1 (MQ=255)

                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    >  1:994666/1‑93 (MQ=255)

                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    >  1:1336842/1‑93 (MQ=255)

                       tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                                  <  2:2336086/93‑1 (MQ=255)

                       gaCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                                  <  1:2733994/92‑1 (MQ=255)

                       aaCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                                  <  1:415122/92‑1 (MQ=255)

                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                                                 >  1:1565226/1‑93 (MQ=255)

                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                                                 >  1:1458985/1‑93 (MQ=255)

                          cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                                               <  1:3276318/93‑1 (MQ=255)

                          cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                                               <  1:4768622/93‑1 (MQ=255)

                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgag                                                              >  2:5132236/1‑93 (MQ=255)

                            gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagg                                                             <  1:439788/93‑1 (MQ=255)

                             aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgaggg                                                            <  1:2159983/93‑1 (MQ=255)

                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGAGACAACGCGAATTTGGgaggggg                                                          >  1:1537806/1‑91 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  2:4755815/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  2:460718/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  1:3638768/1‑93 (MQ=255)

                                 ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                                        <  1:5290903/93‑1 (MQ=255)

                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                                      >  2:4520662/1‑93 (MQ=255)

                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                                      >  2:951486/1‑93 (MQ=255)

                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                                   >  2:1353858/1‑93 (MQ=255)

                                      aaTTTCCTGCGCTTTGCCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                                   <  2:492951/93‑1 (MQ=255)

                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                                  <  1:5791755/93‑1 (MQ=255)

                                       aTTTCATGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                                  <  1:5515707/93‑1 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                                 <  1:4914283/93‑1 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                                 >  2:5344658/1‑93 (MQ=255)

                                          tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                                               >  1:894526/1‑93 (MQ=255)

                                          tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                                               >  1:4026477/1‑93 (MQ=255)

                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                                              >  2:3038201/1‑93 (MQ=255)

                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                                              <  1:3352482/93‑1 (MQ=255)

                                            cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                                             >  1:3249537/1‑93 (MQ=255)

                                            cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                                             <  2:643259/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  2:6049622/1‑93 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  1:4051162/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  1:5409210/1‑93 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  1:2354758/1‑93 (MQ=255)

                                              gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                                           <  2:4973913/93‑1 (MQ=255)

                                               cgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACg                                          <  1:4327603/93‑1 (MQ=255)

                                               cgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACg                                          <  1:583671/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTACCGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  2:2893896/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  1:5008730/93‑1 (MQ=255)

                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       >  2:2192521/1‑93 (MQ=255)

                                                       cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                                  >  1:5647659/1‑93 (MQ=255)

                                                        aaTAAAATCACCCGCATTGACGGCAAGACGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                 <  1:5733339/93‑1 (MQ=255)

                                                         aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaaa                                >  2:3599728/1‑92 (MQ=255)

                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               >  2:4206166/1‑93 (MQ=255)

                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               <  2:3495554/93‑1 (MQ=255)

                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               >  1:2611489/1‑93 (MQ=255)

                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  2:2123036/1‑93 (MQ=255)

                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     >  2:3306475/1‑93 (MQ=255)

                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     <  1:933310/93‑1 (MQ=255)

                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     <  2:5815312/93‑1 (MQ=255)

                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     <  1:2929544/93‑1 (MQ=255)

                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATATTGCAGACGTTCAGCAAACAGGCGCAAACa                     >  1:6009758/1‑93 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  1:1006280/93‑1 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  2:2483163/93‑1 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  2:2403589/93‑1 (MQ=255)

                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   <  2:5074954/93‑1 (MQ=255)

                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  1:2093680/1‑93 (MQ=255)

                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATAGCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  1:5152778/1‑93 (MQ=255)

                                                                        ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                 >  1:4215013/1‑93 (MQ=255)

                                                                         tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                >  1:2557273/1‑93 (MQ=255)

                                                                           aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCaa              <  1:1249775/93‑1 (MQ=255)

                                                                           aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCaa              <  2:791502/93‑1 (MQ=255)

                                                                             ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc            >  2:571430/1‑93 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           <  2:5194458/93‑1 (MQ=255)

                                                                                aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  1:2081037/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:6120487/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:1043310/1‑93 (MQ=255)

                                                                                    cGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGGGGAATATCTTGCAGACGTTCAGCAAACAGGCGAAAACATTGCAAACCGCGCTcc     >  2:1870127/1‑93 (MQ=255)

                                                                                    cGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTcc     >  2:2698299/1‑93 (MQ=255)

                                                                                     gTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    >  1:3037749/1‑93 (MQ=255)

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

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