BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,935	G→A	R106C (CGC→TGC)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,935	0	G	A	98.8%	871.0 / NA	254	R106C (CGC→TGC)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/3); new base A (111/136); total (111/143)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCA  >  NC_000913/308845‑309024

gggcggAAGCAGTAGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCAct                                                                                         >  2:3952946/1‑93 (MQ=255)

 ggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActc                                                                                        <  2:3469287/93‑1 (MQ=255)

   cggAGGCAGTGGCCGCAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtt                                                                                      <  2:2293539/93‑1 (MQ=255)

   cggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtt                                                                                      <  2:2873105/93‑1 (MQ=255)

    ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttt                                                                                     <  2:2940465/93‑1 (MQ=255)

    ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttt                                                                                     <  2:4473255/93‑1 (MQ=255)

    ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttt                                                                                     <  1:1252794/93‑1 (MQ=255)

     gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                    <  2:858526/93‑1 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCCCTCttttt                                                                                   <  2:432101/93‑1 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                 >  2:81564/1‑93 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                 <  1:460202/93‑1 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                 <  2:2626804/93‑1 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                 <  1:5674373/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                               <  1:2489862/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                               <  1:6156681/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                               <  2:4641162/93‑1 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                              >  1:4723314/1‑93 (MQ=255)

             ggCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGc                                                                            >  1:969117/1‑93 (MQ=255)

             ggCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGc                                                                            <  2:5618854/93‑1 (MQ=255)

              gCCACAGCGCTGCGGTTGTCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCg                                                                           <  2:1931013/93‑1 (MQ=255)

              gCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCg                                                                           <  1:1725158/93‑1 (MQ=255)

               ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCgg                                                                          <  1:3719186/93‑1 (MQ=255)

               ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCgg                                                                          <  1:2388779/93‑1 (MQ=255)

                      gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCtt                                                                   >  2:6283665/1‑93 (MQ=255)

                      gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCtt                                                                   >  1:6071814/1‑93 (MQ=255)

                         gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                                >  2:1087111/1‑93 (MQ=255)

                            gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                             <  2:2672355/93‑1 (MQ=255)

                              tGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaaga                                                           >  2:5461322/1‑93 (MQ=255)

                               ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                          <  1:2375693/93‑1 (MQ=255)

                                gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                         >  2:6377620/1‑93 (MQ=255)

                                 cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                        <  1:2875944/93‑1 (MQ=255)

                                 cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                        <  1:6207438/93‑1 (MQ=255)

                                  aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCg                                                       >  1:1623819/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                     >  2:4955178/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                     <  1:3360880/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                     <  1:3238862/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                     >  1:3011981/1‑93 (MQ=255)

                                     aTCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAt                                                    >  2:6203183/1‑93 (MQ=255)

                                      tCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATc                                                   >  2:671002/1‑93 (MQ=255)

                                      tCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTGCATCTGCGGGTCCCTTATAGAAGAAGCGGATc                                                   <  1:3397206/93‑1 (MQ=255)

                                         cgcTGCGCATCACTGAGTGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACt                                                <  2:1523431/93‑1 (MQ=255)

                                         cgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACt                                                <  2:1795058/93‑1 (MQ=255)

                                         cgcTGCGCATCACTGAGGGAGTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACt                                                <  1:5416683/93‑1 (MQ=255)

                                          gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACtt                                               >  2:4738826/1‑93 (MQ=255)

                                          gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACtt                                               <  1:4010488/93‑1 (MQ=255)

                                          gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCCGGACCCTTATAGAAGAAGCGGATCGCtt                                               >  1:1003089/1‑93 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                             >  2:2945508/1‑93 (MQ=255)

                                              cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                           <  1:741020/93‑1 (MQ=255)

                                              cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                           <  1:3925546/93‑1 (MQ=255)

                                               gcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTg                                          <  2:5159159/93‑1 (MQ=255)

                                               gcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTg                                          <  1:4974093/93‑1 (MQ=255)

                                               gcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTg                                          <  2:2368670/93‑1 (MQ=255)

                                                cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTgg                                         <  1:6588099/93‑1 (MQ=255)

                                                cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTgg                                         <  1:3905717/93‑1 (MQ=255)

                                                 aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTggt                                        >  2:1984787/1‑92 (MQ=255)

                                                 aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGc                                        >  2:1313648/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                       <  2:3218481/93‑1 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                       <  2:658064/93‑1 (MQ=255)

                                                   cACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCtt                                      <  1:4064154/93‑1 (MQ=255)

                                                    aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                     >  1:5907774/1‑93 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                    >  2:3605543/1‑93 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                    >  2:6447739/1‑93 (MQ=255)

                                                        aGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCg                                 >  1:1089774/1‑93 (MQ=255)

                                                        aGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCg                                 >  2:420839/1‑93 (MQ=255)

                                                         gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                                >  1:4305664/1‑93 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                               >  2:6651447/1‑93 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                               >  2:944047/1‑93 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                               <  1:892951/93‑1 (MQ=255)

                                                           gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                              >  1:1502465/1‑93 (MQ=255)

                                                             cTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                            >  2:3371170/1‑93 (MQ=255)

                                                               gATTAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          >  2:4131076/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          >  2:522239/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          <  2:5858506/93‑1 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          >  2:5050135/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          >  2:3178852/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                          <  2:256518/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                      <  1:5785705/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                      <  2:376983/93‑1 (MQ=255)

                                                                    aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                     >  2:2132348/1‑93 (MQ=255)

                                                                    aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                     >  1:63925/1‑93 (MQ=255)

                                                                    aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                     >  2:5584474/1‑93 (MQ=255)

                                                                    aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                     >  1:867499/1‑93 (MQ=255)

                                                                     aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                    <  1:2621155/93‑1 (MQ=255)

                                                                     aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGACCCTTATAGAAGAAGCGGATCACTTCGCTGGCGTGGGCGGGTAGCAGCAAGc                    >  2:632031/1‑93 (MQ=255)

                                                                      ccAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCt                   <  1:4193445/93‑1 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGGAGCAGCAAGCTg                  >  2:4938167/1‑93 (MQ=255)

                                                                        aGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTgg                 <  2:2431302/93‑1 (MQ=255)

                                                                         gACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGCAGCAGCAAGCTGGc                >  1:2728745/1‑93 (MQ=255)

                                                                           cAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCgg              <  2:1787612/93‑1 (MQ=255)

                                                                            aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg             <  2:745798/93‑1 (MQ=255)

                                                                               gCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagt          <  2:1572147/93‑1 (MQ=255)

                                                                               gCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagt          <  1:5193797/93‑1 (MQ=255)

                                                                                cGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtg         >  1:5800043/1‑93 (MQ=255)

                                                                                   tagtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGCAGCAGCAAGCTGGCGGgggtgagg      >  1:558622/1‑92 (MQ=255)

                                                                                    agtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTa     <  2:1237292/93‑1 (MQ=255)

                                                                                    agtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTa     <  2:1632094/93‑1 (MQ=255)

                                                                                    agtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTa     <  2:4102347/93‑1 (MQ=255)

                                                                                      tagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGc   <  1:1860297/93‑1 (MQ=255)

                                                                                       agCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCa  <  1:5613898/93‑1 (MQ=255)

                                                                                       agCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCa  >  2:3735049/1‑93 (MQ=255)

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCA  >  NC_000913/308845‑309024

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