BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,505	G→A	R462R (CGG→CGA)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,505	0	G	A	98.5%	460.0 / ‑4.6	136	R462R (CGG→CGA)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (62/72); total (64/72)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 140 total aligned reads displayed.

TGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAT  >  NC_000913/2208415‑2208592

tGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGAcc                                                                                       >  2:1664161/1‑93 (MQ=255)

  ccTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGTAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctc                                                                                     <  2:2036706/93‑1 (MQ=255)

   cTGGCAGATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctct                                                                                    <  1:3108818/93‑1 (MQ=255)

   cTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctct                                                                                    <  1:2956821/93‑1 (MQ=255)

    tGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCtt                                                                                   >  1:321626/1‑93 (MQ=255)

      gCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTaa                                                                                 <  2:3811764/93‑1 (MQ=255)

      gCAAATGTGGAGCGCGCCGTTACGCCAGCCGGGGATTGTGCCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTaa                                                                                 <  2:1756485/93‑1 (MQ=255)

       cAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAAt                                                                                <  1:3453464/93‑1 (MQ=255)

       cAAAGGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAAt                                                                                <  1:160093/93‑1 (MQ=255)

           tGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAgg                                                                            <  1:2203228/93‑1 (MQ=255)

           tGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAgg                                                                            >  2:2946849/1‑93 (MQ=255)

            tttGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                           <  2:3832668/92‑1 (MQ=255)

            gTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                           <  1:281262/93‑1 (MQ=255)

                agggcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGAt                                                                       <  1:1240459/90‑1 (MQ=255)

                aGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGAt                                                                       >  2:2944310/1‑93 (MQ=255)

                 gcgcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATa                                                                      >  2:1678302/1‑93 (MQ=255)

                 gcgcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATa                                                                      >  1:3908466/1‑93 (MQ=255)

                   gcgcCGTTCCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                    <  2:1268676/93‑1 (MQ=255)

                   gcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCTGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                    <  2:1541294/93‑1 (MQ=255)

                   gcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                    <  1:2562985/93‑1 (MQ=255)

                   cccccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                    <  1:1622297/90‑1 (MQ=255)

                      ccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtta                                                                 <  2:4156408/93‑1 (MQ=255)

                      ccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtta                                                                 >  1:1368745/1‑93 (MQ=255)

                       cGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttat                                                                >  1:1499607/1‑93 (MQ=255)

                        gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                               <  2:4677880/93‑1 (MQ=255)

                          tACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattta                                                             <  2:4456326/93‑1 (MQ=255)

                           aCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatttat                                                            <  2:4329369/93‑1 (MQ=255)

                             gCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTc                                                          <  1:66730/93‑1 (MQ=255)

                             gCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTc                                                          <  1:3494099/93‑1 (MQ=255)

                              ccAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCt                                                         <  1:2600985/93‑1 (MQ=255)

                              ccAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCt                                                         <  2:2223602/93‑1 (MQ=255)

                                aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                       >  1:282447/1‑93 (MQ=255)

                                aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                       <  2:1642796/93‑1 (MQ=255)

                                   cGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                    >  1:2586248/1‑93 (MQ=255)

                                   cGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                    >  1:182159/1‑93 (MQ=255)

                                   cGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                    >  2:136801/1‑93 (MQ=255)

                                   cGGGCATTGTGGCTTTGCGGAATTATTTCCGCCAGCTCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                    <  1:3371026/93‑1 (MQ=255)

                                    gggCATTGGGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGt                                                   <  1:3508404/93‑1 (MQ=255)

                                     ggCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTc                                                  >  1:3427556/1‑93 (MQ=255)

                                     ggCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTc                                                  >  2:3015360/1‑93 (MQ=255)

                                      gCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCg                                                 <  1:2997777/93‑1 (MQ=255)

                                      gCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCg                                                 >  2:1384741/1‑93 (MQ=255)

                                          tgtgGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGa                                             <  1:4342753/93‑1 (MQ=255)

                                            tgGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAgt                                           <  2:1048430/93‑1 (MQ=255)

                                             ggCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAgtg                                          >  2:22491/1‑93 (MQ=255)

                                              gCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTATGCCGGTCGCTGAgtgt                                         <  1:475165/93‑1 (MQ=255)

                                               cTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTa                                        <  1:227762/93‑1 (MQ=255)

                                                tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                       >  1:3414325/1‑93 (MQ=255)

                                                tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                       <  1:1375445/93‑1 (MQ=255)

                                                tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                       <  2:4171011/93‑1 (MQ=255)

                                                 ttGCGGGAATATTTCCGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAtt                                      <  2:3378481/93‑1 (MQ=255)

                                                     gggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGc                                  >  1:1977861/1‑93 (MQ=255)

                                                     gggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATGTCGc                                  >  1:2058145/1‑93 (MQ=255)

                                                      ggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCt                                 <  1:4053518/93‑1 (MQ=255)

                                                      ggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCt                                 >  2:148333/1‑93 (MQ=255)

                                                      ggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCt                                 <  1:744921/93‑1 (MQ=255)

                                                         atatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCTGTCGCTGAGTGTATTTCGCTCgg                              >  2:2608368/1‑93 (MQ=255)

                                                         atatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCgg                              <  1:3780704/93‑1 (MQ=255)

                                                         atatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCgg                              <  1:2204895/93‑1 (MQ=255)

                                                          tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                             <  1:733370/93‑1 (MQ=255)

                                                          tatTTACGCCAGCGCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                             <  1:4101415/93‑1 (MQ=255)

                                                           atTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGtt                            <  1:3922633/93‑1 (MQ=255)

                                                           atTTACGCCAGCGTCCCCCCGACTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGtt                            >  2:3990608/1‑93 (MQ=255)

                                                            tttACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTACTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTg                           >  1:4572293/1‑93 (MQ=255)

                                                               aCGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGa                        >  2:1882763/1‑93 (MQ=255)

                                                                cGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAc                       <  2:2566742/93‑1 (MQ=255)

                                                                cGCCAGCGCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAc                       <  1:52440/93‑1 (MQ=255)

                                                                 ccccgcgcCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTTGCTCGGTTGGGAca                      <  2:3564581/89‑1 (MQ=255)

                                                                   cAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagc                    <  2:103186/93‑1 (MQ=255)

                                                                    aGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGGTGGGAcagca                   >  1:1802590/1‑93 (MQ=255)

                                                                     gCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagcag                  >  2:2014741/1‑93 (MQ=255)

                                                                     gCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagcag                  >  1:4680905/1‑93 (MQ=255)

                                                                     gCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagcag                  >  1:1951361/1‑93 (MQ=255)

                                                                      cgcccccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagc                 >  2:1009051/4‑93 (MQ=255)

                                                                      cGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagc                 <  2:571506/93‑1 (MQ=255)

                                                                       gTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgcc                >  2:4392569/1‑93 (MQ=255)

                                                                          cccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgcc             >  2:1809345/1‑93 (MQ=255)

                                                                          cccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgcc             >  2:399698/1‑93 (MQ=255)

                                                                           ccccTCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa            <  2:3097075/93‑1 (MQ=255)

                                                                           ccccGCCTGCATCCGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCATCCGCCa            <  2:3097074/93‑1 (MQ=255)

                                                                           ccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa            <  1:215748/93‑1 (MQ=255)

                                                                            cccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAg           >  1:4678614/1‑93 (MQ=255)

                                                                             ccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGa          >  1:4359957/1‑93 (MQ=255)

                                                                             ccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGCTTGCGACAGCAGCCGCCAGa          >  2:1267415/1‑93 (MQ=255)

                                                                              cGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGAc         <  1:1518705/93‑1 (MQ=255)

                                                                               gCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACg        >  2:3756997/1‑93 (MQ=255)

                                                                               gCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACg        >  2:4088146/1‑93 (MQ=255)

                                                                                 cTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt      >  1:1554759/1‑93 (MQ=255)

                                                                                 cTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt      >  1:2554795/1‑93 (MQ=255)

                                                                                 cTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt      >  1:3245043/1‑93 (MQ=255)

                                                                                 cTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt      <  1:373191/93‑1 (MQ=255)

                                                                                 cTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt      <  1:1185448/93‑1 (MQ=255)

                                                                                  tGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTg     <  2:2555457/93‑1 (MQ=255)

                                                                                  tGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTg     >  1:2694368/1‑93 (MQ=255)

                                                                                   gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg    >  1:2420509/1‑93 (MQ=255)

                                                                                   gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg    >  2:3625926/1‑93 (MQ=255)

                                                                                   gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg    <  1:1156094/93‑1 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt  <  1:2298319/93‑1 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt  <  2:3540352/93‑1 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt  >  2:1293356/1‑93 (MQ=255)

TGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAT  >  NC_000913/2208415‑2208592

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