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	99.1%	1137.4 / NA	335	R106C (CGC→TGC)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (159/170); minor base C (1/2); total (160/175)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

GGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCA  >  NC_000913/308846‑309024

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

 gcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActct                                                                                       <  1:1084736/93‑1 (MQ=255)

 gcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActct                                                                                       >  1:5231999/1‑93 (MQ=255)

  cggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtt                                                                                      >  1:5146417/1‑93 (MQ=255)

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

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

    gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                    >  2:3026782/1‑93 (MQ=255)

    gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                    >  2:3086961/1‑93 (MQ=255)

    gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                    >  2:5870833/1‑93 (MQ=255)

    gAAGCAGTGGCCACAGCGCTGCGGATGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGGAGTAGCACTCtttt                                                                                    >  2:894625/1‑93 (MQ=255)

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

        cAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCAt                                                                                >  2:5906469/1‑93 (MQ=255)

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

           tGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTg                                                                             <  1:1693441/93‑1 (MQ=255)

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

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

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

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

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

               cacaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCggg                                                                         <  2:5547175/93‑1 (MQ=255)

               cacaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCggg                                                                         >  1:4950505/1‑93 (MQ=255)

                ccaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGt                                                                        <  1:6293143/92‑1 (MQ=255)

                  aGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTcc                                                                      >  1:5605537/1‑93 (MQ=255)

                   gcgcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTccc                                                                     <  1:4162568/93‑1 (MQ=255)

                      cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                  >  1:6593190/1‑93 (MQ=255)

                      cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                  <  2:7316424/93‑1 (MQ=255)

                      cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                  <  2:5671036/93‑1 (MQ=255)

                       tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                 <  2:2612110/93‑1 (MQ=255)

                       tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                 <  2:3557837/93‑1 (MQ=255)

                        gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTtata                                                                <  2:4637266/93‑1 (MQ=255)

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

                        gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                                <  2:1777063/93‑1 (MQ=255)

                        gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                                <  1:325291/93‑1 (MQ=255)

                          ggTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATaga                                                              <  1:4737154/93‑1 (MQ=255)

                           gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                             >  1:543244/1‑93 (MQ=255)

                            ttGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaag                                                            <  2:3612789/93‑1 (MQ=255)

                             tGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaaga                                                           <  2:2909609/93‑1 (MQ=255)

                              ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                          >  1:68312/1‑93 (MQ=255)

                               gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                         <  1:7145335/93‑1 (MQ=255)

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

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

                                  ttATCGCGCTGCGCATCTCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCgg                                                      >  2:7040781/1‑93 (MQ=255)

                                   tATCGCCCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                     <  2:5608832/93‑1 (MQ=255)

                                       gcgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                 <  1:1119693/93‑1 (MQ=255)

                                       gcgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                 <  1:4587468/93‑1 (MQ=255)

                                       gcgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                 >  1:2604219/1‑93 (MQ=255)

                                         gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATATAAGAAGCGGATCACtt                                               <  1:2425027/93‑1 (MQ=255)

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

                                          cTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTc                                              <  2:1673736/93‑1 (MQ=255)

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

                                            gcgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGc                                            <  1:5865473/93‑1 (MQ=255)

                                              gcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTg                                          >  1:5524829/1‑93 (MQ=255)

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

                                                 tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                       >  1:7732325/1‑93 (MQ=255)

                                                 tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                       >  1:1243326/1‑93 (MQ=255)

                                                 tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                       >  1:3050631/1‑93 (MQ=255)

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

                                                   aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                     <  2:5431043/93‑1 (MQ=255)

                                                    cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                    <  1:2910350/93‑1 (MQ=255)

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

                                                    cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                    <  2:4321040/93‑1 (MQ=255)

                                                     tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                   <  1:5489066/93‑1 (MQ=255)

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

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

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

                                                           ccTAATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                             <  1:6398900/93‑1 (MQ=255)

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

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

                                                             tGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagc                           >  1:1949720/1‑93 (MQ=255)

                                                             tGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTCGCTTGGGCGGGTagc                           >  2:2105511/1‑93 (MQ=255)

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

                                                               aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                         <  1:1539210/93‑1 (MQ=255)

                                                               aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                         >  2:2771328/1‑93 (MQ=255)

                                                                tCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcagc                        >  1:2751311/1‑93 (MQ=255)

                                                                  aaaGCAGACGATGCGGTGGTAGCGCTTTTTGTGATCGGGGGCTCCCTGATAGAAGAAGCGGATCAGTTCGCTGGCTTGGGCGGGTAGCAGCaa                      <  1:6575554/93‑1 (MQ=255)

                                                                  aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTCTCGAAGAAGCGGGTTCCCTCGCTGGCTTGGGCGGTGAacagcac                      >  1:6099913/1‑92 (MQ=255)

                                                                  aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                      >  2:2368917/1‑93 (MQ=255)

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

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

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

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

                                                                    aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                    >  2:103490/1‑93 (MQ=255)

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

                                                                     ccAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCt                   >  1:4871138/1‑93 (MQ=255)

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

                                                                       aGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTgg                 >  1:7804401/1‑93 (MQ=255)

                                                                        gACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGc                >  2:4137384/1‑93 (MQ=255)

                                                                        gACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGc                <  2:6442133/93‑1 (MQ=255)

                                                                           aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg             >  2:516744/1‑93 (MQ=255)

                                                                             tGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgag           >  2:4066429/1‑93 (MQ=255)

                                                                                ggTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga        <  1:4270420/93‑1 (MQ=255)

                                                                                ggTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga        >  2:4952795/1‑93 (MQ=255)

                                                                                ggTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga        >  1:1752320/1‑93 (MQ=255)

                                                                                 gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag       <  1:3555540/93‑1 (MQ=255)

                                                                                 gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag       <  2:4727470/93‑1 (MQ=255)

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

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

                                                                                    gtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAg    >  1:4012763/1‑93 (MQ=255)

                                                                                     tagCCCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGc   <  2:4542975/93‑1 (MQ=255)

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

GGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCA  >  NC_000913/308846‑309024

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