BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,098,623	T→C	R174R (CGT→CGC)	ycdX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,098,623	0	T	C	96.6%	360.7 / 0.9	119	R174R (CGT→CGC)	ycdX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (59/56); total (62/57)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 120 total aligned reads displayed.

GCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAAGAG  >  NC_000913/1098534‑1098713

gCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAg                                                                                         <  2:2122867/93‑1 (MQ=255)

 cTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAgg                                                                                        >  2:246387/1‑93 (MQ=255)

  tGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAgga                                                                                       >  1:3735611/1‑92 (MQ=255)

  tGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGt                                                                                       <  1:2241299/93‑1 (MQ=255)

  tGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGt                                                                                       <  2:1593859/93‑1 (MQ=255)

  tGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACGCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGt                                                                                       <  2:1191600/93‑1 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                      >  2:1594719/1‑93 (MQ=255)

     agcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGc                                                                                    <  2:2080623/93‑1 (MQ=255)

     agcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGc                                                                                    <  2:1674769/93‑1 (MQ=255)

      gcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGcc                                                                                   <  2:846037/93‑1 (MQ=255)

      gcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGcc                                                                                   <  2:4703618/93‑1 (MQ=255)

         gcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCa                                                                                >  2:148509/1‑93 (MQ=255)

           cgcgAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGc                                                                              <  2:923967/93‑1 (MQ=255)

             cgAAACATCAGGTGGCGCTGGAAATCACTAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACACCTGTCGCGAGGTAGCCGCAGCgg                                                                            <  1:890133/93‑1 (MQ=255)

              gAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGGCGCAGCGGt                                                                           <  2:4456896/93‑1 (MQ=255)

               aaaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTa                                                                          <  2:2346219/93‑1 (MQ=255)

               aaaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTa                                                                          <  1:2863537/93‑1 (MQ=255)

               aaaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTa                                                                          >  1:3856312/1‑93 (MQ=255)

                aaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTAc                                                                         <  2:3299800/93‑1 (MQ=255)

                aaCATCAGGAGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTAc                                                                         >  2:1556200/1‑93 (MQ=255)

                  cATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGt                                                                       >  2:3673811/1‑93 (MQ=255)

                  cATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGt                                                                       <  1:1768150/93‑1 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                     >  2:1834107/1‑93 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                     >  2:1067626/1‑93 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                     >  1:298525/1‑93 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                     <  1:2763381/93‑1 (MQ=255)

                     cAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGAt                                                                    <  1:3481443/93‑1 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                   >  1:3599804/1‑93 (MQ=255)

                       ggTGGCGCTTGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGc                                                                  <  2:4041844/93‑1 (MQ=255)

                       ggTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGc                                                                  <  2:3047717/93‑1 (MQ=255)

                       ggTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGc                                                                  <  2:470571/93‑1 (MQ=255)

                       ggTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGc                                                                  <  1:4466338/93‑1 (MQ=255)

                        gaggCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGTAGACAACTATCGCTAGTGAGCCGCAGCGGTACGtgctgcg                                                                 >  2:4111378/3‑92 (MQ=255)

                          ggCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTgg                                                               <  2:3168456/93‑1 (MQ=255)

                          ggCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTgg                                                               >  1:3562632/1‑93 (MQ=255)

                           gcgcTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGa                                                              >  2:841626/1‑93 (MQ=255)

                             gcTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAgg                                                            >  2:32606/1‑93 (MQ=255)

                              cTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGt                                                           <  1:4539797/93‑1 (MQ=255)

                               tGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGtt                                                          >  2:4529088/1‑93 (MQ=255)

                                ggAAATCAATAATTTCTCATTTTTGCACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTg                                                         <  2:1934629/93‑1 (MQ=255)

                                ggAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTg                                                         <  2:2877079/93‑1 (MQ=255)

                                ggAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTg                                                         >  1:860853/1‑93 (MQ=255)

                                 gAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTgg                                                        >  1:1768728/1‑93 (MQ=255)

                                  aaaTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTggg                                                       >  2:4543509/1‑93 (MQ=255)

                                   aaTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGt                                                      <  2:586378/93‑1 (MQ=255)

                                   aaTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGt                                                      <  1:1405425/93‑1 (MQ=255)

                                    aTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTggggg                                                     >  2:671477/1‑91 (MQ=255)

                                    aTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTg                                                     >  2:1326851/1‑93 (MQ=255)

                                    aTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGATGGGTg                                                     >  1:1681023/1‑93 (MQ=255)

                                      cAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGc                                                   >  1:944726/1‑93 (MQ=255)

                                        ataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAt                                                 >  2:3199154/1‑93 (MQ=255)

                                        ataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAt                                                 >  2:2664829/1‑93 (MQ=255)

                                        ataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAt                                                 >  2:4501141/1‑93 (MQ=255)

                                         taatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAtt                                                >  2:699025/1‑93 (MQ=255)

                                         taatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAtt                                                >  2:4176411/1‑93 (MQ=255)

                                          aatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTa                                               <  1:2248471/93‑1 (MQ=255)

                                           atTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAg                                              >  2:393013/1‑93 (MQ=255)

                                              ccTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCt                                           >  1:2870484/1‑93 (MQ=255)

                                               cTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTTGCATTAGGCTc                                          <  1:2323991/93‑1 (MQ=255)

                                               cTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTc                                          >  1:994195/1‑93 (MQ=255)

                                               cTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTc                                          >  1:3534050/1‑93 (MQ=255)

                                                tCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCg                                         <  2:3619626/93‑1 (MQ=255)

                                                   tttttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGAt                                      >  1:2981945/1‑93 (MQ=255)

                                                   tttttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGAt                                      <  1:3848619/93‑1 (MQ=255)

                                                   tttttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGAt                                      >  1:1641411/1‑93 (MQ=255)

                                                    ttttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTGCGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGAtt                                     <  2:2136416/93‑1 (MQ=255)

                                                          acTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTcaca                               >  1:3606624/1‑93 (MQ=255)

                                                           cTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCacac                              >  1:2073954/1‑93 (MQ=255)

                                                           cTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCacac                              >  1:2590086/1‑93 (MQ=255)

                                                           cTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCacac                              >  2:2536583/1‑93 (MQ=255)

                                                           cTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCacac                              >  1:2590087/1‑93 (MQ=255)

                                                            tCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACAcc                             >  2:2043268/1‑93 (MQ=255)

                                                            tCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACAcc                             >  1:3092175/1‑93 (MQ=255)

                                                             cACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcg                            >  1:4361066/1‑93 (MQ=255)

                                                             cACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcg                            <  1:3417283/93‑1 (MQ=255)

                                                              aCGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcgc                           >  2:572655/1‑93 (MQ=255)

                                                              aCGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcgc                           <  2:4659026/93‑1 (MQ=255)

                                                              aCGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcgc                           >  1:4504372/1‑93 (MQ=255)

                                                                gTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGt                         >  2:3893221/1‑93 (MQ=255)

                                                                 tAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGtt                        >  2:4592187/1‑93 (MQ=255)

                                                                  aaGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGttt                       >  1:1101041/1‑93 (MQ=255)

                                                                    gggCAGGGAAGACAACTGGCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTAc                     <  1:1757787/93‑1 (MQ=255)

                                                                      gCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCa                   <  1:4440263/93‑1 (MQ=255)

                                                                      gCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCa                   <  1:2401886/93‑1 (MQ=255)

                                                                       cAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCAt                  <  1:1341945/93‑1 (MQ=255)

                                                                       cAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGCTTACCAt                  >  2:1500136/1‑93 (MQ=255)

                                                                        aGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTATCACACCGCGTTTACCATg                 <  1:339956/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg                <  1:3091746/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg                >  1:2166218/1‑93 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg                <  1:541328/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg                >  1:4487720/1‑93 (MQ=255)

                                                                          tGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATggg               <  1:1246599/93‑1 (MQ=255)

                                                                            aaGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATggggg             <  1:4246563/93‑1 (MQ=255)

                                                                            aaGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATggggg             <  1:3569443/93‑1 (MQ=255)

                                                                             aGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGa            <  2:4253549/93‑1 (MQ=255)

                                                                             aGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGa            >  1:2368109/1‑93 (MQ=255)

                                                                             aGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGa            >  1:2590509/1‑93 (MQ=255)

                                                                                  aCTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTg       >  2:4130280/1‑93 (MQ=255)

                                                                                   cTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGa      >  2:3992961/1‑93 (MQ=255)

                                                                                     gTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAag    >  1:1150110/1‑93 (MQ=255)

GCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAAGAG  >  NC_000913/1098534‑1098713

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