BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	804,639	C→T	G646G (GGC→GGT)	ybhJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	804,639	0	C	T	99.3%	497.1 / ‑5.8	144	G646G (GGC→GGT)	ybhJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (72/71); total (72/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 148 total aligned reads displayed.

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

gtTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGa                                                                                         >  2:3462817/1‑93 (MQ=255)

 tttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGAt                                                                                        <  2:565955/93‑1 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       <  2:3558573/93‑1 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       <  1:190491/93‑1 (MQ=255)

   tGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATgg                                                                                      <  2:3551458/93‑1 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   >  1:2363970/1‑93 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   >  1:230327/1‑93 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  1:2239250/93‑1 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  1:3699991/1‑93 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 <  2:2630068/93‑1 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 <  1:3987320/93‑1 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  2:3899360/1‑93 (MQ=255)

          aTTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgcg                                                                               >  2:3798679/1‑93 (MQ=255)

           ttAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgcgc                                                                              >  2:643342/1‑93 (MQ=255)

            tAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGcgcg                                                                             <  1:263100/93‑1 (MQ=255)

             aaGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGTGGTGAAACCAGGTGATGGTTCTGCGCGt                                                                            <  2:3958214/93‑1 (MQ=255)

               gCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGa                                                                          >  2:3058822/1‑93 (MQ=255)

               acAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGa                                                                          <  2:1206287/92‑1 (MQ=255)

                cAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGaa                                                                         <  2:4409482/93‑1 (MQ=255)

                  gATTGCTGGTCAGCAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  2:776198/1‑93 (MQ=255)

                     tGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggc                                                                    <  2:3368419/93‑1 (MQ=255)

                      gCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcg                                                                   >  2:3948340/1‑93 (MQ=255)

                       cTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcgg                                                                  >  1:778974/1‑93 (MQ=255)

                       cTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcgg                                                                  >  1:315853/1‑93 (MQ=255)

                        tGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcggc                                                                 <  2:3899606/93‑1 (MQ=255)

                        tGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcggc                                                                 <  1:3288784/93‑1 (MQ=255)

                        tGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAATGGCAGCATGGTCTCTGCGGTGAAACCAGGTGATGGTTCTGTGCGTGAACAggctgg                                                                 >  2:3464657/1‑90 (MQ=255)

                         ggTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGgcggcg                                                                <  2:3651667/93‑1 (MQ=255)

                          gTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGa                                                               <  2:1350508/93‑1 (MQ=255)

                             aGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                                                            >  1:2133265/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  1:3187678/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  2:3686646/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  2:3523773/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  1:1950716/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  2:204305/1‑93 (MQ=255)

                                   atatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGc                                                      >  1:87481/1‑93 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGct                                                     >  2:2598236/1‑92 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     <  1:1855609/93‑1 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     >  1:4739562/1‑93 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     >  1:4739563/1‑93 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    <  2:3143286/93‑1 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    >  2:3297735/1‑93 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    >  2:568446/1‑93 (MQ=255)

                                     atTGATACGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    <  1:4670596/93‑1 (MQ=255)

                                      ttGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg                                                   >  2:1988088/1‑93 (MQ=255)

                                      ttGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg                                                   >  1:4410157/1‑93 (MQ=255)

                                       tGATGCACTGCAAACTGAAATTGACAGCGGGGTCTATGCGGTGAAGCCAGGTGATGGTTCTGCGAGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  <  2:2598470/93‑1 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  >  1:3232501/1‑93 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  <  2:1537125/93‑1 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgg                                                  >  2:4103574/1‑92 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:1485255/93‑1 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:2976130/93‑1 (MQ=255)

                                           ccGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt                                              <  2:1704042/93‑1 (MQ=255)

                                            cGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg                                             <  2:4097398/93‑1 (MQ=255)

                                             gCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg                                            >  1:736494/1‑93 (MQ=255)

                                              cTGCAAACTGAAATTGGCAGCTTGGTCTATGCGGTGAACCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc                                           <  1:2302599/93‑1 (MQ=255)

                                              cTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc                                           >  2:4643565/1‑93 (MQ=255)

                                              cTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc                                           <  2:2963752/93‑1 (MQ=255)

                                              cTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc                                           >  2:4254963/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         >  2:2485446/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         <  1:4424261/93‑1 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt                                      <  1:269672/93‑1 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt                                      <  1:1181301/93‑1 (MQ=255)

                                                     cTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTgg                                    >  1:2503587/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   >  1:2770885/1‑93 (MQ=255)

                                                       gAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGATTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCg                                  <  1:4122662/93‑1 (MQ=255)

                                                         aaTTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGaa                                >  1:1193302/1‑93 (MQ=255)

                                                         aaTTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGaa                                <  1:3022291/93‑1 (MQ=255)

                                                           ttGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAata                              >  2:4536709/1‑93 (MQ=255)

                                                           ttGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAata                              >  1:1260297/1‑93 (MQ=255)

                                                            tGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAatat                             <  1:1573739/93‑1 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            >  1:4122848/1‑93 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            >  1:2513894/1‑93 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            <  2:3715817/93‑1 (MQ=255)

                                                               cagcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGc                          <  2:1705467/93‑1 (MQ=255)

                                                               cagcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGc                          >  1:398632/1‑93 (MQ=255)

                                                                agcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGcc                         >  2:2705188/1‑93 (MQ=255)

                                                                agcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAGTGcc                         >  2:4197383/1‑93 (MQ=255)

                                                                   aTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgag                      >  2:698290/1‑93 (MQ=255)

                                                                   aTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgag                      <  2:2448988/93‑1 (MQ=255)

                                                                   aTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgag                      <  1:2354146/93‑1 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  1:4285046/1‑93 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  2:2349787/1‑93 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    >  2:779584/1‑93 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  2:3649787/93‑1 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  2:234032/93‑1 (MQ=255)

                                                                      gTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgaggag                   >  2:1626896/1‑93 (MQ=255)

                                                                        cTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTa                 <  2:3108606/93‑1 (MQ=255)

                                                                        cTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTa                 <  2:156115/93‑1 (MQ=255)

                                                                        cTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTa                 >  1:104966/1‑93 (MQ=255)

                                                                         tATGCGGTGAAACCAGGTGATGGTTCTGCGCGTAAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAc                <  2:4431016/93‑1 (MQ=255)

                                                                           tGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcgc              <  1:523164/93‑1 (MQ=255)

                                                                            gCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcgcg             >  1:3848864/1‑93 (MQ=255)

                                                                              ggTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGAc           >  2:3710369/1‑93 (MQ=255)

                                                                               gTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcgcgcca          >  1:266400/1‑90 (MQ=255)

                                                                                  aaaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTaaa       <  1:4609156/93‑1 (MQ=255)

                                                                                  aaaCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTaaa       <  1:1090621/93‑1 (MQ=255)

                                                                                    aCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACg     >  1:4364815/1‑93 (MQ=255)

                                                                                      cAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCt   >  1:3066973/1‑93 (MQ=255)

                                                                                       aGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTa  >  1:2153591/1‑93 (MQ=255)

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

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