BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,647,829	T→C	T827A (ACG→GCG)	b2520 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,647,829	0	T	C	84.5%	568.7 / 91.0	226	T827A (ACG→GCG)	b2520	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (15/20); new base C (91/100); total (106/120)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 229 total aligned reads displayed.

GGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACC  >  NC_000913/2647739‑2647918

ggTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCAc                                                                                         >  2:3526946/1‑93 (MQ=255)

 gTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACt                                                                                        >  2:4089672/1‑93 (MQ=255)

  ttttCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACtt                                                                                       >  2:89203/1‑93 (MQ=255)

   tttCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTa                                                                                      <  1:2893378/93‑1 (MQ=255)

   tttCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTa                                                                                      >  1:4499674/1‑93 (MQ=255)

   tttCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTa                                                                                      >  1:1610226/1‑93 (MQ=255)

      cATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAAt                                                                                   <  1:4865958/93‑1 (MQ=255)

       aTCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATa                                                                                  >  2:2383636/1‑93 (MQ=255)

        tCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATAc                                                                                 >  1:443455/1‑93 (MQ=255)

        ccgcgAAGCGGCAGATGCAACTCACCACCCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATAc                                                                                 <  1:4098537/92‑1 (MQ=255)

         cGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACa                                                                                <  2:6624937/93‑1 (MQ=255)

         cGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACa                                                                                <  2:3522614/93‑1 (MQ=255)

          gCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAg                                                                               <  1:5169310/93‑1 (MQ=255)

          gCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAg                                                                               >  1:744285/1‑93 (MQ=255)

            cAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGAt                                                                             <  2:5049884/93‑1 (MQ=255)

            cAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGAt                                                                             >  2:5963845/1‑93 (MQ=255)

             aaGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATc                                                                            <  2:3592995/93‑1 (MQ=255)

             aaGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATc                                                                            <  2:3241134/93‑1 (MQ=255)

             aaGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATc                                                                            <  1:5043593/93‑1 (MQ=255)

              aGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCa                                                                           <  2:5668309/93‑1 (MQ=255)

               gCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGGGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATTCAGATCAt                                                                          <  1:6503799/93‑1 (MQ=255)

                 ggCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGa                                                                        >  2:3961917/1‑93 (MQ=255)

                 ggCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGa                                                                        <  1:456841/93‑1 (MQ=255)

                    aGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGa                                                                     >  2:6464125/1‑93 (MQ=255)

                     gATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCTTGACGAt                                                                    <  2:2908977/93‑1 (MQ=255)

                     gATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGAt                                                                    >  2:5788486/1‑93 (MQ=255)

                      aTGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGAtt                                                                   <  1:6234995/93‑1 (MQ=255)

                       tGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTc                                                                  >  1:579619/1‑93 (MQ=255)

                       tGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTc                                                                  <  1:589454/93‑1 (MQ=255)

                        gCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTcc                                                                 <  2:6360149/93‑1 (MQ=255)

                        gCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTcc                                                                 >  1:2562902/1‑93 (MQ=255)

                         cAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCa                                                                >  2:2446552/1‑93 (MQ=255)

                           acacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAgg                                                              <  2:6414688/93‑1 (MQ=255)

                            cacacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGt                                                             <  1:1427029/93‑1 (MQ=255)

                              cacCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGttt                                                           <  2:713033/93‑1 (MQ=255)

                               cccAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATCCCAGGtttt                                                          <  2:1156387/92‑1 (MQ=255)

                               acCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGtttt                                                          >  1:4985546/1‑93 (MQ=255)

                                  aaCCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATc                                                       >  1:938509/1‑93 (MQ=255)

                                   aCCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCg                                                      >  1:2678393/1‑93 (MQ=255)

                                   aCCGCGCGTTTTGGCGGCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCg                                                      <  1:3911435/93‑1 (MQ=255)

                                    ccGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGa                                                     >  1:1992280/1‑93 (MQ=255)

                                     cgcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATNATGACGATTCCAGGTTTTATCGAc                                                    >  1:1281413/1‑93 (MQ=255)

                                     cgcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGAc                                                    >  1:1087210/1‑93 (MQ=255)

                                      gcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGAcc                                                   >  2:4634842/1‑93 (MQ=255)

                                      gcgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGAcc                                                   <  2:1666628/93‑1 (MQ=255)

                                       cgcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCg                                                  <  1:222015/93‑1 (MQ=255)

                                        gcgTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCgg                                                 >  2:2069644/1‑93 (MQ=255)

                                         cgTTTTGGCGGCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGa                                                <  1:2046187/93‑1 (MQ=255)

                                          gTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGaa                                               <  2:3818675/93‑1 (MQ=255)

                                           ttttGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAAt                                              <  1:1562594/93‑1 (MQ=255)

                                           ttttGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAAt                                              <  1:2784951/93‑1 (MQ=255)

                                           ttttGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAAt                                              <  1:530788/93‑1 (MQ=255)

                                            tttGGCGGCGCGGAACGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATc                                             <  2:1272602/93‑1 (MQ=255)

                                              tGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGt                                           <  1:3648964/93‑1 (MQ=255)

                                              tGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGt                                           <  1:5257802/93‑1 (MQ=255)

                                               ggCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTc                                          <  2:5446293/93‑1 (MQ=255)

                                               ggCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTc                                          <  1:2402316/93‑1 (MQ=255)

                                               ggCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTc                                          <  2:3922580/93‑1 (MQ=255)

                                               ggCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTc                                          <  2:1728907/93‑1 (MQ=255)

                                                gCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCa                                         <  1:794882/93‑1 (MQ=255)

                                                gCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCa                                         >  2:1130617/1‑93 (MQ=255)

                                                     gcgGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATcc                                    <  1:5642405/93‑1 (MQ=255)

                                                     gcgGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATcc                                    >  1:3100316/1‑93 (MQ=255)

                                                     gcgGAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATcc                                    <  2:4482265/93‑1 (MQ=255)

                                                       ggAGCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAg                                  <  1:4618524/93‑1 (MQ=255)

                                                          gCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGccc                               <  2:1847944/93‑1 (MQ=255)

                                                          gCGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGccc                               <  2:2160712/93‑1 (MQ=255)

                                                           cGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGccat                              <  1:6549908/93‑3 (MQ=255)

                                                           cGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCt                              >  2:4255337/1‑93 (MQ=255)

                                                           cGAGATTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCt                              <  1:1692691/93‑1 (MQ=255)

                                                            gagaTTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCtt                             <  1:1021469/93‑1 (MQ=255)

                                                             agaTTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTg                            >  1:2713014/1‑93 (MQ=255)

                                                              gaTTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGa                           >  2:4506056/1‑93 (MQ=255)

                                                              gaTTTATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGa                           >  1:3237818/1‑93 (MQ=255)

                                                                 ttATCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGcc                        >  1:967327/1‑93 (MQ=255)

                                                                  tATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGccc                       <  2:3036074/93‑1 (MQ=255)

                                                                   aTCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCg                      >  2:4919907/1‑93 (MQ=255)

                                                                   aTCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCg                      >  1:3431323/1‑93 (MQ=255)

                                                                   aTCGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCg                      >  2:1421256/1‑93 (MQ=255)

                                                                     cGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGaa                    <  1:1893257/93‑1 (MQ=255)

                                                                     cGCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGaa                    <  1:4029863/93‑1 (MQ=255)

                                                                      gCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAAc                   >  2:2849939/1‑93 (MQ=255)

                                                                      gCCAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAAc                   <  2:6085695/93‑1 (MQ=255)

                                                                       ccAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACa                  <  1:5466279/93‑1 (MQ=255)

                                                                       ccAGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACa                  >  2:4745219/1‑93 (MQ=255)

                                                                         aGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATc                >  1:3302817/1‑93 (MQ=255)

                                                                         aGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATc                >  2:4818619/1‑93 (MQ=255)

                                                                         aGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATc                <  2:4124398/93‑1 (MQ=255)

                                                                         aGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATc                <  1:3535473/93‑1 (MQ=255)

                                                                         aGGACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATc                >  2:2248795/1‑93 (MQ=255)

                                                                           gACGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCaa              >  2:2089438/1‑93 (MQ=255)

                                                                            aCGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAt             <  2:5263597/93‑1 (MQ=255)

                                                                             cGTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAtc            >  1:6009787/1‑93 (MQ=255)

                                                                              gTACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAtct           <  2:3743386/93‑1 (MQ=255)

                                                                               tACCACCAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAAtctc          <  1:5596137/93‑1 (MQ=255)

                                                                                accaccAGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATctct         <  2:3620458/93‑1 (MQ=255)

                                                                                    ccaGCGCACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGcc     >  2:3076025/1‑93 (MQ=255)

                                                                                      aGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGTATCGTCAGATTCAGCCCTTGAGACCGAACATCAATCTCTTGccac   >  1:6643873/1‑93 (MQ=255)

                                                                                       gcgcACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGccacc  <  1:2303257/93‑1 (MQ=255)

                                                                                       gcgcACTCAAATACAGACCATGACGATTCCAGGTTCTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGccacc  <  1:3619422/93‑1 (MQ=255)

GGTTTTCATCGCCAAGCGGCAGATGCAACACACCAACCGCGCGTTTTGGCGTCGCGGAGCGAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACC  >  NC_000913/2647739‑2647918

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