BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	204,556	C→G	P22A (CCG→GCG)	rnhB →	RNAse HII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	204,556	0	C	G	97.5%	378.2 / 1.1	121	P22A (CCG→GCG)	rnhB	RNAse HII
Reads supporting (aligned to +/- strand): ref base C (0/3); new base G (56/62); total (56/65)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

Only 100 of 130 total aligned reads displayed.

GGCACAAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGCCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCG  >  NC_000913/204466‑204643

   acaAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTa                                                                                    >  1:3452873/1‑93 (MQ=255)

   acaAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTa                                                                                    >  1:2680279/1‑93 (MQ=255)

   acaAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTa                                                                                    >  1:2680281/1‑93 (MQ=255)

    caAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAg                                                                                   >  2:2798131/1‑93 (MQ=255)

      aGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGtt                                                                                 >  1:4204217/1‑93 (MQ=255)

       gCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTg                                                                                <  2:120875/93‑1 (MQ=255)

       gCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTg                                                                                >  1:2038343/1‑93 (MQ=255)

       gCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTg                                                                                <  2:1244477/93‑1 (MQ=255)

              tGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCggg                                                                         >  1:508462/1‑92 (MQ=255)

              tGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGgt                                                                         >  1:1115848/1‑93 (MQ=255)

                gAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGgtcg                                                                       >  1:1194223/1‑93 (MQ=255)

                  gTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGgtcgtc                                                                     <  2:3865457/93‑1 (MQ=255)

                  gTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGgtcgtc                                                                     <  2:3525561/93‑1 (MQ=255)

                  gTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGgtcgtc                                                                     <  1:2918680/93‑1 (MQ=255)

                   ttAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCa                                                                    >  2:2087225/1‑93 (MQ=255)

                    tAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCAc                                                                   <  1:4389051/93‑1 (MQ=255)

                     aGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCAcc                                                                  >  1:1803575/1‑93 (MQ=255)

                      gCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCg                                                                 >  1:983851/1‑93 (MQ=255)

                       cacaATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGTGCGTTAGTTGGCGCGGTCGTCACCGc                                                                <  1:4172519/93‑1 (MQ=255)

                       cacaATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGGGGTCGTCACCGc                                                                >  1:1385674/1‑93 (MQ=255)

                       cacaATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGc                                                                <  1:3269359/93‑1 (MQ=255)

                        acaATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGGGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCCGCACCGCt                                                               >  1:411757/1‑93 (MQ=255)

                          caTGCTCGATTTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCGCCGCTGc                                                             <  2:4365519/92‑1 (MQ=255)

                          aaTGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTgg                                                             >  1:608439/1‑92 (MQ=255)

                          aaTGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGc                                                             <  1:836237/93‑1 (MQ=255)

                           aTGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCg                                                            >  2:3465886/1‑93 (MQ=255)

                           aTGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCg                                                            <  2:1422732/93‑1 (MQ=255)

                           aTGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCg                                                            >  1:4305997/1‑93 (MQ=255)

                            tGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCgg                                                           >  2:2269674/1‑93 (MQ=255)

                             gATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGt                                                          <  1:3423699/93‑1 (MQ=255)

                             gATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGt                                                          <  1:3115855/93‑1 (MQ=255)

                              aTCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTg                                                         <  1:4407340/93‑1 (MQ=255)

                              aTCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTg                                                         <  2:653372/93‑1 (MQ=255)

                               tCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGa                                                        >  1:2037415/1‑93 (MQ=255)

                                cGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGAt                                                       <  1:824929/93‑1 (MQ=255)

                                 gAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATc                                                      <  1:550712/93‑1 (MQ=255)

                                   aTTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTTATCCt                                                    >  1:4720805/1‑93 (MQ=255)

                                   aTTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCt                                                    <  2:4033493/93‑1 (MQ=255)

                                   aTTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCt                                                    <  1:894319/93‑1 (MQ=255)

                                      tGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGa                                                 <  1:2683525/93‑1 (MQ=255)

                                       gTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGAc                                                >  2:3251875/1‑93 (MQ=255)

                                       gTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGAc                                                >  2:2454317/1‑93 (MQ=255)

                                       gTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGAc                                                <  1:3080330/93‑1 (MQ=255)

                                        tttATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGAcc                                               >  1:3442944/1‑93 (MQ=255)

                                        tttATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGGCACCGCTGCGGTGATCCTTGAcc                                               >  2:3372083/1‑93 (MQ=255)

                                          tATCAGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCg                                             >  2:3995545/1‑93 (MQ=255)

                                           aTCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCgg                                            >  1:4373992/1‑93 (MQ=255)

                                            tCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgc                                           >  1:2852286/1‑93 (MQ=255)

                                            tCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgc                                           <  2:4129089/93‑1 (MQ=255)

                                            tCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGCCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgc                                           <  1:3871125/93‑1 (MQ=255)

                                             ccGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgcg                                          <  2:2746671/93‑1 (MQ=255)

                                               gCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgcgcg                                        >  1:4389125/1‑93 (MQ=255)

                                               gCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGgcgcg                                        >  1:868827/1‑93 (MQ=255)

                                                 acacGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGcc                                      >  1:4420598/1‑93 (MQ=255)

                                                 acacGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGCCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGcc                                      <  1:4028040/93‑1 (MQ=255)

                                                  cacGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGccc                                     >  1:3804505/1‑93 (MQ=255)

                                                   ccggaGCTGGTTGCGGGTGGGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCg                                    <  2:1047874/89‑1 (MQ=255)

                                                   acGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGTCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCg                                    >  2:4609031/1‑93 (MQ=255)

                                                   acGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCg                                    >  2:3311262/1‑93 (MQ=255)

                                                   acGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCg                                    >  2:3949256/1‑93 (MQ=255)

                                                    cGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGa                                   >  1:1935610/1‑93 (MQ=255)

                                                     gCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGAt                                  >  1:524686/1‑93 (MQ=255)

                                                     gCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGAt                                  <  1:3193766/93‑1 (MQ=255)

                                                     gCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGAt                                  <  1:3330562/93‑1 (MQ=255)

                                                        gCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGc                               >  1:4590553/1‑93 (MQ=255)

                                                          tGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCg                             <  2:2037067/93‑1 (MQ=255)

                                                          tGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCg                             <  2:269289/93‑1 (MQ=255)

                                                          tGGTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCg                             <  2:3024889/93‑1 (MQ=255)

                                                           ggTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCgg                            <  2:4248161/93‑1 (MQ=255)

                                                           ggTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCgg                            <  1:3672311/93‑1 (MQ=255)

                                                            gTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCggg                           >  2:2349604/1‑93 (MQ=255)

                                                            gTTGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCggg                           <  2:363708/93‑1 (MQ=255)

                                                              tGCGGGTGTGGATGAAGTCGGCCGAGGGCCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCt                         <  2:2779757/93‑1 (MQ=255)

                                                              tGCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCt                         <  2:1811687/93‑1 (MQ=255)

                                                               gCGGGTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTg                        <  1:1722082/93‑1 (MQ=255)

                                                                 gggTGTGGGTGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCAGCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGaa                      <  1:3659980/93‑1 (MQ=255)

                                                                 gggTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGaa                      <  1:491494/93‑1 (MQ=255)

                                                                 gggTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGaa                      <  1:4720860/93‑1 (MQ=255)

                                                                  ggTGTGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAAt                     >  1:508783/1‑93 (MQ=255)

                                                                   gtgtGGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATg                    <  2:4592904/93‑1 (MQ=255)

                                                                    tgtgGATGAAGTCGGACGCTGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGa                   <  2:2414597/93‑1 (MQ=255)

                                                                    tgtgGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGa                   <  1:162574/93‑1 (MQ=255)

                                                                     gtgGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGAGTGCCGGGCTGAATGAt                  >  2:4406445/1‑93 (MQ=255)

                                                                      tgGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGAtt                 <  2:461263/93‑1 (MQ=255)

                                                                      tgGATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGAtt                 >  2:2515594/1‑93 (MQ=255)

                                                                       ggATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTc                >  2:4749892/1‑93 (MQ=255)

                                                                        gATGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTcc               <  2:714023/93‑1 (MQ=255)

                                                                         aTGAAGTCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCa              <  2:4376961/93‑1 (MQ=255)

                                                                               tCGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGc        <  1:341413/93‑1 (MQ=255)

                                                                                cGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCt       <  1:4398043/93‑1 (MQ=255)

                                                                                cGGACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCt       <  1:1327143/93‑1 (MQ=255)

                                                                                cGGACGCGGCGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCt       <  2:3292941/93‑1 (MQ=255)

                                                                                 ggACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTg      <  2:1792703/93‑1 (MQ=255)

                                                                                 ggACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTg      <  2:2362395/93‑1 (MQ=255)

                                                                                  gACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATTCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGa     <  2:1222475/93‑1 (MQ=255)

                                                                                  gACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGa     >  2:3761853/1‑93 (MQ=255)

                                                                                  gACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGa     <  2:3050753/93‑1 (MQ=255)

                                                                                  gACGCGGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGa     <  2:622448/93‑1 (MQ=255)

                                                                                   aCGCGGGGCGTTAGTTGGCGCGGTAGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAg    <  1:3949743/93‑1 (MQ=255)

                                                                                     gcgGGGCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCg  <  1:647098/93‑1 (MQ=255)

GGCACAAGCCGTTCTGGAGTTAGCACAATGATCGAATTTGTTTATCCGCACACGCAGCTGGTTGCGGGTGTGGATGAAGTCGGACGCGGGCCGTTAGTTGGCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCG  >  NC_000913/204466‑204643

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