BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,025,676	+C	coding (535/570 nt)	ygfQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,025,675	1	.	C	99.7%	993.3 / ‑4.1	288	coding (534/570 nt)	ygfQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (140/147); total (141/147)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 328 total aligned reads displayed.

GCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑CAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTCGCG  >  NC_000913/3025583‑3025762

gCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCGGCTAACGGTGTGGGTAAGGTTATCAAAAACCCGATTGAAGGCTTg                                                                                          >  1:7805718/1‑93 (MQ=255)

  tCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑cc                                                                                       >  2:3387674/1‑92 (MQ=255)

  tCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑cc                                                                                       >  2:673113/1‑92 (MQ=255)

   cACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑cca                                                                                      >  2:1456613/1‑91 (MQ=255)

    aCGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑ccag                                                                                     >  1:2015699/1‑90 (MQ=255)

    aCGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑ccag                                                                                     <  1:2770752/93‑4 (MQ=255)

       gTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTgg                                                                                   <  2:5243148/93‑1 (MQ=255)

       gTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGCGTGGGTATGGTTATCAAAAACCCGCTTGAAGGCTTGCCAGTgg                                                                                   >  2:6519824/1‑93 (MQ=255)

        tACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGGATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGgc                                                                                  <  2:4224171/93‑1 (MQ=255)

          cAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGgcgc                                                                                >  1:8380072/1‑93 (MQ=255)

          cAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGgcgc                                                                                >  2:743006/1‑93 (MQ=255)

                tcggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg                                                                          >  1:4753/1‑93 (MQ=255)

                taggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg                                                                          <  1:8103558/91‑1 (MQ=255)

                 cggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCgcgtc                                                                         >  2:2044342/1‑89 (MQ=255)

                 cggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGc                                                                         <  2:3527599/93‑1 (MQ=255)

                 cggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGc                                                                         >  2:7645857/1‑93 (MQ=255)

                 cggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGc                                                                         <  1:1044690/93‑1 (MQ=255)

                  ggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCg                                                                        >  1:7744922/1‑93 (MQ=255)

                   gtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGt                                                                       <  2:6318675/93‑1 (MQ=255)

                    tatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGtt                                                                      <  1:7172962/93‑1 (MQ=255)

                      tcggGCTGTTTCTTCTGCTGATTGCTTCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTa                                                                    <  1:2329664/93‑1 (MQ=255)

                      tcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTa                                                                    <  2:7175322/93‑1 (MQ=255)

                        gggCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcc                                                                  <  1:5378550/93‑1 (MQ=255)

                         ggCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                                 >  1:8783517/1‑93 (MQ=255)

                         ggCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                                 >  2:7893005/1‑93 (MQ=255)

                          gCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctc                                                                <  1:853239/93‑1 (MQ=255)

                          gCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctc                                                                <  2:2104763/93‑1 (MQ=255)

                          gCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctc                                                                >  2:6093333/1‑93 (MQ=255)

                           cTGGTTCTGCTGCGGATTGCTGCTAACGGTTTCGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCAGTGCGTTTAcctcc                                                               <  1:7624366/93‑1 (MQ=255)

                            tGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcct                                                              <  2:4310199/93‑1 (MQ=255)

                             gTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCtt                                                             <  2:3044321/93‑1 (MQ=255)

                                 ctgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCg                                                         <  2:3208890/93‑1 (MQ=255)

                                  tgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCgg                                                        <  1:5436951/93‑1 (MQ=255)

                                  tgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCgg                                                        >  2:6822873/1‑93 (MQ=255)

                                  tgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGCGCGTTTACCTCCTTCCCgg                                                        >  2:500591/1‑93 (MQ=255)

                                     tgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtga                                                     >  1:954371/1‑93 (MQ=255)

                                     tgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtga                                                     >  2:3923529/1‑93 (MQ=255)

                                        tgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatga                                                  >  1:6883554/1‑93 (MQ=255)

                                         gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgat                                                 >  2:1771275/1‑93 (MQ=255)

                                         gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgat                                                 >  1:6974254/1‑93 (MQ=255)

                                         gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgat                                                 <  2:1134565/93‑1 (MQ=255)

                                         gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgat                                                 >  1:6403139/1‑93 (MQ=255)

                                          aTTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGCGGCGCTCGGCGCGTTTACCTCCTTCCCGGcgatgatg                                                >  1:1268316/1‑93 (MQ=255)

                                             gctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGc                                             <  2:5799943/93‑1 (MQ=255)

                                              ctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCt                                            <  2:4680029/93‑1 (MQ=255)

                                              ctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCt                                            <  2:8637193/93‑1 (MQ=255)

                                               tgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                                           <  1:753815/93‑1 (MQ=255)

                                               tgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                                           <  2:1966407/93‑1 (MQ=255)

                                               tgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                                           <  1:5702173/93‑1 (MQ=255)

                                                 ctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGc                                         <  1:3098317/93‑1 (MQ=255)

                                                  tAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCt                                        >  1:4985605/1‑93 (MQ=255)

                                                  tAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCt                                        >  2:6184344/1‑93 (MQ=255)

                                                   aaCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTg                                       >  2:7241328/1‑93 (MQ=255)

                                                    aCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTgg                                      >  2:3749555/1‑93 (MQ=255)

                                                     cGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTggg                                     <  1:379695/93‑1 (MQ=255)

                                                       gtgtGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGc                                   >  2:2943121/1‑93 (MQ=255)

                                                         gtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTg                                 >  2:5151198/1‑93 (MQ=255)

                                                         gtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTg                                 <  1:6071464/93‑1 (MQ=255)

                                                          tgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCGGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTgg                                <  2:4838770/93‑1 (MQ=255)

                                                          tgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTgg                                <  1:271148/93‑1 (MQ=255)

                                                           gggTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGc                               >  2:8067319/1‑93 (MQ=255)

                                                            ggTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCg                              <  2:4177928/93‑1 (MQ=255)

                                                             gTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCgg                             <  1:3482403/93‑1 (MQ=255)

                                                             gTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCgg                             >  2:8034845/1‑93 (MQ=255)

                                                              tATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGt                            >  2:8393752/1‑93 (MQ=255)

                                                               aTGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTc                           >  1:6647105/1‑93 (MQ=255)

                                                               aTGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTc                           >  2:5478044/1‑93 (MQ=255)

                                                                tGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCa                          >  1:288069/1‑93 (MQ=255)

                                                                tGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCa                          <  1:531123/93‑1 (MQ=255)

                                                                  gTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATc                        >  1:7963224/1‑93 (MQ=255)

                                                                  gTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGCGCGTTTACCTCCTTACCGGCGATGATGAGCTTGCTGGGGATGGCGGTCATc                        >  1:6083237/1‑93 (MQ=255)

                                                                   ttATCAAAAACCCGATTGAAGGCTTGCCAGTGTCGCTCGGTCCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGATGGCGGTCATCt                       >  2:5251685/1‑93 (MQ=255)

                                                                    tATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCtt                      >  1:7255151/1‑93 (MQ=255)

                                                                    tATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCtt                      >  2:6954541/1‑93 (MQ=255)

                                                                      tCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCg                    <  2:2311065/93‑1 (MQ=255)

                                                                       cAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCgg                   >  2:4020489/1‑93 (MQ=255)

                                                                        aaaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGc                  <  2:1280443/93‑1 (MQ=255)

                                                                        aaaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGc                  <  1:3292470/93‑1 (MQ=255)

                                                                         aaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGcc                 <  1:4052227/93‑1 (MQ=255)

                                                                         aaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGcc                 <  2:6039140/93‑1 (MQ=255)

                                                                          aaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTCTACCTCCTTCCCGGTGCTGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCt                >  1:1798232/1‑93 (MQ=255)

                                                                             cccGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGga             <  1:6367545/93‑1 (MQ=255)

                                                                              ccGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGgag            >  2:844934/1‑93 (MQ=255)

                                                                              ccGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGgag            <  1:1869963/93‑1 (MQ=255)

                                                                                gATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGgagag          >  2:2574164/1‑92 (MQ=255)

                                                                                gATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGaa          >  2:6894805/1‑93 (MQ=255)

                                                                                 aTTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAg         >  1:1485717/1‑93 (MQ=255)

                                                                                 aTTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAg         >  2:4659705/1‑93 (MQ=255)

                                                                                  ttGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgt        >  2:3300001/1‑93 (MQ=255)

                                                                                  ttGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgt        <  1:3649809/93‑1 (MQ=255)

                                                                                  ttGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgt        <  2:7518411/93‑1 (MQ=255)

                                                                                  ttGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGGGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGCAgt        >  1:442806/1‑93 (MQ=255)

                                                                                    gAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGCAgagt      >  1:4016048/1‑93 (MQ=255)

                                                                                    gAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgtgt      <  2:7714542/93‑1 (MQ=255)

                                                                                    gAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgtgt      <  2:5323910/93‑1 (MQ=255)

                                                                                    gAAGACTTGCCAGTGGCGCTCGGCGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAgtgt      <  2:4759216/93‑1 (MQ=255)

                                                                                     aaGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTc     <  1:5826334/93‑1 (MQ=255)

                                                                                     aaGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTc     <  1:818995/93‑1 (MQ=255)

                                                                                     aaGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTc     <  1:1108401/93‑1 (MQ=255)

                                                                                        gCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTcgcg  >  1:510935/1‑93 (MQ=255)

GCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑CAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGGCCTGGAGAAGTGTCGCG  >  NC_000913/3025583‑3025762

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