BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,365,636	A→G	R59G (AGG→GGG)	yhcG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,365,636	0	A	G	98.9%	569.2 / ‑2.7	179	R59G (AGG→GGG)	yhcG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (79/98); total (80/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 190 total aligned reads displayed.

GAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGC  >  NC_000913/3365546‑3365725

gaAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCggg                                                                                         >  2:3254936/1‑93 (MQ=255)

gaAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCggg                                                                                         <  1:2812282/93‑1 (MQ=255)

gaAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCggg                                                                                         <  1:1917605/93‑1 (MQ=255)

  aaCGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGc                                                                                       <  2:668474/93‑1 (MQ=255)

  aaCGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGc                                                                                       >  2:657842/1‑93 (MQ=255)

     ggTACGTAGCGTTAGGGCGTTAATGACCGGGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCtac                                                                                    <  2:1990443/93‑3 (MQ=255)

     ggTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGc                                                                                    >  2:4569708/1‑93 (MQ=255)

     ggTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGc                                                                                    >  2:587760/1‑93 (MQ=255)

      gTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCg                                                                                   >  1:4717150/1‑93 (MQ=255)

      gTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCg                                                                                   >  1:1706107/1‑93 (MQ=255)

        aCGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTa                                                                                 <  1:1051213/93‑1 (MQ=255)

        aCGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTa                                                                                 >  1:848319/1‑93 (MQ=255)

         cGTAGCGTTACCGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTAt                                                                                <  1:2108921/93‑1 (MQ=255)

          gTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATg                                                                               >  2:3187911/1‑93 (MQ=255)

          gTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGACCGGAGCTGCGGATg                                                                               >  2:1190556/1‑93 (MQ=255)

               gTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTgcg                                                                          >  1:4586394/1‑93 (MQ=255)

               gTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTgcg                                                                          >  2:42963/1‑93 (MQ=255)

                ttAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTgcgc                                                                         <  2:2572268/93‑1 (MQ=255)

                     gcgTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTa                                                                    <  1:4514398/93‑1 (MQ=255)

                      cgTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTaa                                                                   <  1:2044827/93‑1 (MQ=255)

                       gTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAAt                                                                  >  2:382196/1‑93 (MQ=255)

                        ttAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATc                                                                 >  1:1147181/1‑93 (MQ=255)

                         tAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCa                                                                <  1:1317529/93‑1 (MQ=255)

                         tAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCa                                                                <  1:4490072/93‑1 (MQ=255)

                          aaTGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCaa                                                               <  1:2353285/93‑1 (MQ=255)

                           aTGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAg                                                              <  1:1808995/93‑1 (MQ=255)

                             gACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCg                                                            <  1:4453369/93‑1 (MQ=255)

                               ccGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGAc                                                          <  2:3006267/93‑1 (MQ=255)

                               ccGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGAc                                                          <  1:4235052/93‑1 (MQ=255)

                               ccGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGAc                                                          <  2:3617290/93‑1 (MQ=255)

                                cgcgACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACt                                                         >  2:3915169/1‑93 (MQ=255)

                                cgcgACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACt                                                         >  1:2727954/1‑93 (MQ=255)

                                 gcgACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTa                                                        <  2:3951157/93‑1 (MQ=255)

                                 gcgACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTa                                                        >  1:3165958/1‑93 (MQ=255)

                                 gcgACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTa                                                        >  1:3179073/1‑93 (MQ=255)

                                   gACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATc                                                      <  1:1864338/93‑1 (MQ=255)

                                    aCGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCa                                                     >  1:501523/1‑93 (MQ=255)

                                    aCGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCa                                                     >  2:3210841/1‑93 (MQ=255)

                                      ctACCCGGAAATTGGGCGGAGGATTGGTGAATTTTGAACAAGG‑GGCGAAGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCaaa                                                   <  1:3669708/92‑1 (MQ=255)

                                       tACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAg                                                  >  1:3897445/1‑93 (MQ=255)

                                       tACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAg                                                  <  1:4603964/93‑1 (MQ=255)

                                       tACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAg                                                  >  2:2418180/1‑93 (MQ=255)

                                        aCCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAgg                                                 >  2:1337286/1‑93 (MQ=255)

                                         ccAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGa                                                <  2:4322558/93‑1 (MQ=255)

                                           aGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAtt                                              >  2:804585/1‑93 (MQ=255)

                                           aGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAtt                                              <  2:101676/93‑1 (MQ=255)

                                            ggAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAttt                                             >  1:3514379/1‑93 (MQ=255)

                                            ggAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAttt                                             <  1:1480757/93‑1 (MQ=255)

                                             gAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTa                                            <  2:4669503/93‑1 (MQ=255)

                                               aaTTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATg                                          <  2:1120922/93‑1 (MQ=255)

                                               aaTTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATg                                          <  1:355136/93‑1 (MQ=255)

                                                 ttGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTc                                        >  2:2320917/1‑93 (MQ=255)

                                                 ttGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTc                                        <  2:2013008/93‑1 (MQ=255)

                                                  tGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCt                                       <  2:3759428/93‑1 (MQ=255)

                                                   ggCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTa                                      >  1:4439112/1‑93 (MQ=255)

                                                   ggCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTa                                      <  2:3369352/93‑1 (MQ=255)

                                                   ggCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTa                                      >  2:1962194/1‑93 (MQ=255)

                                                    gCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTaa                                     <  1:3370108/93‑1 (MQ=255)

                                                    gCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTaa                                     <  2:3539940/93‑1 (MQ=255)

                                                     ccGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAg                                    <  2:1562412/93‑1 (MQ=255)

                                                      cgacgaATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAgg                                   <  1:1053754/93‑1 (MQ=255)

                                                      cgacgaATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAgg                                   <  1:2534292/93‑1 (MQ=255)

                                                        acgtATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGta                                 >  2:4232480/1‑93 (MQ=255)

                                                        acgaATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGta                                 <  2:1839434/93‑1 (MQ=255)

                                                          gaATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGtata                               >  1:3317071/1‑93 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              >  2:902020/1‑93 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              <  2:1753396/93‑1 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              <  1:2826799/93‑1 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              <  1:2826800/93‑1 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              <  1:3053589/93‑1 (MQ=255)

                                                           aaTTGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATaa                              <  2:3485684/93‑1 (MQ=255)

                                                             ttGTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGc                            >  1:1973810/1‑93 (MQ=255)

                                                               gTCGAA‑TTTGAACTAGTTGGCGAGGCCGGGGCGGCGTATGGCGCGCAGCTAATCAAGCGACTATCAAAGGATTTATATCTAAGGTATAAGCgg                          >  2:3901956/1‑92 (MQ=255)

                                                               gTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGt                          <  2:704042/93‑1 (MQ=255)

                                                               gTCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGt                          <  2:4405743/93‑1 (MQ=255)

                                                                tCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTg                         >  1:4244404/1‑93 (MQ=255)

                                                                tCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTg                         <  2:4714701/93‑1 (MQ=255)

                                                                tCGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTg                         <  2:4633635/93‑1 (MQ=255)

                                                                 cGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTgg                        <  1:3762533/93‑1 (MQ=255)

                                                                 cGAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTgg                        <  1:3424309/93‑1 (MQ=255)

                                                                  gAA‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTggg                       >  2:554660/1‑93 (MQ=255)

                                                                   aa‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGt                      <  2:2050312/93‑1 (MQ=255)

                                                                   aa‑TTTGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGt                      >  1:2968806/1‑93 (MQ=255)

                                                                      tttGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTtc                    <  1:420262/93‑1 (MQ=255)

                                                                       ttGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTtct                   <  1:4694322/93‑1 (MQ=255)

                                                                        tGAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTtctc                  <  1:272085/93‑1 (MQ=255)

                                                                         gAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTctct                 >  1:965524/1‑93 (MQ=255)

                                                                         gAACAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTctct                 <  2:1432114/93‑1 (MQ=255)

                                                                           aCAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGc               >  1:2218323/1‑93 (MQ=255)

                                                                           aCAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGc               >  1:1759492/1‑93 (MQ=255)

                                                                            cAAGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCa              >  2:2465464/1‑93 (MQ=255)

                                                                             aaGGTGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCaa             >  1:2354402/1‑93 (MQ=255)

                                                                                 tGGCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCaaaaaa         >  1:2928817/1‑93 (MQ=255)

                                                                                   gCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACt       >  2:1717266/1‑93 (MQ=255)

                                                                                   gCGAGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACt       >  2:3472573/1‑93 (MQ=255)

                                                                                      aGGCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTAc    >  2:4448504/1‑93 (MQ=255)

                                                                                       ggCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACg   >  2:746642/1‑93 (MQ=255)

                                                                                       ggCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACg   >  2:2808225/1‑93 (MQ=255)

                                                                                        gCCGGGGCTGCTTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGc  <  2:2725658/93‑1 (MQ=255)

                                                                                        gCCGGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGc  >  2:1750496/1‑93 (MQ=255)

GAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAA‑TTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGC  >  NC_000913/3365546‑3365725

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