BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,530	T→C	P312P (CCT→CCC)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,530	0	T	C	99.4%	1015.4 / ‑3.6	325	P312P (CCT→CCC)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (158/165); total (158/167)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 330 total aligned reads displayed.

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

taAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCgg                                                                                         >  1:43771/1‑93 (MQ=255)

  aGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTc                                                                                       <  1:1825035/93‑1 (MQ=255)

   gATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCt                                                                                      <  2:3508559/93‑1 (MQ=255)

     tCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTgg                                                                                    >  2:3260042/1‑92 (MQ=255)

      cGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAg                                                                                   <  1:2516545/93‑1 (MQ=255)

       gATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGt                                                                                  >  1:5943866/1‑93 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTa                                                                                 <  1:4117087/93‑1 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTa                                                                                 >  1:3488449/1‑93 (MQ=255)

          aaGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               >  1:2630138/1‑93 (MQ=255)

           aGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCg                                                                              >  2:2234482/1‑93 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             <  1:2432324/93‑1 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             <  2:5338041/93‑1 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             <  1:2839581/93‑1 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             >  1:6289040/1‑93 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             <  2:3789400/93‑1 (MQ=255)

             tCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTg                                                                            <  1:1014710/93‑1 (MQ=255)

             tCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTg                                                                            >  1:4863161/1‑93 (MQ=255)

             gcATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTg                                                                            <  1:7163069/92‑1 (MQ=255)

              cATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTgg                                                                           <  2:2374470/93‑1 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          >  1:7206598/1‑93 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          <  1:6182741/93‑1 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          <  2:7350219/93‑1 (MQ=255)

                tttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCg                                                                         >  2:4454714/1‑93 (MQ=255)

                   gtgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGaa                                                                      <  2:5353903/93‑1 (MQ=255)

                   gtgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGaa                                                                      <  2:3878478/93‑1 (MQ=255)

                   gtgCAAAGGGCGATTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGaa                                                                      <  1:993899/93‑1 (MQ=255)

                    tgCGAATGGCGAGTGTCGACGAAATCACTGGTCATATTGTGGACAATATTGTCGAACGAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAg                                                                     <  1:5891685/93‑1 (MQ=255)

                    tgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAg                                                                     <  1:1331193/93‑1 (MQ=255)

                       aaaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGTTCTGAGTATCGTGGCGGAAgggg                                                                  >  2:4316404/1‑93 (MQ=255)

                       aaaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgggg                                                                  >  1:4564939/1‑93 (MQ=255)

                       aaaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgggg                                                                  >  1:5761216/1‑93 (MQ=255)

                       aaaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgggg                                                                  >  2:6996710/1‑93 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgggga                                                                 >  1:159938/1‑92 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGt                                                                 <  1:7475264/93‑1 (MQ=255)

                           ggCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGa                                                              <  2:6893719/93‑1 (MQ=255)

                           ggCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGa                                                              <  2:1269437/93‑1 (MQ=255)

                            gCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGaa                                                             >  2:3413059/1‑93 (MQ=255)

                             cGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGaaa                                                            >  1:3717632/1‑93 (MQ=255)

                              gAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAc                                                           >  2:7790269/1‑93 (MQ=255)

                               aGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAcc                                                          >  2:1617161/1‑93 (MQ=255)

                                gtgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         <  2:6095553/93‑1 (MQ=255)

                                atgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         <  1:5376394/92‑1 (MQ=255)

                                  gtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAg                                                       <  2:4698255/93‑1 (MQ=255)

                                   tCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAgg                                                      >  2:3628896/1‑93 (MQ=255)

                                    cgacgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGa                                                     <  2:4593938/93‑1 (MQ=255)

                                       cgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCa                                                  <  1:117637/93‑1 (MQ=255)

                                        gaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAg                                                 >  1:7431707/1‑93 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAgg                                                <  2:7930106/93‑1 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAgg                                                >  1:1209448/1‑93 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAgg                                                <  1:7713565/93‑1 (MQ=255)

                                          aaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGc                                               >  1:5200459/1‑93 (MQ=255)

                                          aaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGc                                               <  1:401217/93‑1 (MQ=255)

                                          aaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGc                                               >  1:726246/1‑93 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             >  2:6495329/1‑93 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             >  2:5193316/1‑93 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           <  1:1691852/93‑1 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           >  2:2537057/1‑93 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGCAGAAACCCACGAGGAGGCGGAt                                           >  2:7696540/1‑93 (MQ=255)

                                                  gTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTaa                                       <  2:1925203/93‑1 (MQ=255)

                                                  gTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTaa                                       <  1:6710506/93‑1 (MQ=255)

                                                  gTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTaa                                       >  1:4754163/1‑93 (MQ=255)

                                                   tCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAAt                                      >  1:1889508/1‑93 (MQ=255)

                                                   tCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAACCAGGAGCAGGCGGATTTAAt                                      >  2:4934033/1‑93 (MQ=255)

                                                     atatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGa                                    <  1:7909913/93‑1 (MQ=255)

                                                     atatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGa                                    <  2:6365891/93‑1 (MQ=255)

                                                        ttGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCg                                 >  1:6230499/1‑93 (MQ=255)

                                                        ttGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCg                                 >  1:2191905/1‑93 (MQ=255)

                                                        ttGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCg                                 >  1:5182208/1‑93 (MQ=255)

                                                         tgtgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCgg                                <  2:1502724/93‑1 (MQ=255)

                                                           tgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTa                              >  2:5540280/1‑93 (MQ=255)

                                                           tgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTa                              >  2:4665118/1‑93 (MQ=255)

                                                            ggACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaa                             >  1:1570126/1‑93 (MQ=255)

                                                             gACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaaa                            <  2:7643851/93‑1 (MQ=255)

                                                             gACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaaa                            <  2:1779590/93‑1 (MQ=255)

                                                             gACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaaa                            <  2:7437159/93‑1 (MQ=255)

                                                              aCAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAg                           <  1:5271195/93‑1 (MQ=255)

                                                               cAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATTGTGGCGGAAGGGGGAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAgg                          >  2:5992904/1‑93 (MQ=255)

                                                               cAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAgg                          >  2:3731911/1‑93 (MQ=255)

                                                                   atTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGtt                      >  2:1839688/1‑93 (MQ=255)

                                                                    ttGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTc                     >  2:139599/1‑93 (MQ=255)

                                                                     tGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCa                    <  1:7669222/93‑1 (MQ=255)

                                                                         gAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACttt                >  2:3053846/1‑93 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtttt               >  2:6770055/1‑93 (MQ=255)

                                                                           aCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACttttt              <  1:6394378/93‑1 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGTGTTCACTTTTTGc            >  1:7046293/1‑93 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGc            <  2:2489376/93‑1 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGc            <  2:2778353/93‑1 (MQ=255)

                                                                              aGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCa           <  1:836749/93‑1 (MQ=255)

                                                                              aGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCa           <  1:3074696/93‑1 (MQ=255)

                                                                                cgcgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAgg         >  1:5537109/1‑93 (MQ=255)

                                                                                 gcgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAggg        <  2:2925064/93‑1 (MQ=255)

                                                                                  cgTAATCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc       <  2:6459814/93‑1 (MQ=255)

                                                                                  cgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc       >  2:7125589/1‑93 (MQ=255)

                                                                                  cgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc       >  2:7221266/1‑93 (MQ=255)

                                                                                   gTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCt      <  1:5334108/93‑1 (MQ=255)

                                                                                   gTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCt      >  1:6485592/1‑93 (MQ=255)

                                                                                    tAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTa     <  1:3346240/93‑1 (MQ=255)

                                                                                     aaGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAt    >  2:1200810/1‑93 (MQ=255)

                                                                                     aaGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAt    <  1:4237259/93‑1 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  >  1:7581890/1‑93 (MQ=255)

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

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