BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,382,443	A→G	T357A (ACG→GCG)	yjeM →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,382,443	0	A	G	99.5%	1209.4 / ‑4.5	394	T357A (ACG→GCG)	yjeM	putative transport
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (209/181); total (210/184)
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.68e-01

Only 100 of 404 total aligned reads displayed.

GCGCGTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGACGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGGGTTGGTTAC  >  NC_000913/4382353‑4382534

 cgcgTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGgcgc                                                                                          <  2:1186879/93‑1 (MQ=255)

   cgTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCg                                                                                        <  2:7550688/93‑1 (MQ=255)

    gTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGACAGTCATCATCCA‑GGGGGCGCCGa                                                                                       >  2:774408/1‑93 (MQ=255)

    gTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGa                                                                                       >  2:2557407/1‑93 (MQ=255)

     tATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGaa                                                                                      >  2:5157367/1‑93 (MQ=255)

     tATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGaa                                                                                      >  2:6648069/1‑93 (MQ=255)

      aTTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGaaa                                                                                     >  2:5555181/1‑93 (MQ=255)

       ttACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAg                                                                                    >  2:6699515/1‑93 (MQ=255)

         aCTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCa                                                                                  <  2:6180541/93‑1 (MQ=255)

          cTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCAt                                                                                 <  2:4811321/93‑1 (MQ=255)

          cTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCAt                                                                                 >  1:3396198/1‑93 (MQ=255)

           tGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCAtt                                                                                >  2:2605149/1‑93 (MQ=255)

            ggACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATtg                                                                               >  2:4108202/1‑93 (MQ=255)

            ggACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATtg                                                                               <  2:2771774/93‑1 (MQ=255)

             gACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATtgt                                                                              <  1:6043724/93‑1 (MQ=255)

              aCTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATtgtg                                                                             <  1:7420679/93‑1 (MQ=255)

                tttCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGc                                                                           >  1:5082506/1‑93 (MQ=255)

                tttCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATCGTGGc                                                                           >  2:3849723/1‑93 (MQ=255)

                  tCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCg                                                                         >  2:3094548/1‑93 (MQ=255)

                    gATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGa                                                                       >  1:6720697/1‑93 (MQ=255)

                      tGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAAc                                                                     >  1:7121323/1‑93 (MQ=255)

                       gTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAAcc                                                                    >  2:738543/1‑93 (MQ=255)

                         tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGa                                                                  <  2:2493775/93‑1 (MQ=255)

                            tCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGa                                                               <  2:1814245/93‑1 (MQ=255)

                             cGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAc                                                              <  2:6099609/93‑1 (MQ=255)

                              gCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcg                                                             >  1:662041/1‑93 (MQ=255)

                               ccTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgc                                                            >  2:4030008/1‑93 (MQ=255)

                                cTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgcg                                                           >  1:4048384/1‑93 (MQ=255)

                                 tataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACgcgc                                                          <  2:4493061/93‑1 (MQ=255)

                                   taCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCt                                                        >  1:7671203/1‑93 (MQ=255)

                                   taCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCt                                                        >  1:955382/1‑93 (MQ=255)

                                   taCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCt                                                        >  1:2734711/1‑93 (MQ=255)

                                      cGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGaa                                                     <  2:1013703/93‑1 (MQ=255)

                                       ggTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAAt                                                    <  2:568707/93‑1 (MQ=255)

                                        gTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATg                                                   >  1:5664912/1‑93 (MQ=255)

                                        gTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATg                                                   >  2:4664196/1‑93 (MQ=255)

                                        gTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATg                                                   <  1:3220162/93‑1 (MQ=255)

                                         tGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGCACCGATGACGCGCCTGAATgg                                                  >  2:2218959/1‑92 (MQ=255)

                                         tGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGc                                                  >  1:5853689/1‑93 (MQ=255)

                                          gCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCg                                                 >  2:3261470/1‑93 (MQ=255)

                                           ctttCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGa                                                <  1:2919968/91‑1 (MQ=255)

                                           cGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGa                                                <  2:6934551/93‑1 (MQ=255)

                                            gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                               >  2:2364297/1‑93 (MQ=255)

                                            gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                               >  2:5376600/1‑93 (MQ=255)

                                            gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                               >  2:83261/1‑93 (MQ=255)

                                            gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                               <  1:830260/93‑1 (MQ=255)

                                            gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                               <  1:1989421/93‑1 (MQ=255)

                                             ttCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATg                                              <  1:6817808/93‑1 (MQ=255)

                                                tttACGCTTTGCTATTCACAGTTGAAAGCCATCATCCA‑GGGGCCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATgggg                                           <  2:638698/93‑1 (MQ=255)

                                                 ttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGa                                          >  1:1070279/1‑93 (MQ=255)

                                                 ttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGa                                          <  2:3618258/93‑1 (MQ=255)

                                                 ttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGa                                          >  1:4930559/1‑93 (MQ=255)

                                                   aCGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATg                                        >  1:3948712/1‑93 (MQ=255)

                                                   aCGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATg                                        >  2:3190864/1‑93 (MQ=255)

                                                    cGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGc                                       >  1:477384/1‑93 (MQ=255)

                                                      cTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCg                                     <  2:193167/93‑1 (MQ=255)

                                                      cTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCg                                     >  1:6232515/1‑93 (MQ=255)

                                                       tttGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCTCCTGAATGCGATGGGGATGCCGt                                    >  2:4364787/1‑93 (MQ=255)

                                                       tttGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGt                                    >  1:5440242/1‑93 (MQ=255)

                                                        ttGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTc                                   <  1:4324882/93‑1 (MQ=255)

                                                          gCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTa                                 >  1:5786490/1‑93 (MQ=255)

                                                          gCTATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTa                                 <  1:2302681/93‑1 (MQ=255)

                                                            tATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                               >  1:5907271/1‑93 (MQ=255)

                                                            tATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                               <  1:298797/93‑1 (MQ=255)

                                                            tATTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                               <  2:6537650/93‑1 (MQ=255)

                                                             aTTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCg                              <  1:7821823/93‑1 (MQ=255)

                                                             aTTCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCg                              >  2:1083754/1‑93 (MQ=255)

                                                              ttCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGc                             >  1:3706000/1‑93 (MQ=255)

                                                              ttCACCGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGc                             <  2:6821007/93‑1 (MQ=255)

                                                                  ccGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtg                         >  1:3791033/1‑93 (MQ=255)

                                                                  ccGTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtg                         <  2:2046072/93‑1 (MQ=255)

                                                                    gTTGAAAGCCATCATCCAGGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgt                        <  2:3707816/93‑1 (MQ=255)

                                                                    gTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgtg                       >  1:2792661/1‑93 (MQ=255)

                                                                    gTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgtg                       <  1:2575082/93‑1 (MQ=255)

                                                                    gTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgtg                       >  1:7276831/1‑93 (MQ=255)

                                                                    gTTGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgtg                       >  2:493606/1‑93 (MQ=255)

                                                                     ttGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGCTGCCGTCTATCGCCATGTgg                      >  1:556783/1‑93 (MQ=255)

                                                                     ttGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTgg                      >  2:5842373/1‑93 (MQ=255)

                                                                      tGAAAGCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGa                     <  2:1515342/93‑1 (MQ=255)

                                                                           gCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg                <  1:3822696/93‑1 (MQ=255)

                                                                           gCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg                <  2:1685637/93‑1 (MQ=255)

                                                                           gCCATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg                >  1:6938987/1‑93 (MQ=255)

                                                                            ccATCATCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGt               <  2:708281/93‑1 (MQ=255)

                                                                             catcatCCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTg              >  1:2200278/1‑93 (MQ=255)

                                                                              atcatcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCATTGc             >  2:3533951/1‑93 (MQ=255)

                                                                              atcatcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGc             <  1:5141517/93‑1 (MQ=255)

                                                                               tcatcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCg            >  2:290980/1‑93 (MQ=255)

                                                                               tcatcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCg            <  1:4522888/93‑1 (MQ=255)

                                                                                 atcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCggg          >  1:1421924/1‑93 (MQ=255)

                                                                                 atcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCggg          <  2:4234865/93‑1 (MQ=255)

                                                                                  tcCA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggt         >  1:1631041/1‑93 (MQ=255)

                                                                                    cA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg       <  1:7188240/93‑1 (MQ=255)

                                                                                    cA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg       >  2:2986298/1‑93 (MQ=255)

                                                                                    cA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg       >  1:24141/1‑93 (MQ=255)

                                                                                    cA‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg       >  2:7650805/1‑93 (MQ=255)

                                                                                     a‑GGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttgg      <  1:167544/93‑1 (MQ=255)

                                                                                       ggggGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggt     >  2:7280941/1‑93 (MQ=255)

                                                                                       ggggGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggt     <  2:5307239/93‑1 (MQ=255)

                                                                                        ggggCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggtt    >  1:2447523/1‑93 (MQ=255)

                                                                                        ggggCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggtt    >  2:4734934/1‑93 (MQ=255)

GCGCGTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCA‑GGGGACGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGGGTTGGTTAC  >  NC_000913/4382353‑4382534

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