BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,971,922	T→C	D476G (GAC→GGC)	cheA ←	CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,971,922	0	T	C	90.4%	611.7 / 42.0	208	D476G (GAC→GGC)	cheA	CheA
Reads supporting (aligned to +/- strand): ref base T (12/8); new base C (94/94); total (106/102)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 213 total aligned reads displayed.

GGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTTCGTCGTCGC  >  NC_000913/1971832‑1972012

ggCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGccc                                                                                          <  1:5923229/93‑1 (MQ=255)

  cAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCAt                                                                                        <  2:1476274/93‑1 (MQ=255)

  cAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCAt                                                                                        >  2:2954541/1‑93 (MQ=255)

  cAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCAt                                                                                        <  1:577493/93‑1 (MQ=255)

   aGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATg                                                                                       >  1:6470316/1‑93 (MQ=255)

   aGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATg                                                                                       <  2:3457866/93‑1 (MQ=255)

   aGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATg                                                                                       <  2:3984222/93‑1 (MQ=255)

   aGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATg                                                                                       <  1:6068750/93‑1 (MQ=255)

     tAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGcc                                                                                     <  2:50579/93‑1 (MQ=255)

     tAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGcc                                                                                     >  1:5746438/1‑93 (MQ=255)

      aaaaTGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCg                                                                                    <  2:3576625/93‑1 (MQ=255)

       aaaTGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGa                                                                                   >  1:2530906/1‑93 (MQ=255)

       aaaTGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGa                                                                                   <  2:563730/93‑1 (MQ=255)

             ggATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGAcgccgc                                                                             <  1:6116600/93‑1 (MQ=255)

              gATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCcgcg                                                                            >  2:1574166/1‑93 (MQ=255)

                tCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGcc                                                                          <  2:2891721/93‑1 (MQ=255)

                  gTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCg                                                                        >  2:2873204/1‑93 (MQ=255)

                    aGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGga                                                                      <  1:4860492/93‑1 (MQ=255)

                     gTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGgag                                                                     >  2:6314013/1‑93 (MQ=255)

                     gTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGgag                                                                     >  2:5785887/1‑93 (MQ=255)

                        ccAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACg                                                                  >  2:6119986/1‑93 (MQ=255)

                         cAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGt                                                                 >  2:3621144/1‑93 (MQ=255)

                          aGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTc                                                                <  2:3062222/93‑1 (MQ=255)

                           gTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCg                                                               <  2:292603/93‑1 (MQ=255)

                            tACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCgg                                                              >  2:3346614/1‑93 (MQ=255)

                            tACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCgg                                                              <  2:6465477/93‑1 (MQ=255)

                             tcccTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGt                                                             <  1:3168773/92‑1 (MQ=255)

                             aCCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGt                                                             <  1:1452279/93‑1 (MQ=255)

                              cccTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTg                                                            >  2:5571138/1‑93 (MQ=255)

                               ccTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGa                                                           >  1:2138412/1‑93 (MQ=255)

                                cTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGAc                                                          >  1:5722461/1‑93 (MQ=255)

                                cTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCACCCGGAGACGTCGGTGAc                                                          >  1:1112244/1‑93 (MQ=255)

                                 tGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGAcc                                                         <  1:6331263/93‑1 (MQ=255)

                                 cgCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGAcc                                                         <  2:4196199/92‑1 (MQ=255)

                                  gCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCt                                                        >  1:1553662/1‑93 (MQ=255)

                                   cTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTg                                                       >  2:6588332/1‑93 (MQ=255)

                                   cTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTg                                                       >  2:1167688/1‑93 (MQ=255)

                                      cGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGctc                                                    <  1:2792161/93‑1 (MQ=255)

                                      cGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGctc                                                    <  1:679304/93‑1 (MQ=255)

                                       gACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGctct                                                   >  1:11540/1‑93 (MQ=255)

                                        aCTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTg                                                  >  2:6102772/1‑93 (MQ=255)

                                        aCTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTg                                                  >  1:1925454/1‑93 (MQ=255)

                                        aCTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTg                                                  <  1:3040726/93‑1 (MQ=255)

                                         cTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGc                                                 <  2:1883985/93‑1 (MQ=255)

                                         cTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGc                                                 <  2:810211/93‑1 (MQ=255)

                                            gATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGt                                              >  2:1339582/1‑93 (MQ=255)

                                            gATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGt                                              <  1:1767751/93‑1 (MQ=255)

                                             aTTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTg                                             <  1:881590/93‑1 (MQ=255)

                                             aTTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTg                                             <  1:6190379/93‑1 (MQ=255)

                                                tCGACATGAGCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGgag                                          >  2:874054/1‑93 (MQ=255)

                                                tCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGgag                                          >  1:1340363/1‑93 (MQ=255)

                                                 cGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGgaga                                         <  2:4846252/93‑1 (MQ=255)

                                                    cATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGGGAAGc                                      >  2:2985836/1‑93 (MQ=255)

                                                    cATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGc                                      >  2:4343848/1‑93 (MQ=255)

                                                    cATGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGc                                      <  1:3629399/93‑1 (MQ=255)

                                                     aTGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGcc                                     >  2:1264452/1‑93 (MQ=255)

                                                     aTGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGcc                                     <  2:4266960/93‑1 (MQ=255)

                                                      tGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCa                                    <  2:266723/93‑1 (MQ=255)

                                                      tGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCa                                    <  1:3310009/93‑1 (MQ=255)

                                                      tGACCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCa                                    >  2:5861750/1‑93 (MQ=255)

                                                        aCCGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAgg                                  <  1:4121896/93‑1 (MQ=255)

                                                         ccGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGt                                 <  1:51706/93‑1 (MQ=255)

                                                          cGCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTg                                >  1:4621759/1‑93 (MQ=255)

                                                           gCCCATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGc                               <  1:1639416/93‑1 (MQ=255)

                                                           ccccATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGc                               <  2:6246638/92‑1 (MQ=255)

                                                             ccATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCaa                             >  2:5707965/1‑93 (MQ=255)

                                                             ccATCTTCTGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCaa                             >  1:3022617/1‑93 (MQ=255)

                                                               aTCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAat                           >  1:3563247/1‑93 (MQ=255)

                                                                 cttctGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAatat                         <  2:1471168/93‑1 (MQ=255)

                                                                   tctGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCa                       <  2:3661960/93‑1 (MQ=255)

                                                                     tGGATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAga                     >  1:42087/1‑92 (MQ=255)

                                                                      ggATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCa                    <  2:2211179/93‑1 (MQ=255)

                                                                       gATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCAt                   <  2:1722876/93‑1 (MQ=255)

                                                                       gATATTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCAt                   <  1:5026173/93‑1 (MQ=255)

                                                                         tatTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATcg                 <  1:2619207/93‑1 (MQ=255)

                                                                         tatTACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATcg                 <  1:3382453/93‑1 (MQ=255)

                                                                          atTACGTTTAACGACGCCCATGCCTACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATcgc                <  1:4366120/93‑1 (MQ=255)

                                                                           ttACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATcgcg               >  2:2754523/1‑93 (MQ=255)

                                                                            tACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGa              >  1:3102797/1‑93 (MQ=255)

                                                                            tACGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGa              <  2:665074/93‑1 (MQ=255)

                                                                             aCGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGAc             <  1:2366118/93‑1 (MQ=255)

                                                                             aCGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGAc             >  1:5124259/1‑93 (MQ=255)

                                                                              cGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACt            <  1:4177152/93‑1 (MQ=255)

                                                                              cGTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACt            >  1:4737211/1‑93 (MQ=255)

                                                                               gTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCCGGTGCAAATATCAGCATCGCGACtt           <  1:2308294/93‑1 (MQ=255)

                                                                               gTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACtt           <  2:5172528/93‑1 (MQ=255)

                                                                               gTTTAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACtt           <  2:2369369/93‑1 (MQ=255)

                                                                                tttAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtc          >  2:3566216/1‑93 (MQ=255)

                                                                                 ttAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcg         >  1:2717214/1‑93 (MQ=255)

                                                                                  tAACGACGCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgt        <  1:4010477/93‑1 (MQ=255)

                                                                                     cgacgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgt     <  1:4511125/93‑1 (MQ=255)

                                                                                     cgacgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgt     >  2:360185/1‑93 (MQ=255)

                                                                                      gacgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtc    <  1:997594/93‑1 (MQ=255)

                                                                                      gacgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtc    <  1:1646534/93‑1 (MQ=255)

                                                                                       acgTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtcg   >  1:2975113/1‑93 (MQ=255)

                                                                                       acgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtcg   >  2:3561278/1‑93 (MQ=255)

                                                                                       acgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtcg   >  1:3737268/1‑93 (MQ=255)

                                                                                       acgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTtcgtcgtcg   >  2:947232/1‑93 (MQ=255)

                                                                                       acgCCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGAGCAAATATCAGCATCGCGACTtcgtcgtcg   >  2:4711341/1‑93 (MQ=255)

                                                                                        cgTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTTCGTCGTCGc  <  1:3933256/93‑1 (MQ=255)

GGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGATTTCGACATGACCGCCCATCTTCTGGATATTACGTTTAACGACGTCCATGCCGACGCCGCGCCCGGAGACGTCGGTGACCTGCTCTGCCGTGGAGAAGCCAGGTGCAAATATCAGCATCGCGACTTCGTCGTCGC  >  NC_000913/1971832‑1972012

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