BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,516,689	T→C	100%	T382A (ACG→GCG)	yfeA ←	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,516,689	0	T	C	100.0%	515.3 / ‑6.7	147	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (82/64); minor base G (1/0); total (83/64)

Only 100 of 147 total aligned reads displayed.

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784

aGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACCGCAGTTCACTACCCGGAAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAc                                                                                               >  2:470589/1‑100 (MQ=255)

aGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAc                                                                                               <  2:1684771/100‑1 (MQ=255)

aGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAc                                                                                               >  1:197929/1‑100 (MQ=255)

 gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              <  1:1796392/100‑1 (MQ=255)

  tCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACgg                                                                                             <  2:247150/100‑1 (MQ=255)

  tCGCCCTTCCGGTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACGTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACgg                                                                                             <  2:699913/100‑1 (MQ=255)

   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            <  2:2440951/100‑1 (MQ=255)

      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         >  2:3242495/1‑100 (MQ=255)

      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         <  2:3524767/100‑1 (MQ=255)

        ttCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAtt                                                                                       >  2:2416757/1‑100 (MQ=255)

        ttCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAtt                                                                                       >  2:3649635/1‑100 (MQ=255)

          ccGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGa                                                                                     >  1:919886/1‑100 (MQ=255)

           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  2:1811442/1‑100 (MQ=255)

            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   >  2:1269936/1‑100 (MQ=255)

             tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCTGCAGCGCACGGCAAATTGAGCg                                                                                  <  2:68583/100‑1 (MQ=255)

             tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  1:787581/100‑1 (MQ=255)

                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               >  2:624487/1‑100 (MQ=255)

                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               >  1:393986/1‑100 (MQ=255)

                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               <  1:2549549/100‑1 (MQ=255)

                ccGGCCCGCGCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTCTCGTTTTCCTGCATCAACCGCAGCAGCGCAATGCGAATTGAGCGGAt                                                                               >  1:1661019/1‑100 (MQ=255)

                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  2:1039635/1‑100 (MQ=255)

                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:2708260/1‑100 (MQ=255)

                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             >  2:3678678/1‑100 (MQ=255)

                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             <  2:1490910/100‑1 (MQ=255)

                   gCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                            >  1:3592852/1‑100 (MQ=255)

                    cccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAg                                                                           <  1:3210175/100‑1 (MQ=255)

                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                         <  1:308647/100‑1 (MQ=255)

                       gCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGa                                                                        >  1:1546709/1‑100 (MQ=255)

                          cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                     >  1:216175/1‑100 (MQ=255)

                          cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                     >  1:215147/1‑100 (MQ=255)

                             cACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACgcgc                                                                  >  1:2935381/1‑100 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  2:1536374/100‑1 (MQ=255)

                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcat                                                                >  1:856710/1‑100 (MQ=255)

                                tAACAGCATTGCACTACCCGGCAATGGATACTACTTTTCGTTTTCCTGCATCAACGGCAGCAGCTCACGGCAAATTGAGCTGATACAGTGAGCGCGcatc                                                               <  1:1306642/100‑1 (MQ=255)

                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  1:1044052/1‑100 (MQ=255)

                                  aCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatcat                                                             >  2:1107047/1‑100 (MQ=255)

                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            >  1:2955594/1‑100 (MQ=255)

                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            >  1:1292317/1‑100 (MQ=255)

                                    agcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                           <  1:150307/100‑1 (MQ=255)

                                     gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          >  2:3734388/1‑100 (MQ=255)

                                     gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          <  1:1152184/100‑1 (MQ=255)

                                     gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAAAGGCAGCAGCGCACAGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          >  1:639658/1‑100 (MQ=255)

                                      cagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAg                                                         <  1:2610234/100‑1 (MQ=255)

                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  1:1066626/100‑1 (MQ=255)

                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  1:304459/100‑1 (MQ=255)

                                         ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:1921589/100‑1 (MQ=255)

                                          tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                     <  2:1879423/100‑1 (MQ=255)

                                          tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                     >  1:1204489/1‑100 (MQ=255)

                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    <  2:2343517/100‑1 (MQ=255)

                                              tACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATg                                                 <  2:1431610/100‑1 (MQ=255)

                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  2:1262687/1‑100 (MQ=255)

                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               <  1:470587/100‑1 (MQ=255)

                                                 ccGGCAATTTATACAAATTTTCGTTTTCCTGCATCAACGGCAGCAGCAGACGACAACTTGAGCGGATACAGTGTACGCGCATCAATAAGCCGTAATgacg                                              >  1:1489320/1‑100 (MQ=255)

                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              <  1:2356056/100‑1 (MQ=255)

                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              >  2:1946097/1‑100 (MQ=255)

                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                             <  1:3716679/100‑1 (MQ=255)

                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                             >  2:3034101/1‑100 (MQ=255)

                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                             >  1:2794393/1‑100 (MQ=255)

                                                    gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  2:2866620/100‑1 (MQ=255)

                                                      agttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTTAGCGGATACAGTGAACGCGCATCATTAATGCGTAATGACGACtga                                         >  1:23714/3‑98 (MQ=255)

                                                      aaTTGATACAACTTTTCGTTTTCCTGGATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  1:1770620/98‑1 (MQ=255)

                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                       <  2:3107104/100‑1 (MQ=255)

                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                     >  1:2776249/1‑100 (MQ=255)

                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  2:3504657/100‑1 (MQ=255)

                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  2:377823/100‑1 (MQ=255)

                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    >  1:1170701/1‑100 (MQ=255)

                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:1995553/1‑100 (MQ=255)

                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  2:1724274/1‑100 (MQ=255)

                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 <  2:2108421/100‑1 (MQ=255)

                                                               aaCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCaa                                >  1:586765/1‑100 (MQ=255)

                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               >  2:3753117/1‑100 (MQ=255)

                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  2:3300065/100‑1 (MQ=255)

                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              >  2:1327061/1‑100 (MQ=255)

                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                             <  2:478203/100‑1 (MQ=255)

                                                                    ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTa                           <  2:569149/100‑1 (MQ=255)

                                                                      cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  2:2131933/100‑1 (MQ=255)

                                                                       gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        >  2:2453221/1‑100 (MQ=255)

                                                                        ttttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCaa                       >  2:2678918/1‑100 (MQ=255)

                                                                          ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  2:776709/100‑1 (MQ=255)

                                                                          ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  1:1658611/100‑1 (MQ=255)

                                                                           tCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgc                    <  2:820244/100‑1 (MQ=255)

                                                                            ccTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcg                   >  2:2868450/1‑100 (MQ=255)

                                                                            ccTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcg                   >  2:292221/1‑100 (MQ=255)

                                                                              ttcATCAACGGCAGCAGCGCACGGCAAATGGAGTGGATACAGTGATCTCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 <  2:363285/98‑1 (MQ=255)

                                                                               gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                >  1:3672530/1‑100 (MQ=255)

                                                                                cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg               >  1:2870963/1‑100 (MQ=255)

                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  1:3338434/1‑100 (MQ=255)

                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  1:810458/1‑100 (MQ=255)

                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             >  1:989903/1‑100 (MQ=255)

                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             >  2:3514844/1‑100 (MQ=255)

                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             <  2:2633392/100‑1 (MQ=255)

                                                                                   cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  1:3542176/100‑1 (MQ=255)

                                                                                   cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            >  1:408961/1‑100 (MQ=255)

                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           <  2:2928903/100‑1 (MQ=255)

                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  1:1462322/1‑100 (MQ=255)

                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  1:3429599/1‑100 (MQ=255)

                                                                                           gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    <  1:1933506/100‑1 (MQ=255)

                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  >  2:2366989/1‑100 (MQ=255)

                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  >  2:2367117/1‑100 (MQ=255)

                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    <  1:60698/98‑1 (MQ=255)

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784

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