BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	657,814	G→A	A113T (GCC→ACC)	ybeM →	putative amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,814	0	G	A	99.5%	758.7 / NA	216	A113T (GCC→ACC)	ybeM	putative amidase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (88/127); minor base C (0/1); total (88/128)
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.99e-01

Only 100 of 231 total aligned reads displayed.

GCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCC  >  NC_000913/657724‑657903

gcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTAcc                                                                                         >  1:3114814/1‑93 (MQ=255)

  acaGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGc                                                                                       >  1:3506443/1‑93 (MQ=255)

  acaGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGc                                                                                       <  2:2815220/93‑1 (MQ=255)

    aGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTc                                                                                     >  1:25933/1‑93 (MQ=255)

     ggCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCg                                                                                    >  2:4403872/1‑93 (MQ=255)

       cATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTg                                                                                  <  2:3669373/93‑1 (MQ=255)

        aTTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgc                                                                                 >  1:2309333/1‑93 (MQ=255)

         ttACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcg                                                                                >  2:5250172/1‑93 (MQ=255)

         ttACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcg                                                                                >  1:807047/1‑93 (MQ=255)

           aCAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCt                                                                              <  2:5437632/93‑1 (MQ=255)

           aCAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCt                                                                              <  2:5683195/93‑1 (MQ=255)

           aCAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCt                                                                              <  2:5870585/93‑1 (MQ=255)

            cAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTg                                                                             <  1:6078880/93‑1 (MQ=255)

             aGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTgg                                                                            >  2:4529066/1‑93 (MQ=255)

              ttgAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           <  1:1723315/91‑1 (MQ=255)

              gggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           <  1:4648159/93‑1 (MQ=255)

              gggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           <  1:2800538/93‑1 (MQ=255)

               ggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGAt                                                                          >  1:1072277/1‑93 (MQ=255)

                   cTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATacca                                                                      <  2:4281226/93‑1 (MQ=255)

                     gAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATaccacc                                                                    <  2:390973/93‑1 (MQ=255)

                       aaTCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTg                                                                  >  1:4516566/1‑93 (MQ=255)

                           cTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTtat                                                              <  1:6140149/93‑1 (MQ=255)

                           cTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTtat                                                              >  1:5286477/1‑93 (MQ=255)

                             ggTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTatat                                                            <  2:4130616/93‑1 (MQ=255)

                                tcTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  1:2305293/92‑1 (MQ=255)

                                aCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  1:5885810/93‑1 (MQ=255)

                                aCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         >  2:866505/1‑93 (MQ=255)

                                aCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  2:5199646/93‑1 (MQ=255)

                                aCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  1:3980156/93‑1 (MQ=255)

                                 cTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTCCCGCTCGTGCGCTGGATACCACCTGTTATAtggtg                                                        <  1:278649/93‑1 (MQ=255)

                                 cTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggtg                                                        <  2:561858/93‑1 (MQ=255)

                                 cTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggtg                                                        <  2:500704/93‑1 (MQ=255)

                                   tCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggc                                                      <  1:5576993/93‑1 (MQ=255)

                                     cTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcgg                                                    <  1:2121756/93‑1 (MQ=255)

                                      tGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcggc                                                   <  1:4752440/93‑1 (MQ=255)

                                        ccgccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGcggcgg                                                 <  2:1376196/93‑1 (MQ=255)

                                         cgccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggg                                                <  2:898330/93‑1 (MQ=255)

                                          gccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggg                                               >  1:2845148/1‑93 (MQ=255)

                                          gccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggg                                               <  1:754354/93‑1 (MQ=255)

                                           ccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                              <  2:3607911/93‑1 (MQ=255)

                                            cTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGa                                             <  2:350472/93‑1 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGGTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            >  1:4441004/1‑93 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            >  2:3932897/1‑93 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            >  1:1051540/1‑93 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            <  1:888789/93‑1 (MQ=255)

                                              gggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATTGTGGCGGCGGGGGAGt                                           >  2:767193/1‑93 (MQ=255)

                                               ggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTg                                          >  1:2239351/1‑93 (MQ=255)

                                               ggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTg                                          >  1:255521/1‑93 (MQ=255)

                                                gTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGc                                         <  2:3799858/93‑1 (MQ=255)

                                                 ttCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCg                                        <  2:1359982/93‑1 (MQ=255)

                                                 ttCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCg                                        <  1:1287790/93‑1 (MQ=255)

                                                  tCGCGAGCGGCTCAAAGAGCATCGCTGGTCAACGTTACTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCgg                                       <  2:6305253/93‑1 (MQ=255)

                                                    gcgGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGa                                     <  1:3256225/93‑1 (MQ=255)

                                                    gcgGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGa                                     >  2:376432/1‑93 (MQ=255)

                                                    gcgGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGa                                     <  1:4868915/93‑1 (MQ=255)

                                                     cgGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGaa                                    <  2:477218/93‑1 (MQ=255)

                                                      gggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGCGCGGGAAc                                   >  2:491979/1‑93 (MQ=255)

                                                       ggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGGGCGGGAACa                                  >  1:5491094/1‑93 (MQ=255)

                                                         ccGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACaaa                                <  2:4363104/93‑1 (MQ=255)

                                                          cGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGCGGCGGCGGGGGACTGCGGCAACaaaa                               >  2:2154027/1‑93 (MQ=255)

                                                           gCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACaaaaa                              <  1:2371954/93‑1 (MQ=255)

                                                             tCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAata                            >  1:2428616/1‑93 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           <  2:327474/93‑1 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           >  1:1852665/1‑93 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           <  1:1505668/93‑1 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           <  2:5936710/93‑1 (MQ=255)

                                                               aaaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATc                          <  2:876840/93‑1 (MQ=255)

                                                               aaaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATc                          >  1:5766822/1‑93 (MQ=255)

                                                                aaGAGCTTCGCTGGTCATGGGTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGTCAGCGGGGGAGTGCGGGAACAAAAATATCg                         <  2:1269322/93‑1 (MQ=255)

                                                                aaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCg                         <  2:2933915/93‑1 (MQ=255)

                                                                aaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCg                         >  1:3543979/1‑93 (MQ=255)

                                                                aaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCg                         >  2:2187549/1‑93 (MQ=255)

                                                                 agagCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCgg                        <  2:2541028/93‑1 (MQ=255)

                                                                  gagCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGt                       >  1:1470495/1‑93 (MQ=255)

                                                                  gagCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGt                       >  2:2053224/1‑93 (MQ=255)

                                                                   agCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTc                      >  2:3150965/1‑93 (MQ=255)

                                                                    gCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCa                     <  1:3074296/93‑1 (MQ=255)

                                                                     cATCCCTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaa                    <  2:2212248/93‑1 (MQ=255)

                                                                     cATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaa                    <  1:568009/93‑1 (MQ=255)

                                                                      aTCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaaa                   <  1:5085228/93‑1 (MQ=255)

                                                                      aTCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaaa                   >  2:3209364/1‑93 (MQ=255)

                                                                       tCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAg                  <  2:6028072/93‑1 (MQ=255)

                                                                        cACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGc                 >  1:560725/1‑93 (MQ=255)

                                                                         aCTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGcc                >  2:5535133/1‑93 (MQ=255)

                                                                          cTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCg               >  2:754265/1‑93 (MQ=255)

                                                                           tGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCgg              >  1:4327226/1‑93 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             <  1:4998013/93‑1 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             <  2:1785845/93‑1 (MQ=255)

                                                                             gTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAt            >  1:6223527/1‑93 (MQ=255)

                                                                               cAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATta          <  1:3270129/93‑1 (MQ=255)

                                                                               cAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATta          <  2:2614842/93‑1 (MQ=255)

                                                                               cAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATta          <  1:5167495/93‑1 (MQ=255)

                                                                                 aCGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATtata        <  2:1960414/93‑1 (MQ=255)

                                                                                 aCGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATtata        >  2:6139496/1‑93 (MQ=255)

                                                                                  cGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAg       <  1:249269/93‑1 (MQ=255)

                                                                                   tttgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGa      <  2:462481/92‑1 (MQ=255)

                                                                                   gTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGa      >  1:341739/1‑93 (MQ=255)

                                                                                   gTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGa      <  1:1603075/93‑1 (MQ=255)

                                                                                      gcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATcc   <  1:1222084/93‑1 (MQ=255)

                                                                                      gcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATcc   <  1:3634637/93‑1 (MQ=255)

GCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCC  >  NC_000913/657724‑657903

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