BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,298,112	T→C	G96G (GGA→GGG)	yhbS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,298,112	0	T	C	99.2%	817.5 / ‑4.1	239	G96G (GGA→GGG)	yhbS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (117/119); total (119/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 248 total aligned reads displayed.

GCCGCCAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGTCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCTTC  >  NC_000913/3298023‑3298201

 ccgccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCtt                                                                                       <  1:525754/93‑1 (MQ=255)

  cgccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTc                                                                                      >  1:2096652/1‑93 (MQ=255)

  cgccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTc                                                                                      <  1:16191/93‑1 (MQ=255)

   gccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCa                                                                                     >  2:12697/1‑93 (MQ=255)

   gccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCa                                                                                     >  2:3821412/1‑93 (MQ=255)

    ccAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCAt                                                                                    <  1:4098766/93‑1 (MQ=255)

     cAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATa                                                                                   >  2:45616/1‑93 (MQ=255)

     cAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCCAGCCCTTCATa                                                                                   >  1:312736/1‑93 (MQ=255)

       gTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATCGa                                                                                 >  1:1479722/1‑93 (MQ=255)

       gTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGa                                                                                 <  1:648010/93‑1 (MQ=255)

        ttCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGAc                                                                                <  1:1257362/93‑1 (MQ=255)

          cAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCa                                                                              <  2:2665827/93‑1 (MQ=255)

            aaGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGt                                                                            <  2:2870635/93‑1 (MQ=255)

            aaGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGt                                                                            <  2:2615470/93‑1 (MQ=255)

             aGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGtt                                                                           <  1:5927862/93‑1 (MQ=255)

             aGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGtt                                                                           >  2:4475647/1‑93 (MQ=255)

              gCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTg                                                                          <  1:4152245/93‑1 (MQ=255)

              gCCGAAACGACTGTACAGCACCGGAACGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAAACGAGCCCTTCATAGACCAGTTg                                                                          <  1:5649703/93‑1 (MQ=255)

               ccGAAACGACTTTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTgg                                                                         <  1:4980219/93‑1 (MQ=255)

                cGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGCACAGTTGGc                                                                        >  2:2154856/1‑93 (MQ=255)

                cGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGc                                                                        >  1:5846707/1‑93 (MQ=255)

                 gAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCg                                                                       >  1:4261086/1‑93 (MQ=255)

                    aCGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGc                                                                    >  1:5189348/1‑93 (MQ=255)

                      gACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCa                                                                  <  1:2821157/93‑1 (MQ=255)

                      gACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCa                                                                  <  2:64978/93‑1 (MQ=255)

                      gACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCa                                                                  >  2:2813122/1‑93 (MQ=255)

                           tACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGcccc                                                             <  1:917139/93‑1 (MQ=255)

                            aCAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCt                                                            >  2:5180466/1‑93 (MQ=255)

                            aCAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCt                                                            >  2:2935944/1‑93 (MQ=255)

                              aGCGCCGGATCG‑CCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGGCCTTCATAAAACAGTTGGGGTGCCAAGCCCCGTc                                                         >  2:4461112/1‑93 (MQ=255)

                                cgcCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTcc                                                        <  2:5863563/93‑1 (MQ=255)

                                cgcCGGATCGCCCACCGTCA‑CACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCa                                                       <  1:1900165/93‑1 (MQ=255)

                                 gcCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCa                                                       >  2:777217/1‑93 (MQ=255)

                                  ccGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCAc                                                      <  1:713196/93‑1 (MQ=255)

                                  ccGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCAc                                                      >  2:4059887/1‑93 (MQ=255)

                                  ccGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCAc                                                      >  1:3898712/1‑93 (MQ=255)

                                    ggATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACgg                                                    <  2:212461/93‑1 (MQ=255)

                                    ggATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACgg                                                    >  1:3082037/1‑93 (MQ=255)

                                    ggATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGACCCTGTCCACgg                                                    >  1:2522718/1‑93 (MQ=255)

                                      aTCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTa                                                  <  2:3637521/93‑1 (MQ=255)

                                      aTCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTa                                                  <  1:4004824/93‑1 (MQ=255)

                                        cGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGCCCACGGCAtt                                                >  2:5725436/1‑93 (MQ=255)

                                         gCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTAttt                                               <  2:5994963/93‑1 (MQ=255)

                                           ccAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTAttttt                                             <  1:542919/93‑1 (MQ=255)

                                            cAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtc                                            >  1:289221/1‑93 (MQ=255)

                                            cAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtc                                            >  2:4643163/1‑93 (MQ=255)

                                             aGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtca                                           >  2:5434737/1‑93 (MQ=255)

                                             aGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtca                                           >  2:5711857/1‑93 (MQ=255)

                                              gCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtcat                                          >  2:3351367/1‑93 (MQ=255)

                                              gCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtcat                                          >  2:5101055/1‑93 (MQ=255)

                                               cGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTtcatc                                         <  1:2634672/93‑1 (MQ=255)

                                                gTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCTCCTGTCCACGGTATGTTtcatca                                        >  1:4278880/1‑93 (MQ=255)

                                                 tCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCaa                                       >  2:5240399/1‑93 (MQ=255)

                                                 tCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCaa                                       <  2:5286058/93‑1 (MQ=255)

                                                   accacTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAAcc                                     <  1:4892164/93‑1 (MQ=255)

                                                   accacTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAAcc                                     >  2:3478696/1‑93 (MQ=255)

                                                   accacTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAAcc                                     <  1:1139922/93‑1 (MQ=255)

                                                     cacTGCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGc                                   >  1:1711465/1‑93 (MQ=255)

                                                       cTGCGGAATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGCCTTTTTCATCAACCGCCa                                 >  1:2562156/1‑93 (MQ=255)

                                                         gCGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAAc                               >  2:1689050/1‑93 (MQ=255)

                                                          cGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACg                              <  1:3925708/93‑1 (MQ=255)

                                                          cGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACg                              >  2:3126740/1‑93 (MQ=255)

                                                          cGGCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGCATTTTTCATCAAACGCCAACg                              >  1:4522558/1‑93 (MQ=255)

                                                           ggCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACgg                             <  2:1989195/93‑1 (MQ=255)

                                                            gCATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGt                            >  1:878000/1‑93 (MQ=255)

                                                             cATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGtc                           >  2:4930808/1‑92 (MQ=255)

                                                             cATAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTg                           <  1:1213092/93‑1 (MQ=255)

                                                              aTAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGc                          >  2:3214636/1‑93 (MQ=255)

                                                               tAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGcc                         <  1:3713470/93‑1 (MQ=255)

                                                               tAGCCGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGcc                         <  1:54542/93‑1 (MQ=255)

                                                                   cGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGc                     <  2:826665/93‑1 (MQ=255)

                                                                   cGAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGc                     <  2:2379901/93‑1 (MQ=255)

                                                                    gAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGcc                    >  2:3541395/1‑93 (MQ=255)

                                                                    gAACTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCAAGGTATTTTTCATCAACCGCCAACGGTGCCATGcc                    >  2:3996355/1‑93 (MQ=255)

                                                                     aaCTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACTGCCAACGGTGCCATGCCg                   >  1:1957557/1‑93 (MQ=255)

                                                                     aaCTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCg                   >  2:4805487/1‑93 (MQ=255)

                                                                      aCTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGa                  <  1:3261507/93‑1 (MQ=255)

                                                                      aCTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGGGCCATGCCgc                  >  1:209892/1‑92 (MQ=255)

                                                                       cTCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGAc                 >  1:4718646/1‑93 (MQ=255)

                                                                        tCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGAcc                >  1:2148863/1‑93 (MQ=255)

                                                                        tCATTAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGAcc                >  1:2529323/1‑93 (MQ=255)

                                                                           ttAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCAt             <  2:5360244/93‑1 (MQ=255)

                                                                            tAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCAtt            >  1:1387122/1‑93 (MQ=255)

                                                                            tAAGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGGATTTTTCATCAAACGCCAACGGTGCCATGCAGACCCAtt            >  1:1574080/1‑93 (MQ=255)

                                                                             aaGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTg           >  2:2452165/1‑93 (MQ=255)

                                                                             aaGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTg           >  1:977919/1‑93 (MQ=255)

                                                                              aGCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGGCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGc          <  2:4176806/93‑1 (MQ=255)

                                                                               gCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCa         >  1:3983986/1‑93 (MQ=255)

                                                                               gCGAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCa         <  1:5728326/93‑1 (MQ=255)

                                                                                 gAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAgg       <  1:3397675/93‑1 (MQ=255)

                                                                                 gAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAgg       <  1:4796042/93‑1 (MQ=255)

                                                                                 gAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAgg       <  2:1143862/93‑1 (MQ=255)

                                                                                 gAATCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAgg       <  2:1905809/93‑1 (MQ=255)

                                                                                   aTCGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTc     <  2:1495990/93‑1 (MQ=255)

                                                                                   aTCAAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTc     <  1:1817009/93‑1 (MQ=255)

                                                                                     cGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCtt   <  1:2176825/93‑1 (MQ=255)

                                                                                     cGAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCtt   <  2:2699922/93‑1 (MQ=255)

                                                                                      gAGGCCTTCATAGACCAGTTGGCGTGCCATCCCCTGTCCATGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCTTc  <  2:6064779/93‑1 (MQ=255)

                                                                                      gAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCTTc  <  1:4466713/93‑1 (MQ=255)

                                                                                      gAGCCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCTTc  <  2:4226489/93‑1 (MQ=255)

GCCGCCAGTTCAAAGCCGAAACGACTGTACAGCGCCGGATCGCCCAGCGTCACCACTGCGGCATAGCCGAACTCATTAAGCGAATCGAGTCCTTCATAGACCAGTTGGCGTGCCAGCCCCTGTCCACGGTATTTTTCATCAACCGCCAACGGTGCCATGCCGACCCATTGCAGGTCTTC  >  NC_000913/3298023‑3298201

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