BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,338,481	C→T	A189V (GCC→GTC)	yrbG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,338,481	0	C	T	95.7%	1132.9 / 16.7	346	A189V (GCC→GTC)	yrbG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (5/10); new base T (122/208); total (128/218)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 409 total aligned reads displayed.

ACCAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGCCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATCGGA  >  NC_000913/3338390‑3338572

gccAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGcc                                                                                            <  1:6502749/92‑1 (MQ=255)

 ccAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCa                                                                                           <  2:5015671/93‑1 (MQ=255)

  cAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGGATTGCGCTTATCATCATGCCAGTGGTCAc                                                                                          >  1:2528508/1‑93 (MQ=255)

   agagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGGGCTTATCATCATGCCCGTGATCAcg                                                                                         >  2:5240330/1‑93 (MQ=255)

   agagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCAGCATTGCGCTTATCATCATGCCAGTGGTCAcg                                                                                         <  1:9604895/93‑1 (MQ=255)

    gagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGGATTGCGCTTATCATCATGCCAGTGGGCAcga                                                                                        >  1:1425247/1‑92 (MQ=255)

          agcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtgg                                                                                  <  2:3137642/93‑1 (MQ=255)

          agcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtgg                                                                                  >  2:3402913/1‑93 (MQ=255)

          agcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtgg                                                                                  <  2:8098130/93‑1 (MQ=255)

           gcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGGGCTTATCATCATGCCAGTGGTCACGCGGAtggg                                                                                 >  1:9852292/1‑92 (MQ=255)

           gcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGGtggg                                                                                 >  1:3031950/1‑92 (MQ=255)

           gcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtggt                                                                                 <  2:2612431/93‑1 (MQ=255)

            cTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtggtg                                                                                <  1:2674177/93‑1 (MQ=255)

              tGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATggtggt                                                                              >  1:6569828/1‑93 (MQ=255)

              tGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATggtggt                                                                              <  2:1553263/93‑1 (MQ=255)

               gCAGAGCTGCCGCGTGACGGCGGATTTCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGGTGGGGGtt                                                                             >  2:4963526/1‑93 (MQ=255)

               gCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGtt                                                                             <  2:7362028/93‑1 (MQ=255)

               gCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGtt                                                                             <  1:8280852/93‑1 (MQ=255)

               gCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCCGGATTGTGCTTTTCATCTTGTCCGTGGTCAAGGCGAtggggggg                                                                             >  2:8491300/1‑91 (MQ=255)

                cAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCAATGCGCTTATCATCATGCCAGTGGTCACGCGGCTGGTGGttt                                                                            >  2:5616293/1‑92 (MQ=255)

                  gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGAt                                                                          <  2:9364239/93‑1 (MQ=255)

                  gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCCGATggtgggtggt                                                                          >  1:4783260/1‑91 (MQ=255)

                   agCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATa                                                                         <  1:8496034/93‑1 (MQ=255)

                   agCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATa                                                                         <  2:2903024/93‑1 (MQ=255)

                    gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCCTTGCGCTTATCATCATTCCAGTGGTCACGCGGATGGTGGTTGATaa                                                                        >  1:104442/1‑93 (MQ=255)

                    gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATaa                                                                        <  1:3278512/93‑1 (MQ=255)

                    gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATaa                                                                        <  1:7939026/93‑1 (MQ=255)

                     cTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAAc                                                                       <  1:8964464/93‑1 (MQ=255)

                       gCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGGTGGTGGGTGATAACGc                                                                     >  1:5303639/1‑93 (MQ=255)

                       gCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGc                                                                     <  1:1313574/93‑1 (MQ=255)

                        ccGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGcc                                                                    <  2:1409577/93‑1 (MQ=255)

                        ccGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGcc                                                                    <  2:5753123/93‑1 (MQ=255)

                          gcgTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCAc                                                                  <  1:9701324/93‑1 (MQ=255)

                          gcgTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGGGGTTGATAACGCCAc                                                                  >  1:550120/1‑93 (MQ=255)

                            gTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAACGCCACgg                                                                >  1:9468647/1‑93 (MQ=255)

                             tGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAACGCCACCGt                                                               >  1:6126206/1‑93 (MQ=255)

                              gACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCATGGGGCACGCGGAGCGTGGTTGATAACGCCACgggg                                                              >  2:6293134/1‑91 (MQ=255)

                               aCGGCGGATTGCCCGTCGCGTTTTTATGGCTCCGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGGGGGTGATAACGCCACGGTGc                                                             >  2:105549/1‑93 (MQ=255)

                                 ggcggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAACGCCACGGTGCTg                                                           >  2:714272/1‑93 (MQ=255)

                                  gcggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGGGGTTGATAACGCCACGGTGCTgg                                                          >  1:741146/1‑93 (MQ=255)

                                    ggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGGTTATCATCATGCCAGTGGTGACGCGGATGGTGGTTGATAAAGGCACGGGGCTggga                                                        >  2:4659184/1‑91 (MQ=255)

                                       ttGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAAt                                                     <  2:9818789/93‑1 (MQ=255)

                                         gCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTa                                                   <  1:5158961/93‑1 (MQ=255)

                                         gCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGGGGTCACGCGGATGGTGGGTGATAACGCCACGGTGCTGGCGAATTa                                                   >  1:3489037/1‑93 (MQ=255)

                                          cccGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTAc                                                  <  1:4338719/93‑1 (MQ=255)

                                           ccGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACt                                                 >  2:238750/1‑93 (MQ=255)

                                           ccGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACt                                                 <  1:5638692/93‑1 (MQ=255)

                                           ccGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAAAGCCACGGTGCTGGCGGATTACt                                                 >  2:8806350/1‑93 (MQ=255)

                                            cGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACtt                                                >  2:4570737/1‑93 (MQ=255)

                                            cGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACtt                                                >  2:3736266/1‑93 (MQ=255)

                                            cGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACtt                                                <  1:4909435/93‑1 (MQ=255)

                                              tCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTg                                              >  2:2166235/1‑93 (MQ=255)

                                              tCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTg                                              >  2:1212548/1‑93 (MQ=255)

                                                gcgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGcc                                            <  2:8268288/93‑1 (MQ=255)

                                                 cgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCa                                           <  2:2485208/93‑1 (MQ=255)

                                                 cgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCa                                           <  2:5695149/93‑1 (MQ=255)

                                                 cgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCa                                           <  1:7888880/93‑1 (MQ=255)

                                                   tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                         <  2:65911/93‑1 (MQ=255)

                                                   tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                         <  2:9473701/93‑1 (MQ=255)

                                                          ggCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGt                                  <  2:4214025/93‑1 (MQ=255)

                                                             tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                               <  2:6894548/93‑1 (MQ=255)

                                                             tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                               <  1:2253283/93‑1 (MQ=255)

                                                             tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                               >  2:9776165/1‑93 (MQ=255)

                                                             tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                               <  1:7550878/93‑1 (MQ=255)

                                                                gCATTGCGCTTGTCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGa                            <  2:4387760/93‑1 (MQ=255)

                                                                gCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGa                            <  2:5853558/93‑1 (MQ=255)

                                                                gCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGa                            >  1:3802537/1‑93 (MQ=255)

                                                                  aTTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATg                          >  2:2835525/1‑93 (MQ=255)

                                                                  aTTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATg                          >  1:2384556/1‑93 (MQ=255)

                                                                  aTTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATg                          >  1:5878393/1‑93 (MQ=255)

                                                                   ttGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATgg                         >  1:213541/1‑93 (MQ=255)

                                                                   ttGCGCTTATCATCATGCCAGTAGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATgg                         >  2:9666615/1‑93 (MQ=255)

                                                                   ttGCGCTTATCATCATGCCAGGGGTCACGCGGATGGGGGGTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGCTGAGGAggg                         >  1:7733918/1‑93 (MQ=255)

                                                                    tGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATggg                        <  1:10001073/93‑1 (MQ=255)

                                                                     gcgcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGt                       >  1:8748725/1‑93 (MQ=255)

                                                                      cgcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTc                      <  2:1333182/93‑1 (MQ=255)

                                                                      cgcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTc                      <  1:8103366/93‑1 (MQ=255)

                                                                       gcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCt                     >  2:3513189/1‑93 (MQ=255)

                                                                         ttATCATCATTCCAGCGGTCACGCGGATTGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGa                   <  1:2763926/93‑1 (MQ=255)

                                                                         ttATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGa                   <  1:4051851/93‑1 (MQ=255)

                                                                          tATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGAc                  <  2:4317656/93‑1 (MQ=255)

                                                                             catcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGc               >  1:6168835/1‑93 (MQ=255)

                                                                               tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACTATGGGTCTGACGGCaa             >  1:6826317/1‑93 (MQ=255)

                                                                               tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCaa             <  2:7503854/93‑1 (MQ=255)

                                                                               tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCaa             >  2:5095450/1‑93 (MQ=255)

                                                                               tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCaa             <  1:628993/93‑1 (MQ=255)

                                                                               tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCaa             >  2:1667873/1‑93 (MQ=255)

                                                                                catGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAAt            <  2:2234771/93‑1 (MQ=255)

                                                                                 atGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAAtt           <  2:5380371/93‑1 (MQ=255)

                                                                                  tGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTg          >  2:8925502/1‑93 (MQ=255)

                                                                                   gCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGc         >  1:1767776/1‑93 (MQ=255)

                                                                                    ccAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCt        >  1:42454/1‑93 (MQ=255)

                                                                                     cAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTa       >  1:994099/1‑93 (MQ=255)

                                                                                     cAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTa       >  1:4515979/1‑93 (MQ=255)

                                                                                     cAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTa       >  1:7138139/1‑93 (MQ=255)

                                                                                     cAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTa       >  2:412922/1‑93 (MQ=255)

                                                                                      aGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTAt      <  2:8688734/93‑1 (MQ=255)

                                                                                       gTGGTCCCGCGGATGGTGGTTGATAACGCCACGGTGCTGGCCAATTACTTTGCCATCAGCGAGTTGACGACGGGTCTGACGGCAATTGCTATc     <  1:1828934/93‑1 (MQ=255)

                                                                                       gTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATc     <  1:3686092/93‑1 (MQ=255)

                                                                                       gTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATc     >  2:6056884/1‑93 (MQ=255)

ACCAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGCCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATCGGA  >  NC_000913/3338390‑3338572

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