BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,338,481	C→T	A189V (GCC→GTC)	yrbG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,338,481	0	C	T	99.5%	803.1 / NA	219	A189V (GCC→GTC)	yrbG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (65/153); minor base G (1/0); total (66/153)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

CCAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGCCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATCG  >  NC_000913/3338391‑3338570

ccAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCa                                                                                         <  2:1201731/93‑1 (MQ=255)

 cAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAc                                                                                        <  2:6288717/93‑1 (MQ=255)

 cAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAc                                                                                        <  2:1165074/93‑1 (MQ=255)

 cAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAc                                                                                        >  1:5070886/1‑93 (MQ=255)

  agagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCCTTGCGCTTTTCATCATGCCCGTGGTCAcc                                                                                       >  1:2091244/1‑92 (MQ=255)

  agagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAcg                                                                                       <  1:5430271/93‑1 (MQ=255)

   gagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAcgc                                                                                      <  1:1445543/93‑1 (MQ=255)

   gagagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCAcgc                                                                                      >  2:3240742/1‑93 (MQ=255)

     gagCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCgg                                                                                    <  1:6025555/93‑1 (MQ=255)

      agcagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCCTTTCGCTTATTATCATTCCCGTGGTGAcgccgg                                                                                   >  2:1093241/1‑90 (MQ=255)

       gcagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCgggt                                                                                  >  2:4972622/1‑91 (MQ=255)

       gcagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAt                                                                                  <  2:6064369/93‑1 (MQ=255)

        cagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATTCTGGCAGGGGGCACGCGGGtg                                                                                 >  1:6156239/1‑93 (MQ=255)

        cagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtg                                                                                 <  2:1748771/93‑1 (MQ=255)

        cagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtg                                                                                 <  2:3331970/93‑1 (MQ=255)

        cagcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtg                                                                                 <  1:6069206/93‑1 (MQ=255)

         agcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtgg                                                                                <  1:4229111/93‑1 (MQ=255)

         agcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtgg                                                                                <  1:3458353/93‑1 (MQ=255)

          gcTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGGCTTGCGCTTATTATCATGACAGTGGTCCCGCGGAgggg                                                                               >  2:1544221/1‑92 (MQ=255)

            ttGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGAtggtgg                                                                             <  2:1007156/93‑1 (MQ=255)

              gCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGtt                                                                           <  1:2731173/93‑1 (MQ=255)

               cAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTg                                                                          <  1:1649127/93‑1 (MQ=255)

                agagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGa                                                                         >  2:3904725/1‑93 (MQ=255)

                 gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGGATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGAt                                                                        >  1:4943807/1‑93 (MQ=255)

                 gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGAt                                                                        <  1:752124/93‑1 (MQ=255)

                 gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGAt                                                                        >  2:4258837/1‑93 (MQ=255)

                 gagCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCCGCCTTTCGCTTATCATCATGCCCGTGGTCACGCCGATGGTGGGTGAt                                                                        >  2:2361867/1‑93 (MQ=255)

                   gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATaa                                                                      <  2:298991/93‑1 (MQ=255)

                   gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATaa                                                                      <  1:151042/93‑1 (MQ=255)

                   gCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCCGCATTTCGCTTATCATCATGCCACTGGTCACGCGGGTGGTGGGTGATaa                                                                      >  2:5492925/1‑93 (MQ=255)

                      gCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGc                                                                   <  2:4120882/93‑1 (MQ=255)

                       ccccgTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGcc                                                                  <  2:5857433/90‑1 (MQ=255)

                        cgcgTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCa                                                                 >  1:3993808/1‑93 (MQ=255)

                         gcgTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCAc                                                                <  2:5402389/93‑1 (MQ=255)

                           gTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGGATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGGGGTTGATAAAGCCACgg                                                              >  2:3418956/1‑93 (MQ=255)

                              aCGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGc                                                           <  2:4608867/93‑1 (MQ=255)

                               cggcggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCt                                                          <  1:2659058/93‑1 (MQ=255)

                                ggcggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTg                                                         <  1:1709753/93‑1 (MQ=255)

                                ggcggATTGCCCGTCGCGTTTTTATGGCTCCGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAACGCCACGGTGCTg                                                         >  2:3557761/1‑93 (MQ=255)

                                 gcggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTgg                                                        >  2:2861447/1‑93 (MQ=255)

                                   ggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATTATCATGCCCGTGGTCACGCGGGTGGTGGGTGATAAAGCCACCGTGCTcgcg                                                      >  2:6003727/1‑93 (MQ=255)

                                   ggATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCg                                                      <  2:1471712/93‑1 (MQ=255)

                                    gATTGCCCGTCGCGTTTTTATGGCTCGGCATTGGGCTTTTTATCATTCCAGTGGTCACCCGGCTGGTGGTTGATAACGTCACCGTGCTGGCGa                                                     >  2:4761807/1‑93 (MQ=255)

                                    gATTGCCCGTCGCGTTTTTATGGCTCCGCATTGCGCTTTTCATCATATCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGa                                                     >  2:5814279/1‑93 (MQ=255)

                                     aTTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGaa                                                    <  2:4500277/93‑1 (MQ=255)

                                     aTTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGaa                                                    <  1:3231716/93‑1 (MQ=255)

                                     aTTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGaa                                                    >  2:1887194/1‑93 (MQ=255)

                                      ttGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAAt                                                   <  2:6204737/93‑1 (MQ=255)

                                      ttGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAAt                                                   >  1:4679484/1‑93 (MQ=255)

                                       tGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAAtt                                                  <  1:1661919/93‑1 (MQ=255)

                                       tGCCCGTCGCGTTTTTATGGCTCCGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACCGTGCTCGCGAAtt                                                  >  1:2057010/1‑93 (MQ=255)

                                        gCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTa                                                 <  1:5343250/93‑1 (MQ=255)

                                        gCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTa                                                 <  1:1360767/93‑1 (MQ=255)

                                          ccGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTTGCGAATTACt                                               >  1:3482779/1‑93 (MQ=255)

                                           cGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACtt                                              <  1:3943441/93‑1 (MQ=255)

                                           cGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACtt                                              <  1:3266159/93‑1 (MQ=255)

                                            gTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACttt                                             <  1:1863315/93‑1 (MQ=255)

                                            gTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACttt                                             >  2:5215256/1‑93 (MQ=255)

                                            gTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGGTGATAACGCCACGGTGCTGGCGAATTACttt                                             >  1:287146/1‑93 (MQ=255)

                                             tCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTg                                            <  2:2177322/93‑1 (MQ=255)

                                               gcgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTggc                                          >  2:1520889/1‑91 (MQ=255)

                                                cgTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCa                                         <  1:1745233/93‑1 (MQ=255)

                                                 gTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCAt                                        <  2:2554731/93‑1 (MQ=255)

                                                  tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                       <  2:1389303/93‑1 (MQ=255)

                                                  tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                       >  2:1630600/1‑93 (MQ=255)

                                                  tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                       <  2:1939208/93‑1 (MQ=255)

                                                  tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                       >  1:5672690/1‑93 (MQ=255)

                                                  tttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATc                                       >  2:781559/1‑93 (MQ=255)

                                                   ttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCa                                      >  2:1254189/1‑93 (MQ=255)

                                                   ttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCa                                      <  2:352904/93‑1 (MQ=255)

                                                   ttttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCa                                      <  1:4750109/93‑1 (MQ=255)

                                                     ttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGc                                    >  2:4316407/1‑93 (MQ=255)

                                                     ttATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGc                                    <  1:2012049/93‑1 (MQ=255)

                                                      tATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCg                                   <  1:453806/93‑1 (MQ=255)

                                                       aTGGCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGa                                  <  2:5614864/93‑1 (MQ=255)

                                                          gCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGtt                               <  1:2283369/93‑1 (MQ=255)

                                                          gCTCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGtt                               <  1:3678367/93‑1 (MQ=255)

                                                            tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                             <  2:813631/93‑1 (MQ=255)

                                                            tCGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGa                             <  1:3293249/93‑1 (MQ=255)

                                                             cGGCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGAc                            <  1:105316/93‑1 (MQ=255)

                                                               gCATTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGa                          >  1:4313750/1‑93 (MQ=255)

                                                                 aTTGCGCTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATg                        <  2:4022468/93‑1 (MQ=255)

                                                                    gcgcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGCCGATGGGt                     >  1:1246590/1‑93 (MQ=255)

                                                                     cgcTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTc                    >  2:2460475/1‑93 (MQ=255)

                                                                       cTTATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTg                  <  2:755900/93‑1 (MQ=255)

                                                                        ttATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGa                 <  1:25052/93‑1 (MQ=255)

                                                                         tATCATCATGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGAc                >  1:381897/1‑93 (MQ=255)

                                                                              tcatGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCaa           >  2:1328374/1‑93 (MQ=255)

                                                                                atGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAAtt         >  2:3823480/1‑93 (MQ=255)

                                                                                atGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAAtt         <  2:2721821/93‑1 (MQ=255)

                                                                                 tGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTg        >  2:3703644/1‑93 (MQ=255)

                                                                                 tGCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTg        <  1:5336776/93‑1 (MQ=255)

                                                                                  gCCAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCGGCGAGTTGACGATGGGTCTGACGGCAATTGc       <  1:30769/93‑1 (MQ=255)

                                                                                   ccAGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCt      >  2:5233605/1‑93 (MQ=255)

                                                                                     aGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTAt    >  1:4114094/1‑93 (MQ=255)

                                                                                     aGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTAt    >  1:109794/1‑93 (MQ=255)

                                                                                     aGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTAt    >  1:903829/1‑93 (MQ=255)

                                                                                     aGTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTCTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTAt    <  1:4428789/93‑1 (MQ=255)

                                                                                      gTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATc   <  1:234636/93‑1 (MQ=255)

                                                                                      gTGGTCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATc   >  1:4795608/1‑93 (MQ=255)

CCAGAGAGCAGCTTGCAGAGCTGCCGCGTGACGGCGGATTGCCCGTCGCGTTTTTATGGCTCGGCATTGCGCTTATCATCATGCCAGTGGCCACGCGGATGGTGGTTGATAACGCCACGGTGCTGGCGAATTACTTTGCCATCAGCGAGTTGACGATGGGTCTGACGGCAATTGCTATCG  >  NC_000913/3338391‑3338570

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