BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	628,816	A→G	V98V (GTA→GTG)	ybdB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,816	0	A	G	100.0%	533.0 / NA	169	V98V (GTA→GTG)	ybdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (89/80); total (89/80)

Only 100 of 180 total aligned reads displayed.

GCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628727‑628905

gCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTgcgc                                                                                        <  2:4663488/93‑1 (MQ=255)

gCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTgcgc                                                                                        >  2:1096344/1‑93 (MQ=255)

gCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTgcgc                                                                                        <  2:2471185/93‑1 (MQ=255)

   ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     >  2:2625392/1‑93 (MQ=255)

    gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    >  1:792752/1‑93 (MQ=255)

     aTTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGGGGGAAAGGTGCGCGGCGt                                                                                   >  1:2512196/1‑93 (MQ=255)

      tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  >  2:2169034/1‑93 (MQ=255)

      tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  >  2:397199/1‑93 (MQ=255)

      tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  <  2:3549453/93‑1 (MQ=255)

      tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  <  1:3687358/93‑1 (MQ=255)

       ttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCt                                                                                 >  1:495289/1‑93 (MQ=255)

       tgatgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCt                                                                                 <  1:3131828/91‑1 (MQ=255)

        tATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTg                                                                                >  2:2323666/1‑93 (MQ=255)

         atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               <  2:2881855/93‑1 (MQ=255)

         atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               >  1:4247579/1‑93 (MQ=255)

          tgacgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGTTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              <  2:2335371/93‑1 (MQ=255)

           gatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCa                                                                             >  1:1615281/1‑93 (MQ=255)

           gatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCa                                                                             <  2:3712076/93‑1 (MQ=255)

                cccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGc                                                                        <  2:3556335/93‑1 (MQ=255)

                 ccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCt                                                                       >  1:3440569/1‑93 (MQ=255)

                 ccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCt                                                                       <  2:3752981/93‑1 (MQ=255)

                  cgcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTg                                                                      >  1:1215450/1‑93 (MQ=255)

                    cgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCa                                                                    >  2:4142375/1‑93 (MQ=255)

                        ggACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCtt                                                                <  1:2082686/93‑1 (MQ=255)

                        ggACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCtt                                                                >  1:2106656/1‑93 (MQ=255)

                         gACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTg                                                               >  2:3360341/1‑93 (MQ=255)

                          aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              <  2:4413106/93‑1 (MQ=255)

                            aGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTc                                                            <  2:3529476/93‑1 (MQ=255)

                            aGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTc                                                            >  2:2837939/1‑93 (MQ=255)

                                 gtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaa                                                       >  2:824448/1‑93 (MQ=255)

                                  tgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaaa                                                      >  2:3570560/1‑93 (MQ=255)

                                   ggTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaaaa                                                     <  2:2209666/93‑1 (MQ=255)

                                    gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  1:193928/1‑93 (MQ=255)

                                    gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  1:4541823/1‑93 (MQ=255)

                                     tAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATc                                                   <  2:2624373/93‑1 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  2:3504863/93‑1 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                >  1:72898/1‑93 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  1:511605/93‑1 (MQ=255)

                                         cacaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               <  2:1539259/93‑1 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              <  2:1505644/93‑1 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              <  1:3061790/93‑1 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  1:592806/1‑93 (MQ=255)

                                            aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            >  1:4304320/1‑93 (MQ=255)

                                             gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           <  1:2383350/93‑1 (MQ=255)

                                               aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         >  1:462693/1‑93 (MQ=255)

                                                cTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCAGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaa                                        >  1:292653/1‑93 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  1:1681828/1‑93 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       <  2:2718381/93‑1 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       <  2:4005121/93‑1 (MQ=255)

                                                  tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTGTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  2:4137010/93‑1 (MQ=255)

                                                  tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  2:4082529/93‑1 (MQ=255)

                                                  tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  2:902123/93‑1 (MQ=255)

                                                  tAAGGCAACACCCCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  2:2693923/93‑1 (MQ=255)

                                                   aaTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtc                                     >  2:4552954/1‑93 (MQ=255)

                                                     tGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgt                                   <  2:3923696/93‑1 (MQ=255)

                                                      gCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  >  2:1106466/1‑93 (MQ=255)

                                                      gCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  >  2:3140821/1‑93 (MQ=255)

                                                      gCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  >  1:3937873/1‑93 (MQ=255)

                                                       cAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtcg                                 <  2:4359614/93‑1 (MQ=255)

                                                        aaCACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                >  1:3383474/1‑93 (MQ=255)

                                                         acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               <  1:3164903/93‑1 (MQ=255)

                                                         acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               >  1:3933400/1‑93 (MQ=255)

                                                          tgcacCATCGCTCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  1:1523808/91‑1 (MQ=255)

                                                          cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  1:1766635/93‑1 (MQ=255)

                                                           acacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtttt                             <  1:3416220/93‑1 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGGAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  2:2788145/1‑93 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           <  1:3797138/93‑1 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  1:1361274/1‑93 (MQ=255)

                                                              ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  1:2166472/93‑1 (MQ=255)

                                                              ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:4106612/93‑1 (MQ=255)

                                                              ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:3154134/93‑1 (MQ=255)

                                                               cATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         >  2:2999058/1‑93 (MQ=255)

                                                                aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        <  2:3220608/93‑1 (MQ=255)

                                                                aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        >  2:43672/1‑93 (MQ=255)

                                                                 tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       <  1:2714494/93‑1 (MQ=255)

                                                                 tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       >  1:102408/1‑93 (MQ=255)

                                                                 tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       <  1:3993696/93‑1 (MQ=255)

                                                                  cGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGaa                      <  2:1607854/93‑1 (MQ=255)

                                                                  cGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGaa                      <  2:1247947/93‑1 (MQ=255)

                                                                  cGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGaa                      >  1:4079870/1‑93 (MQ=255)

                                                                    cccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGcaca                    >  1:3431346/1‑93 (MQ=255)

                                                                     ccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAg                   <  1:1062837/93‑1 (MQ=255)

                                                                      cGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  >  2:3617626/1‑93 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGTCGTCTGCCAGCCGCTGCATCTTGGTCGTCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 <  2:5445/93‑1 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAtgg                 >  1:4480292/1‑93 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 <  1:244924/93‑1 (MQ=255)

                                                                        gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:333149/93‑1 (MQ=255)

                                                                            cTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcggc            >  1:4216849/1‑93 (MQ=255)

                                                                             tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCTTCTTTGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTTGGTGAACCGGGgaggcg           >  2:3442024/1‑93 (MQ=255)

                                                                             tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  1:2242979/1‑93 (MQ=255)

                                                                             tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  2:651592/1‑93 (MQ=255)

                                                                              gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  2:1268883/1‑93 (MQ=255)

                                                                              gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  2:3668319/1‑93 (MQ=255)

                                                                                gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAAAGGGGGGGGCGTtg        >  2:193702/1‑93 (MQ=255)

                                                                                 ggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgc       >  2:4184433/1‑93 (MQ=255)

                                                                                 ggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgc       <  1:3445891/93‑1 (MQ=255)

                                                                                   aaaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctg     >  2:2059949/1‑93 (MQ=255)

                                                                                   aaaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGGtgctg     >  2:1770349/1‑93 (MQ=255)

                                                                                    aaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctgc    <  1:133351/93‑1 (MQ=255)

                                                                                      ggTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAc  <  2:790961/93‑1 (MQ=255)

GCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628727‑628905

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