BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,542,788	A→G	S107G (AGC→GGC)	yhgH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,542,788	0	A	G	81.7%	665.3 / 157.3	323	S107G (AGC→GGC)	yhgH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (30/29); new base G (121/143); total (151/172)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 324 total aligned reads displayed.

TGGCAAAGACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTA  >  NC_000913/3542698‑3542880

tGGCAAAGACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGgc                                                                                            >  1:1101810/1‑93 (MQ=255)

     aaGACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCt                                                                                       >  2:7117730/1‑93 (MQ=255)

     aaGACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCt                                                                                       >  2:1466410/1‑93 (MQ=255)

       gACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGt                                                                                     <  2:6512869/93‑1 (MQ=255)

        aCTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTc                                                                                    <  1:876593/93‑1 (MQ=255)

         cTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCa                                                                                   >  2:704228/1‑93 (MQ=255)

           ggTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACg                                                                                 >  1:4498065/1‑93 (MQ=255)

            gTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGt                                                                                >  2:6665120/1‑93 (MQ=255)

             tCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTc                                                                               <  2:261413/93‑1 (MQ=255)

             tCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTc                                                                               <  2:4144913/93‑1 (MQ=255)

             tCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTc                                                                               <  1:6272764/93‑1 (MQ=255)

               aCGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTg                                                                             >  1:5894048/1‑93 (MQ=255)

               aCGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTg                                                                             >  1:1228009/1‑93 (MQ=255)

               aCGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTg                                                                             <  1:1170110/93‑1 (MQ=255)

                cGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGt                                                                            >  2:5652120/1‑93 (MQ=255)

                cGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGt                                                                            <  1:259197/93‑1 (MQ=255)

                 ggTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGtt                                                                           >  2:3501820/1‑93 (MQ=255)

                  gTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTa                                                                          <  2:7486028/93‑1 (MQ=255)

                  gTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTa                                                                          >  2:4633594/1‑93 (MQ=255)

                   ttGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCGCCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTAc                                                                         <  2:24655/93‑1 (MQ=255)

                   ttGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCGCCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTAc                                                                         <  1:2238687/93‑1 (MQ=255)

                      ccGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACttt                                                                      <  2:434362/93‑3 (MQ=255)

                      ccGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTAt                                                                      <  2:632001/93‑1 (MQ=255)

                      ccGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTAt                                                                      <  2:3890436/93‑1 (MQ=255)

                        gACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTg                                                                    >  2:3160822/1‑93 (MQ=255)

                         aCTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTgg                                                                   >  1:7030065/1‑93 (MQ=255)

                         aCTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTgg                                                                   <  2:1883171/93‑1 (MQ=255)

                         aCTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTgg                                                                   >  1:2084367/1‑93 (MQ=255)

                          cTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGa                                                                  <  2:2356881/93‑1 (MQ=255)

                          cTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGa                                                                  >  1:5449724/1‑93 (MQ=255)

                            aTGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAg                                                                <  2:5735943/93‑1 (MQ=255)

                             tGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGt                                                               <  2:4487510/93‑1 (MQ=255)

                             tGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGt                                                               <  1:3570893/93‑1 (MQ=255)

                              gcgcCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTc                                                              <  1:5191534/93‑1 (MQ=255)

                                gccgccGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCtt                                                            <  1:6103639/93‑1 (MQ=255)

                                 ccgccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTa                                                           <  1:5524561/93‑1 (MQ=255)

                                 ccgccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTa                                                           <  2:1693189/93‑1 (MQ=255)

                                 ccgccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTa                                                           <  2:2060821/93‑1 (MQ=255)

                                  cgccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTac                                                          <  2:2765941/93‑1 (MQ=255)

                                   gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTaca                                                         <  2:3645826/93‑1 (MQ=255)

                                   gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTaca                                                         <  1:96388/93‑1 (MQ=255)

                                   gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTaca                                                         >  1:3019886/1‑93 (MQ=255)

                                   gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTaca                                                         <  2:6936404/93‑1 (MQ=255)

                                   gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTaca                                                         <  1:6782969/93‑1 (MQ=255)

                                       ttAAGTCCGCTTATCCGCCAGCTTGAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCt                                                     >  2:7359424/1‑93 (MQ=255)

                                       ttAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCt                                                     >  1:392833/1‑93 (MQ=255)

                                       ttAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCt                                                     >  1:5784588/1‑93 (MQ=255)

                                       ttAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCt                                                     >  2:3279237/1‑93 (MQ=255)

                                        tAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtc                                                    <  2:5287201/93‑1 (MQ=255)

                                        tAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtc                                                    >  2:1551303/1‑93 (MQ=255)

                                        tAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtc                                                    >  1:796242/1‑93 (MQ=255)

                                         aaGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGACGTCTTACACGCtcg                                                   >  2:4378651/1‑93 (MQ=255)

                                          aGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtcgt                                                  <  2:3195396/93‑1 (MQ=255)

                                          aGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtcgt                                                  <  1:2006822/93‑1 (MQ=255)

                                          aGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCtcgt                                                  <  1:596410/93‑1 (MQ=255)

                                             ccGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGc                                               >  2:2940566/1‑93 (MQ=255)

                                             ccGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGc                                               >  1:7383920/1‑93 (MQ=255)

                                              cGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGca                                              >  1:3253167/1‑93 (MQ=255)

                                                  tATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGcaccac                                          <  2:5921950/93‑1 (MQ=255)

                                                  tATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGcaccac                                          >  2:4075720/1‑93 (MQ=255)

                                                  tATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGcaccac                                          >  2:810732/1‑93 (MQ=255)

                                                    tCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCg                                        <  1:1322874/93‑1 (MQ=255)

                                                    tCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCg                                        <  1:1527081/93‑1 (MQ=255)

                                                     ccaccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCgg                                       <  1:4615355/93‑1 (MQ=255)

                                                      caccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCggg                                      <  2:2902380/93‑1 (MQ=255)

                                                      caccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCggg                                      >  2:6096012/1‑93 (MQ=255)

                                                       accaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGt                                     <  1:2717944/93‑1 (MQ=255)

                                                        ccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGtt                                    <  1:5668648/93‑1 (MQ=255)

                                                         caGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTg                                   >  2:4697171/1‑93 (MQ=255)

                                                          aGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGc                                  <  2:2252376/93‑1 (MQ=255)

                                                          aGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGc                                  >  2:368374/1‑93 (MQ=255)

                                                          aGCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGc                                  >  2:108597/1‑93 (MQ=255)

                                                           gCTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCa                                 >  2:6827719/1‑93 (MQ=255)

                                                            cTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCaa                                <  1:4211853/93‑1 (MQ=255)

                                                            cTTAAATTTTCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCaa                                >  1:6838903/1‑93 (MQ=255)

                                                                 aTTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCCCGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGcc                           >  2:4952392/1‑93 (MQ=255)

                                                                  ttttCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCg                          >  1:1038886/1‑93 (MQ=255)

                                                                    ttCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGa                        >  2:6040199/1‑93 (MQ=255)

                                                                    ttCCCGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGa                        <  2:1394169/93‑1 (MQ=255)

                                                                      cccGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATc                      <  1:3727242/93‑1 (MQ=255)

                                                                        cGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGc                    <  2:3946907/93‑1 (MQ=255)

                                                                        cGGCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGc                    >  1:2346770/1‑93 (MQ=255)

                                                                         ggCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCa                   >  1:5403593/1‑93 (MQ=255)

                                                                         ggCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCa                   <  2:1246914/93‑1 (MQ=255)

                                                                         ggCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCa                   >  2:5089158/1‑93 (MQ=255)

                                                                         ggCGCAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCa                   >  1:6305115/1‑93 (MQ=255)

                                                                          gcgcAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCAt                  >  2:2548701/1‑93 (MQ=255)

                                                                          gcgcAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCAt                  <  1:3522833/93‑1 (MQ=255)

                                                                          gcgcAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCAt                  <  1:4897620/93‑1 (MQ=255)

                                                                            gcAGTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCg                <  1:3825609/93‑1 (MQ=255)

                                                                               gTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCa             >  1:6548904/1‑93 (MQ=255)

                                                                               gTGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCa             <  2:1832639/93‑1 (MQ=255)

                                                                                tGAAATCGCCGGCGCCCTGTCACGTCTGTTATTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAg            <  1:162900/93‑1 (MQ=255)

                                                                                tGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAg            <  2:1130727/93‑1 (MQ=255)

                                                                                tGAAATCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAg            <  2:5514774/93‑1 (MQ=255)

                                                                                  aaaTCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCg          >  2:3805974/1‑93 (MQ=255)

                                                                                  aaaTCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCg          >  2:7220925/1‑93 (MQ=255)

                                                                                   aaTCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGt         >  1:4081588/1‑93 (MQ=255)

                                                                                   aaTCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGt         <  1:7517877/93‑1 (MQ=255)

                                                                                     tCGCCGGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTc       <  1:5557640/93‑1 (MQ=255)

TGGCAAAGACTGGTCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTA  >  NC_000913/3542698‑3542880

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