BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,052,844	T→C	I5M (ATA→ATG)	glnG ←	GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,052,844	0	T	C	99.6%	779.3 / ‑4.8	244	I5M (ATA→ATG)	glnG	GlnG
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (114/128); total (115/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 244 total aligned reads displayed.

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATT  >  NC_000913/4052754‑4052933

gctctcAAACGTCGTACAGGTTAAACCTCCCCTAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTTCCCAGACCAt                                                                                         <  2:1970909/91‑1 (MQ=255)

gTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCAt                                                                                         >  1:1667183/1‑93 (MQ=255)

 ttCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATc                                                                                        >  2:3814863/1‑93 (MQ=255)

 ttCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATc                                                                                        >  2:3199776/1‑93 (MQ=255)

 ttCTCAAACGTCGTACAGGATAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATc                                                                                        >  2:6128914/1‑93 (MQ=255)

   ctcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATccc                                                                                      >  1:3666940/1‑93 (MQ=255)

   ctcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATccc                                                                                      >  1:1919639/1‑93 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    >  2:3518322/1‑93 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    >  2:5531975/1‑93 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    <  2:2978228/93‑1 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    <  1:2035475/93‑1 (MQ=255)

      aaaCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCg                                                                                   <  2:2665809/93‑1 (MQ=255)

      aaaCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCg                                                                                   >  2:1056830/1‑93 (MQ=255)

        aCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGtt                                                                                 >  2:2759312/1‑93 (MQ=255)

         cgtcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTg                                                                                <  2:1653769/93‑1 (MQ=255)

             gtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATa                                                                            <  1:2235943/93‑1 (MQ=255)

              tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCAtga                                                                           >  2:2263319/1‑91 (MQ=255)

               aCAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaaa                                                                          <  2:639443/93‑1 (MQ=255)

                cAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAAc                                                                         <  2:24137/93‑1 (MQ=255)

                cAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAAc                                                                         >  2:4116780/1‑93 (MQ=255)

                 aGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACg                                                                        <  1:3258645/93‑1 (MQ=255)

                  ggTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACACAACGGATGGAACTAACGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGt                                                                       <  1:6298552/93‑1 (MQ=255)

                   gTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTc                                                                      <  2:2556467/93‑1 (MQ=255)

                    ttAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCa                                                                     <  2:3213573/93‑1 (MQ=255)

                      aaaCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCAcc                                                                   <  1:46267/93‑1 (MQ=255)

                       aaCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATACACGTCACCt                                                                  >  1:3282748/1‑93 (MQ=255)

                        aCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCtt                                                                 <  1:6143512/93‑1 (MQ=255)

                        aCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCtt                                                                 >  1:5302563/1‑93 (MQ=255)

                          cTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTa                                                               <  1:2657007/93‑1 (MQ=255)

                          cTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTa                                                               >  2:5132873/1‑93 (MQ=255)

                            gCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTAtt                                                             <  2:2205095/93‑1 (MQ=255)

                            gCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTAtt                                                             >  1:247753/1‑93 (MQ=255)

                              cccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTc                                                           >  1:3960216/1‑93 (MQ=255)

                              cccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTc                                                           >  2:695682/1‑93 (MQ=255)

                               ccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTct                                                          >  2:604796/1‑92 (MQ=255)

                                cAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCt                                                         >  1:2403874/1‑93 (MQ=255)

                                 agcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTg                                                        >  2:1853737/1‑93 (MQ=255)

                                 agcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTg                                                        <  1:3905616/93‑1 (MQ=255)

                                  gcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGa                                                       >  2:5628761/1‑93 (MQ=255)

                                   cgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGAt                                                      <  1:3083532/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAAGTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  2:1266406/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  2:2445168/93‑1 (MQ=255)

                                      gcgcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAgg                                                   <  2:1452260/93‑1 (MQ=255)

                                      gcgcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAgg                                                   >  1:4052938/1‑93 (MQ=255)

                                       cgcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGc                                                  <  1:3745776/93‑1 (MQ=255)

                                       cgcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGc                                                  >  2:4541159/1‑93 (MQ=255)

                                          aCGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAgg                                               >  2:3638707/1‑93 (MQ=255)

                                           cGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGt                                              >  1:148580/1‑93 (MQ=255)

                                            gTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTa                                             >  2:4822406/1‑93 (MQ=255)

                                             ttCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTaa                                            >  2:1006820/1‑93 (MQ=255)

                                               cAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAc                                          <  2:3894048/93‑1 (MQ=255)

                                               cAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAc                                          <  2:622066/93‑1 (MQ=255)

                                                aaGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAcc                                         <  1:4760622/93‑1 (MQ=255)

                                                   cACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgag                                      <  1:2608945/93‑1 (MQ=255)

                                                    aCCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCCGATAGGCAGGTAAACCgaga                                     >  2:4195051/1‑93 (MQ=255)

                                                     cccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGaa                                    <  1:717381/93‑1 (MQ=255)

                                                     cccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGaa                                    >  2:5040672/1‑93 (MQ=255)

                                                      ccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAc                                   >  1:1071633/1‑93 (MQ=255)

                                                      ccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAc                                   <  1:6074516/93‑1 (MQ=255)

                                                      ccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAc                                   >  1:5939571/1‑93 (MQ=255)

                                                       cAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACt                                  <  1:4588029/93‑1 (MQ=255)

                                                        aaCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTc                                 <  2:2925335/93‑1 (MQ=255)

                                                         aCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCg                                <  2:162866/93‑1 (MQ=255)

                                                         aCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCg                                <  2:4911528/93‑1 (MQ=255)

                                                         aCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCg                                >  1:5406545/1‑93 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             >  2:19111/1‑93 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             >  1:4277660/1‑93 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             <  1:421245/93‑1 (MQ=255)

                                                             aTGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTAt                            <  1:135036/93‑1 (MQ=255)

                                                               ggAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGc                          >  1:5187732/1‑93 (MQ=255)

                                                               ggAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGc                          >  1:2158871/1‑93 (MQ=255)

                                                                gAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGcc                         <  1:1225579/93‑1 (MQ=255)

                                                                gAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGcc                         <  2:1727277/93‑1 (MQ=255)

                                                                 aaCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGccc                        <  2:3302443/93‑1 (MQ=255)

                                                                  aCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCt                       <  1:2676710/93‑1 (MQ=255)

                                                                  aCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCt                       <  1:316463/93‑1 (MQ=255)

                                                                      tCGTCATCGACTACCCAGATCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGcc                   >  2:6107004/1‑93 (MQ=255)

                                                                        gTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCaa                 >  1:6019603/1‑93 (MQ=255)

                                                                        gTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCaa                 <  2:5206783/93‑1 (MQ=255)

                                                                          cATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACt               <  2:1338322/93‑1 (MQ=255)

                                                                           aTCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTg              >  2:172254/1‑93 (MQ=255)

                                                                           aTCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTg              >  1:4898363/1‑93 (MQ=255)

                                                                            tCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTgg             <  2:761959/93‑1 (MQ=255)

                                                                            tCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTgg             <  1:2869669/93‑1 (MQ=255)

                                                                              gACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATTCCCTGGCCAACTggca           >  2:5283022/1‑91 (MQ=255)

                                                                              gACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTa           <  1:766337/93‑1 (MQ=255)

                                                                               aCTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa          >  2:6230149/1‑93 (MQ=255)

                                                                               aCTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa          >  1:1446704/1‑93 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        >  2:5014433/1‑93 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        >  1:6128303/1‑93 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        >  2:407222/1‑93 (MQ=255)

                                                                                  aCCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAtt       >  1:5618815/1‑93 (MQ=255)

                                                                                    ccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCa     >  1:585462/1‑93 (MQ=255)

                                                                                    ccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCa     >  2:879267/1‑93 (MQ=255)

                                                                                     cAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCaa    <  1:4116780/93‑1 (MQ=255)

                                                                                     cAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCaa    >  2:5867157/1‑93 (MQ=255)

                                                                                      aGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAt   >  1:4452941/1‑93 (MQ=255)

                                                                                      aGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAt   >  1:6201103/1‑93 (MQ=255)

                                                                                       gACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAtt  <  2:5801740/93‑1 (MQ=255)

                                                                                       gACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAtt  >  1:1600282/1‑93 (MQ=255)

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATT  >  NC_000913/4052754‑4052933

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