BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,409,738	T→C	A130A (GCA→GCG)	yfbT ←	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,409,738	0	T	C	99.5%	652.1 / NA	209	A130A (GCA→GCG)	yfbT	putative phosphatase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (87/121); total (88/121)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

CGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTCAA  >  NC_000913/2409648‑2409827

cgcAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGgcgc                                                                                         <  1:1483248/93‑1 (MQ=255)

  cAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCgg                                                                                       >  2:2716580/1‑93 (MQ=255)

   aaGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCggg                                                                                      >  1:5645524/1‑93 (MQ=255)

   aaGCCCCAGCAGCTGCGAGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCAGTTACAAACACCTCTGGCGCggg                                                                                      <  1:6075683/93‑1 (MQ=255)

     gCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGaa                                                                                    <  2:527618/93‑1 (MQ=255)

     gCCCCAGCAGCCGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCAATCGCTCAGCGGTTACAAACACCTCTGGCGCGGGaa                                                                                    <  1:1383604/93‑1 (MQ=255)

      ccccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAg                                                                                   <  2:6016440/93‑1 (MQ=255)

       cccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGc                                                                                  <  1:3543103/93‑1 (MQ=255)

       cccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGc                                                                                  >  2:1908673/1‑93 (MQ=255)

       cccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGc                                                                                  <  2:5839403/93‑1 (MQ=255)

        ccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGcc                                                                                 >  2:3255651/1‑93 (MQ=255)

        ccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGcc                                                                                 >  1:1826740/1‑93 (MQ=255)

        ccAGCAGCTGCGCACCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGcc                                                                                 <  1:3028082/93‑1 (MQ=255)

         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGccc                                                                                >  1:5725764/1‑93 (MQ=255)

         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGccc                                                                                >  2:426677/1‑93 (MQ=255)

          agcagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCa                                                                               <  1:5517779/93‑1 (MQ=255)

          agcagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCa                                                                               <  1:4610764/93‑1 (MQ=255)

          agcagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCa                                                                               >  2:1178421/1‑93 (MQ=255)

          agcagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCa                                                                               <  2:1564162/93‑1 (MQ=255)

          agcagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCa                                                                               <  2:4698795/93‑1 (MQ=255)

            cagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCGCTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGc                                                                             <  2:5594682/93‑1 (MQ=255)

            cagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGc                                                                             <  2:7055433/93‑1 (MQ=255)

            cagcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGc                                                                             <  2:7522052/93‑1 (MQ=255)

             agcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCt                                                                            >  1:6183815/1‑93 (MQ=255)

              gcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTa                                                                           >  1:7089847/1‑93 (MQ=255)

              gcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTa                                                                           <  1:2716293/93‑1 (MQ=255)

                tGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTAtt                                                                         <  1:7577352/93‑1 (MQ=255)

                 gcgcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAATCCCAGCTAttt                                                                        <  1:5028666/93‑1 (MQ=255)

                 gcgcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTAttt                                                                        >  2:358481/1‑93 (MQ=255)

                 gcgcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCACTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTAttt                                                                        <  2:7007858/93‑1 (MQ=255)

                  cgcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTAtttt                                                                       <  1:1066740/93‑1 (MQ=255)

                   gcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACCCCTCTGGCGCGGGAAGCCCAGCTATTTTa                                                                      <  1:5497430/93‑1 (MQ=255)

                   gcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTa                                                                      <  1:1646865/93‑1 (MQ=255)

                   gcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTa                                                                      >  1:3169112/1‑93 (MQ=255)

                   gcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTa                                                                      <  2:6472077/93‑1 (MQ=255)

                    cgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTAt                                                                     <  2:3121854/93‑1 (MQ=255)

                    cgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTAt                                                                     <  2:1982390/93‑1 (MQ=255)

                       cTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgc                                                                  >  1:4704901/1‑93 (MQ=255)

                       cTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgc                                                                  >  1:3277969/1‑93 (MQ=255)

                       cTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgc                                                                  <  2:4306462/93‑1 (MQ=255)

                        tAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcg                                                                 >  2:2056843/1‑93 (MQ=255)

                           cAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcgc                                                              <  1:2942788/93‑1 (MQ=255)

                           cAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcgc                                                              >  1:6027158/1‑93 (MQ=255)

                            aGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgct                                                             >  2:2796095/1‑93 (MQ=255)

                            aGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgct                                                             <  1:1109575/93‑1 (MQ=255)

                             gATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTAGCGCGGGAAGCCCAGCTATTTTATGCtgcgctc                                                            >  1:2351080/1‑93 (MQ=255)

                                aCGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcg                                                         <  2:7781084/93‑1 (MQ=255)

                                aCGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcg                                                         >  2:3984719/1‑93 (MQ=255)

                                 cGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcgc                                                        <  1:4552357/93‑1 (MQ=255)

                                 cGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcgc                                                        <  1:3510735/93‑1 (MQ=255)

                                   cATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTa                                                      >  1:3768704/1‑93 (MQ=255)

                                   cATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTa                                                      >  1:4542778/1‑93 (MQ=255)

                                    aTCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAc                                                     >  2:6159729/1‑93 (MQ=255)

                                    aTCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAc                                                     <  1:7640002/93‑1 (MQ=255)

                                    aTCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAc                                                     <  2:4479381/93‑1 (MQ=255)

                                     tCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAcc                                                    >  1:7306675/1‑93 (MQ=255)

                                         gTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCa                                                <  1:175949/93‑1 (MQ=255)

                                           tCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATg                                              >  2:3909355/1‑93 (MQ=255)

                                           tCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATg                                              <  2:5079351/93‑1 (MQ=255)

                                           tCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATg                                              <  2:1432797/93‑1 (MQ=255)

                                             tGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAATCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGa                                            <  2:4035425/93‑1 (MQ=255)

                                               gTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGc                                          <  1:3632432/93‑1 (MQ=255)

                                               gTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGc                                          >  1:3836291/1‑93 (MQ=255)

                                                ttttCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGcc                                         <  2:1689077/93‑1 (MQ=255)

                                                  ttCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAg                                       >  1:2460335/1‑93 (MQ=255)

                                                     cgcgCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGTGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAg                                    <  1:4327926/93‑1 (MQ=255)

                                                     cgcgCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAg                                    >  2:7637094/1‑93 (MQ=255)

                                                        gcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTca                                 <  1:4076832/93‑1 (MQ=255)

                                                        gcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGAGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTca                                 >  2:3843886/1‑93 (MQ=255)

                                                         cTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTcac                                <  1:3331484/93‑1 (MQ=255)

                                                          ttCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTcaca                               >  2:36848/1‑93 (MQ=255)

                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                              <  2:6962798/93‑1 (MQ=255)

                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                              <  1:687030/93‑1 (MQ=255)

                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                              >  2:2242922/1‑93 (MQ=255)

                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                              <  1:1873147/93‑1 (MQ=255)

                                                             aCTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATg                            <  1:4600064/93‑1 (MQ=255)

                                                             aCTCGCTCAGCGGTTACAAAAACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATg                            <  2:4157985/93‑1 (MQ=255)

                                                               tCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGc                          >  2:6375213/1‑93 (MQ=255)

                                                               tCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGc                          >  2:5247129/1‑93 (MQ=255)

                                                                 gCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGccc                        <  2:3765061/93‑1 (MQ=255)

                                                                 gCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGccc                        <  1:2285209/93‑1 (MQ=255)

                                                                  cTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCa                       <  1:774862/93‑1 (MQ=255)

                                                                  cTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCa                       <  2:59128/93‑1 (MQ=255)

                                                                   tCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCAc                      >  2:3340759/1‑93 (MQ=255)

                                                                    cAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACg                     >  2:103767/1‑93 (MQ=255)

                                                                     aGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACgg                    >  1:1591238/1‑93 (MQ=255)

                                                                        ggTTACAACCACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAAt                 <  1:1446225/93‑1 (MQ=255)

                                                                         gTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATa                <  2:6657917/93‑1 (MQ=255)

                                                                         gTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATa                <  1:3906489/93‑1 (MQ=255)

                                                                          ttACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATAc               >  1:3691719/1‑93 (MQ=255)

                                                                          gtACAAACCCCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATAc               <  2:2182733/92‑1 (MQ=255)

                                                                           tACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATAcc              <  1:743629/93‑1 (MQ=255)

                                                                            aCAAACACCTCTGGCGCGGGAATCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCt             <  2:7803047/93‑1 (MQ=255)

                                                                              aagcacCTCTGGCGCGGGTGTCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGc           <  1:6080845/90‑1 (MQ=255)

                                                                              aaaCACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGc           >  1:6950003/1‑93 (MQ=255)

                                                                                   ccTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTAt      <  2:5074680/93‑1 (MQ=255)

                                                                                   ccTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTAt      >  1:4492165/1‑93 (MQ=255)

                                                                                   ccTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTAt      <  2:1953763/93‑1 (MQ=255)

                                                                                   ccTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTAt      <  1:5207943/93‑1 (MQ=255)

                                                                                       tGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTCaa  >  1:2566760/1‑93 (MQ=255)

CGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTCAA  >  NC_000913/2409648‑2409827

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