BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	804,639	C→T	G646G (GGC→GGT)	ybhJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	804,639	0	C	T	99.3%	519.9 / ‑5.4	152	G646G (GGC→GGT)	ybhJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (79/72); total (80/72)
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.95e-01

Only 100 of 158 total aligned reads displayed.

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

gtTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGa                                                                                         <  1:4496123/93‑1 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       >  1:1087667/1‑93 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       <  1:3040821/93‑1 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       >  1:2097198/1‑93 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  2:2695863/93‑1 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  1:2598280/93‑1 (MQ=255)

        gccttAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 <  2:964650/90‑1 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  1:2250929/1‑93 (MQ=255)

         cATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgc                                                                                <  1:4220377/93‑1 (MQ=255)

         cATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgc                                                                                <  2:3961564/93‑1 (MQ=255)

         cATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgc                                                                                <  2:2091165/93‑1 (MQ=255)

          aTTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgcg                                                                               >  2:3190821/1‑93 (MQ=255)

              aGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTg                                                                           >  1:2139281/1‑93 (MQ=255)

                 aGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGCTCTGCGCGTGAAc                                                                        >  2:1891124/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  2:3952421/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       <  1:245479/93‑1 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       <  1:3386323/93‑1 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  2:1108319/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  1:2896837/1‑93 (MQ=255)

                   aTTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTAGCAGCATGGCCTATGCGGTGAAACCAGGTGATGGTGCTGCGCGTGAACAg                                                                      >  2:2348141/1‑93 (MQ=255)

                    ttGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAgg                                                                     >  1:1751753/1‑93 (MQ=255)

                     tGCTGGTCAGGAGGATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  1:900269/93‑1 (MQ=255)

                     tGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  2:3001513/93‑1 (MQ=255)

                     tGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  2:1466728/93‑1 (MQ=255)

                     tGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  1:903418/93‑1 (MQ=255)

                     tGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  2:1466729/93‑1 (MQ=255)

                           tCAGGAGCATATTGATCCGCTGCAAAGTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAg                                                              >  2:4208877/1‑93 (MQ=255)

                           tCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAg                                                              >  2:2853725/1‑93 (MQ=255)

                            cAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGc                                                             <  2:4310744/93‑1 (MQ=255)

                             aGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                                                            <  1:3565377/93‑1 (MQ=255)

                                  cATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAg                                                       >  1:1901702/1‑93 (MQ=255)

                                   atatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGc                                                      >  1:2912275/1‑93 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     >  2:2184/1‑93 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     >  1:2142208/1‑93 (MQ=255)

                                      ttGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg                                                   >  2:2172367/1‑93 (MQ=255)

                                      ttGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg                                                   >  1:3924114/1‑93 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  <  1:3456122/93‑1 (MQ=255)

                                         aTCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa                                                <  1:1957781/93‑1 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               >  2:4176671/1‑93 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               >  2:4112573/1‑93 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:943119/93‑1 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:3203907/93‑1 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               >  1:3327285/1‑93 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               >  1:3524442/1‑93 (MQ=255)

                                          tCCGCTGCAAACTGAAAATGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:3959068/93‑1 (MQ=255)

                                           ccGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt                                              <  1:4485496/93‑1 (MQ=255)

                                           ccGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt                                              <  1:4587688/93‑1 (MQ=255)

                                            cGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg                                             <  2:443200/93‑1 (MQ=255)

                                             gCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg                                            >  1:3829288/1‑93 (MQ=255)

                                             gCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg                                            >  2:3199656/1‑93 (MQ=255)

                                              cTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc                                           <  1:470177/93‑1 (MQ=255)

                                               tGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCg                                          <  1:3442800/93‑1 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         >  2:2964204/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         >  2:1257993/1‑93 (MQ=255)

                                                  aaaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc                                       <  2:4597939/93‑1 (MQ=255)

                                                  aaaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc                                       >  1:4117901/1‑93 (MQ=255)

                                                  aaaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc                                       >  1:1451051/1‑93 (MQ=255)

                                                     cTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTgg                                    >  2:3327297/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTTGCGGTCTGGc                                   >  1:1053681/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   <  2:1191426/93‑1 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   <  1:731713/93‑1 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   >  2:260335/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   <  2:2515900/93‑1 (MQ=255)

                                                        aaaTTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGa                                 >  1:2342704/1‑93 (MQ=255)

                                                        aaaTTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGa                                 <  2:2518772/93‑1 (MQ=255)

                                                         aaTTGGCAGCATGGCCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGaa                                >  2:2816742/1‑93 (MQ=255)

                                                            tGGCAGCATGGTCTATGCGGTGAACCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAatat                             <  1:4465/93‑1 (MQ=255)

                                                            tGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAatat                             <  2:7882/93‑1 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            <  2:3150105/93‑1 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            <  2:134758/93‑1 (MQ=255)

                                                                agcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGcg                         >  1:1481275/1‑92 (MQ=255)

                                                                 gcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGGGGGATTGGCGGTCTGGCGAATATTGCCg                        >  2:3749089/1‑93 (MQ=255)

                                                                  cagggTCTATGCGGTGAACCCAGGTGATGGTTCTGCGCGTGAACACGCGGCGAGCTGCCAGCGTGTGATTCGCGGTCTGGCGAATATTGCCga                       <  1:2106551/90‑1 (MQ=255)

                                                                  caTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCga                       >  1:1399906/1‑93 (MQ=255)

                                                                  caTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCga                       <  1:531323/93‑1 (MQ=255)

                                                                   aTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgag                      <  2:1617679/93‑1 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  2:733461/1‑93 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  1:290720/1‑93 (MQ=255)

                                                                     tgTCTATACGGTGAATCCAGGTGATGGTTCTGCCCGCGAACAGGCGGCGACCTGCCAGCGTGTGATTGGCAGTCTGGCGAATATTGCCgagga                    <  2:660900/92‑1 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGGGAATATTGCCgaggg                    >  2:821921/1‑92 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  2:1720564/93‑1 (MQ=255)

                                                                      gTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgaggag                   >  1:3834939/1‑93 (MQ=255)

                                                                        cTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTa                 <  1:1251963/93‑1 (MQ=255)

                                                                        cTATGCGGGGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTa                 <  1:3596011/93‑1 (MQ=255)

                                                                          aTGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcg               >  1:3031951/1‑93 (MQ=255)

                                                                          aTGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGGCTGGCGAATATTGCCGAGGAGTAcg               >  1:2453755/1‑93 (MQ=255)

                                                                           tgtggTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcgc              <  2:3092630/90‑1 (MQ=255)

                                                                           tGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcgc              >  2:1039110/1‑93 (MQ=255)

                                                                              ggTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGAc           <  2:1998303/93‑1 (MQ=255)

                                                                                tGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTa         >  1:2199502/1‑93 (MQ=255)

                                                                                  aaaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTaaa       <  2:1121182/93‑1 (MQ=255)

                                                                                  aaaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGGACGCGACTaaa       >  2:4258907/1‑93 (MQ=255)

                                                                                   aaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGACTATTGCCGAGGAGTACGCGACTAAAc      >  1:629138/1‑93 (MQ=255)

                                                                                   aaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAAc      <  1:4294706/93‑1 (MQ=255)

                                                                                   aaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAAc      >  1:2869180/1‑93 (MQ=255)

                                                                                   aaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAAc      >  2:265408/1‑93 (MQ=255)

                                                                                     ccAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGc    <  1:4156961/93‑1 (MQ=255)

                                                                                      cAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCt   >  2:537842/1‑93 (MQ=255)

                                                                                      cAGGTGATGGTTCTGCGCGTGAACACGCGGCGAGCCGCCAGCGTGTGATTCGCGGTCTGGCGACTATTGCCGAGGAGTACGCGGCTAAACGCt   >  2:3624329/1‑93 (MQ=255)

                                                                                       aGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTa  <  1:472015/93‑1 (MQ=255)

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

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