BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,977,174	G→A	K44K (AAG→AAA)	lysR →	LysR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,977,174	0	G	A	99.6%	1071.3 / NA	277	K44K (AAG→AAA)	lysR	LysR
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (140/135); total (140/137)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GGTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATC  >  NC_000913/2977084‑2977263

ggTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGt                                                                                         >  1:5639012/1‑93 (MQ=255)

 gTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTg                                                                                        >  2:4531557/1‑93 (MQ=255)

 gTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTg                                                                                        <  1:2871232/93‑1 (MQ=255)

 gTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTg                                                                                        <  2:5692243/93‑1 (MQ=255)

   aaTGACCGCCGGAAGCCTGACTGACGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGAt                                                                                      <  1:5777445/93‑1 (MQ=255)

      gACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCgg                                                                                   <  2:1172517/93‑1 (MQ=255)

      gACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCgg                                                                                   >  2:1945133/1‑93 (MQ=255)

      gACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCgg                                                                                   >  1:3141999/1‑93 (MQ=255)

       aCCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCggg                                                                                  <  2:3126510/93‑1 (MQ=255)

       aCCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCggg                                                                                  >  2:1797169/1‑93 (MQ=255)

       aCCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCggg                                                                                  >  1:2400060/1‑93 (MQ=255)

        ccgccgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGc                                                                                 <  1:1638742/93‑1 (MQ=255)

         cgccgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCt                                                                                <  2:581452/93‑1 (MQ=255)

         cgccgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCt                                                                                <  1:613669/93‑1 (MQ=255)

          gccgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTg                                                                               >  1:2172451/1‑93 (MQ=255)

          gccgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTg                                                                               >  2:209595/1‑93 (MQ=255)

            cgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGaa                                                                             >  2:6066834/1‑93 (MQ=255)

            cgGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGaa                                                                             >  2:5443542/1‑93 (MQ=255)

             ggAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGaaa                                                                            >  2:2662487/1‑93 (MQ=255)

             ggAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGaaa                                                                            >  2:535054/1‑93 (MQ=255)

              gAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAAt                                                                           <  1:5025435/93‑1 (MQ=255)

              gAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAAt                                                                           >  2:4297127/1‑93 (MQ=255)

              gAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAAt                                                                           >  1:1595841/1‑93 (MQ=255)

              gAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAAt                                                                           <  2:376282/93‑1 (MQ=255)

               aaGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAAtt                                                                          >  2:3511239/1‑93 (MQ=255)

                aGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTg                                                                         >  1:4141042/1‑93 (MQ=255)

                aGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGATGAAATTg                                                                         >  1:5804429/1‑93 (MQ=255)

                  ccTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGtt                                                                       <  1:3503828/93‑1 (MQ=255)

                  ccTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGtt                                                                       <  2:3412562/93‑1 (MQ=255)

                  ccTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGtt                                                                       >  1:5146630/1‑93 (MQ=255)

                  ccTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGtt                                                                       <  1:1760526/93‑1 (MQ=255)

                   ctgactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGttt                                                                      <  1:4686068/93‑1 (MQ=255)

                    tgactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTg                                                                     >  2:3262161/1‑93 (MQ=255)

                    tgactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTg                                                                     >  1:3057928/1‑93 (MQ=255)

                     gactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGa                                                                    <  1:5361949/93‑1 (MQ=255)

                     gactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGa                                                                    <  1:6000256/93‑1 (MQ=255)

                     gactgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGa                                                                    >  2:1064052/1‑93 (MQ=255)

                      actgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAg                                                                   >  1:327553/1‑93 (MQ=255)

                      actgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAg                                                                   >  1:6421947/1‑93 (MQ=255)

                       ctgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAgc                                                                  >  1:762634/1‑93 (MQ=255)

                       ctgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAgc                                                                  <  1:4579250/93‑1 (MQ=255)

                       ctgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAgc                                                                  <  1:3557963/93‑1 (MQ=255)

                        tgaGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAgcg                                                                 >  1:2365966/1‑93 (MQ=255)

                           ggcggcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGt                                                              <  2:70182/93‑1 (MQ=255)

                            gcggcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTa                                                             <  1:4192091/93‑1 (MQ=255)

                            gcggcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTa                                                             >  2:477200/1‑93 (MQ=255)

                              ggcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACg                                                           >  2:5170522/1‑93 (MQ=255)

                               gcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGt                                                          <  1:4602251/93‑1 (MQ=255)

                               gcACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACCTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGt                                                          <  1:6378689/93‑1 (MQ=255)

                                cacaCCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTg                                                         <  1:4432175/93‑1 (MQ=255)

                                 acacCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTgg                                                        >  2:5950431/1‑93 (MQ=255)

                                  cacCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTggg                                                       >  1:5433545/1‑93 (MQ=255)

                                  cacCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTggg                                                       >  1:2750602/1‑93 (MQ=255)

                                     cTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGa                                                    <  2:4921774/93‑1 (MQ=255)

                                       gCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGAtt                                                  >  1:3089844/1‑93 (MQ=255)

                                        cTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTa                                                 >  1:3734323/1‑93 (MQ=255)

                                        cTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTa                                                 <  1:3873458/93‑1 (MQ=255)

                                          acacacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACa                                               >  2:4812871/1‑93 (MQ=255)

                                          acacacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACa                                               >  1:3936348/1‑93 (MQ=255)

                                          acacacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACa                                               >  1:5198814/1‑93 (MQ=255)

                                             cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATcc                                            >  1:1083013/1‑93 (MQ=255)

                                              cccTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCt                                           <  1:4774873/92‑1 (MQ=255)

                                                  cacaGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCg                                       <  2:1813270/93‑1 (MQ=255)

                                                   acaGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGt                                      >  2:4017765/1‑93 (MQ=255)

                                                   acaGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGt                                      <  2:2457303/93‑1 (MQ=255)

                                                    caGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTg                                     >  1:3086571/1‑93 (MQ=255)

                                                        cAACCGTCAGCCGCGAACTTGCGCGCGTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAg                                 <  2:785390/93‑1 (MQ=255)

                                                          aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGa                               >  1:4294861/1‑93 (MQ=255)

                                                           ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGAc                              <  2:522126/93‑1 (MQ=255)

                                                           ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGAc                              <  2:91482/93‑1 (MQ=255)

                                                           ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGAc                              <  1:3814725/93‑1 (MQ=255)

                                                             gTCAGCCGCGAACTTGCGCGCTGTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTg                            <  1:1104841/93‑1 (MQ=255)

                                                              tCAGCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGc                           <  2:5104897/93‑1 (MQ=255)

                                                                 gCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTc                        >  1:3591319/1‑93 (MQ=255)

                                                                 gCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTc                        >  2:5613605/1‑93 (MQ=255)

                                                                 gCCGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTc                        <  2:1458782/93‑1 (MQ=255)

                                                                  ccGCGAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCt                       <  2:4436801/93‑1 (MQ=255)

                                                                   cgcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                      <  1:5558091/93‑1 (MQ=255)

                                                                   cgcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                      <  1:225650/93‑1 (MQ=255)

                                                                   cgcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTg                      <  2:2501720/93‑1 (MQ=255)

                                                                    gcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGt                     >  2:1337349/1‑93 (MQ=255)

                                                                    gcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGt                     <  1:645482/93‑1 (MQ=255)

                                                                    gcgAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGt                     <  1:2522800/93‑1 (MQ=255)

                                                                      gAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGttt                   >  2:35815/1‑93 (MQ=255)

                                                                      gAACTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGttt                   >  2:2212941/1‑93 (MQ=255)

                                                                       aaCTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTg                  >  1:901804/1‑93 (MQ=255)

                                                                       aaCTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTg                  <  2:5781522/93‑1 (MQ=255)

                                                                        aCTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                 >  1:2100094/1‑93 (MQ=255)

                                                                        aCTTGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTga                 <  2:2849835/93‑1 (MQ=255)

                                                                          ttGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTgaag               <  1:3461566/93‑1 (MQ=255)

                                                                          ttGCGCGCTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTgaag               >  1:1463405/1‑93 (MQ=255)

                                                                            gcgcgcTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa             <  1:4148234/93‑1 (MQ=255)

                                                                            gcgcgcTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa             <  2:4230868/93‑1 (MQ=255)

                                                                            gcgcgcTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa             >  2:6376465/1‑93 (MQ=255)

                                                                                gcTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGc         <  1:4888986/93‑1 (MQ=255)

                                                                                gcTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGc         <  1:1365222/93‑1 (MQ=255)

                                                                                 cTTTGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCa        >  2:5140670/1‑93 (MQ=255)

                                                                                  tttGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa       <  2:3020013/93‑1 (MQ=255)

                                                                                    tGAGAAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACg     <  1:5433143/93‑1 (MQ=255)

                                                                                     gagaAAGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGa    >  1:1949888/1‑93 (MQ=255)

GGTAATGACCGCCGGAAGCCTGACTGAGGCGGCACACCTGCTACACACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATC  >  NC_000913/2977084‑2977263

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