BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,236,190	C→T	G326G (GGC→GGT)	yjbH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,236,190	0	C	T	85.3%	842.1 / 97.8	333	G326G (GGC→GGT)	yjbH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (18/31); new base T (139/145); total (157/176)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-02

Only 100 of 334 total aligned reads displayed.

GATGCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTG  >  NC_000913/4236098‑4236280

gATGCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGc                                                                                            <  2:855267/93‑1 (MQ=255)

   gCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGAt                                                                                         >  2:3551785/1‑93 (MQ=255)

   gCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGAt                                                                                         <  1:3494056/93‑1 (MQ=255)

   gCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGAt                                                                                         <  1:5779178/93‑1 (MQ=255)

   gCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGAt                                                                                         >  1:1373745/1‑93 (MQ=255)

   gCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAAGCGAAAGGTGAt                                                                                         <  2:3733583/93‑1 (MQ=255)

     cATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATAc                                                                                       >  2:6799440/1‑93 (MQ=255)

       ttttGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGc                                                                                     >  1:1572455/1‑93 (MQ=255)

           gcagcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTa                                                                                 <  2:5823651/93‑1 (MQ=255)

           gcagcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTa                                                                                 <  1:6739388/93‑1 (MQ=255)

           gcagcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTa                                                                                 >  1:4977526/1‑93 (MQ=255)

            cagcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTAt                                                                                <  1:7477834/93‑1 (MQ=255)

             agcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATg                                                                               <  2:6762157/93‑1 (MQ=255)

             agcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATg                                                                               <  2:3843848/93‑1 (MQ=255)

              gcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGt                                                                              >  2:6373882/1‑93 (MQ=255)

              gcaTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGt                                                                              <  2:1190818/93‑1 (MQ=255)

                aTTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTa                                                                            >  2:3420014/1‑93 (MQ=255)

                 ttCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTAc                                                                           <  2:305236/93‑1 (MQ=255)

                  tCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTAcc                                                                          <  1:5326201/93‑1 (MQ=255)

                    ggtggtGGCGAATCAGTTAACGCTGGTGAAATACAATACCGGACTTGCCGATCCACAGATCCAGGTGAAAGGCGATACGCTGTATGTTACCgg                                                                        >  1:757813/1‑93 (MQ=255)

                     gtggtgGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGc                                                                       >  1:1788523/1‑93 (MQ=255)

                         tggCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGc                                                                   <  1:4640234/93‑1 (MQ=255)

                         nngCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGc                                                                   <  1:2894482/91‑1 (MQ=255)

                          ggCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCTGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCa                                                                  <  1:3590088/93‑1 (MQ=255)

                          ggCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCa                                                                  >  2:4566189/1‑93 (MQ=255)

                          ggCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCa                                                                  <  1:3182939/93‑1 (MQ=255)

                          ggCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCa                                                                  <  1:2618233/93‑1 (MQ=255)

                             gAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGt                                                               >  2:6294349/1‑93 (MQ=255)

                             gAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGt                                                               <  1:6021242/93‑1 (MQ=255)

                              aaTCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTg                                                              >  1:237681/1‑93 (MQ=255)

                               aTCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGa                                                             >  2:1878421/1‑93 (MQ=255)

                               aTCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGa                                                             >  1:2440503/1‑93 (MQ=255)

                                 cAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGaaa                                                           <  2:2217474/93‑1 (MQ=255)

                                  aGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAat                                                          >  2:1285578/1‑93 (MQ=255)

                                  aGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAat                                                          >  1:5449499/1‑93 (MQ=255)

                                   gTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAata                                                         >  2:2110650/1‑93 (MQ=255)

                                   gTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAata                                                         >  1:6579386/1‑93 (MQ=255)

                                   gTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAata                                                         <  2:4330894/93‑1 (MQ=255)

                                   gTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAata                                                         >  2:7463940/1‑93 (MQ=255)

                                     tAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATc                                                       <  2:3624841/93‑1 (MQ=255)

                                     tAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATc                                                       >  1:5621598/1‑93 (MQ=255)

                                       aCGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGt                                                     <  2:14971/93‑1 (MQ=255)

                                       aCGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGt                                                     <  2:33105/93‑1 (MQ=255)

                                         gCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGa                                                   <  2:2008187/93‑1 (MQ=255)

                                            gTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTc                                                <  1:5269289/93‑1 (MQ=255)

                                            gTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTc                                                >  2:2636693/1‑93 (MQ=255)

                                            gTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTc                                                <  1:7399371/93‑1 (MQ=255)

                                            gTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTc                                                <  2:3803961/93‑1 (MQ=255)

                                              tGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTcgc                                              >  1:889334/1‑93 (MQ=255)

                                              tGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTcgc                                              <  2:4121559/93‑1 (MQ=255)

                                              tGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTcgc                                              <  2:163151/93‑1 (MQ=255)

                                                aaaTACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTcgcgc                                            <  1:4878025/93‑1 (MQ=255)

                                                   tACAATGCCGGACTTTCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCTGCGAGCAGGTGAAATATCGTGATTCGCGCGaa                                         >  1:7059418/1‑93 (MQ=255)

                                                   tACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGaa                                         >  1:2020466/1‑93 (MQ=255)

                                                    aCAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAg                                        <  1:7144116/93‑1 (MQ=255)

                                                     cAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAgg                                       >  1:2524965/1‑93 (MQ=255)

                                                      aaTGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAggg                                      >  2:2469517/1‑93 (MQ=255)

                                                      aaTGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAggg                                      <  1:4152090/93‑1 (MQ=255)

                                                       gtGCCGGCCTAGCCGATCCACAGATCGAGGCGAAAGGTGATACGCTGTATGTTCCCGGCTAGCAGGTGAAATATCGTGATTCGCGCGAAGGGa                                     <  1:3845143/92‑1 (MQ=255)

                                                       aTGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGa                                     >  1:5971639/1‑93 (MQ=255)

                                                         gCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGatc                                   <  2:4785186/93‑1 (MQ=255)

                                                          ccGGACTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGatca                                  <  1:3889429/93‑1 (MQ=255)

                                                            ggACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGatcatc                                <  1:5354735/93‑1 (MQ=255)

                                                              aCTTGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCg                              <  1:274305/93‑1 (MQ=255)

                                                                ttGCCGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTg                            >  1:5723190/1‑93 (MQ=255)

                                                                  gCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCg                          <  2:5219754/93‑1 (MQ=255)

                                                                    cGATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGaa                        <  1:4451663/93‑1 (MQ=255)

                                                                     gATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAAt                       >  1:3618086/1‑93 (MQ=255)

                                                                     gATCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAAt                       >  1:3129151/1‑93 (MQ=255)

                                                                     gATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAAt                       >  1:3737638/1‑93 (MQ=255)

                                                                     gATCCACAGATCCAGGAGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAAt                       <  2:5598116/93‑1 (MQ=255)

                                                                       tCCACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCg                     <  1:2142902/93‑1 (MQ=255)

                                                                        ccACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCgg                    >  1:5174516/1‑93 (MQ=255)

                                                                        ccACAGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCgg                    >  2:4677373/1‑93 (MQ=255)

                                                                        ccACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCgg                    >  1:1262001/1‑93 (MQ=255)

                                                                         cacaGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGa                   <  1:6081200/93‑1 (MQ=255)

                                                                         cacaGATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGa                   <  1:4088181/93‑1 (MQ=255)

                                                                          acaGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGAt                  <  2:5131475/93‑1 (MQ=255)

                                                                             gATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGt               <  2:6040305/93‑1 (MQ=255)

                                                                             gATCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGt               >  2:1946054/1‑93 (MQ=255)

                                                                              aTCCAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtg              >  1:112456/1‑93 (MQ=255)

                                                                                ccAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtgat            <  1:2265850/93‑1 (MQ=255)

                                                                                ccAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtgat            <  2:4273653/93‑1 (MQ=255)

                                                                                 cAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtgatg           <  2:7051716/93‑1 (MQ=255)

                                                                                 cAGGCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtgatg           >  2:6196508/1‑93 (MQ=255)

                                                                                  aGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGtgatga          >  1:6896137/1‑93 (MQ=255)

                                                                                   ggCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGaa         >  1:4809691/1‑93 (MQ=255)

                                                                                   ggCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGaa         >  2:2179662/1‑93 (MQ=255)

                                                                                   ggCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGaa         >  2:7233507/1‑93 (MQ=255)

                                                                                   ggCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGaa         <  2:1725265/93‑1 (MQ=255)

                                                                                   ggCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGaa         <  1:3556600/93‑1 (MQ=255)

                                                                                    gCGAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAAc        <  2:1061964/93‑1 (MQ=255)

                                                                                      gAAAGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGa      <  1:1203264/93‑1 (MQ=255)

                                                                                        aaGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATc    <  2:5340394/93‑1 (MQ=255)

                                                                                        aaGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATc    >  2:2792508/1‑93 (MQ=255)

                                                                                        aaGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATc    <  2:3082472/93‑1 (MQ=255)

                                                                                         aGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCt   >  2:6594791/1‑93 (MQ=255)

                                                                                         aGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCt   <  1:825172/93‑1 (MQ=255)

                                                                                         aGGTGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCt   <  2:7334102/93‑1 (MQ=255)

                                                                                         aGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCt   >  1:6529636/1‑93 (MQ=255)

GATGCCATTTTGCAGCATTCGGTGGTGGCGAATCAGTTAACGCTGTTGAAATACAATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTG  >  NC_000913/4236098‑4236280

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