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	97.2%	1218.9 / 13.2	398	I5M (ATA→ATG)	glnG	GlnG
Reads supporting (aligned to +/- strand): ref base T (1/10); new base C (203/184); total (204/194)
Fisher's exact test for biased strand distribution p-value = 4.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 408 total aligned reads displayed.

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTT  >  NC_000913/4052754‑4052934

cttCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCCCCCAACGGATGGAACTATCGTCATCGACTACCCAGACCAt                                                                                          <  2:8550008/92‑1 (MQ=255)

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

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

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

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

   ctcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCTACTATCCAGCCCATccc                                                                                       >  2:8702052/1‑93 (MQ=255)

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

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

       aaCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGt                                                                                   <  1:5323321/93‑1 (MQ=255)

       aaCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGt                                                                                   >  1:1695826/1‑93 (MQ=255)

          gtcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGc                                                                                <  1:7734487/93‑1 (MQ=255)

           tcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCa                                                                               >  2:9868614/1‑93 (MQ=255)

              tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaa                                                                            >  2:4726667/1‑93 (MQ=255)

              tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaa                                                                            >  1:2665777/1‑93 (MQ=255)

              tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaa                                                                            >  2:4909030/1‑93 (MQ=255)

               aCAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaga                                                                           >  1:1918978/1‑91 (MQ=255)

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

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

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

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

                 aGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACCGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACg                                                                         >  2:39862/1‑93 (MQ=255)

                    ttAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTcc                                                                      >  1:7490869/1‑92 (MQ=255)

                     tAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCAc                                                                     >  1:3412641/1‑93 (MQ=255)

                      aaaCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTCCCCAGACCATCCCTCGTTGCATAAACGTCAcc                                                                    <  2:2700912/93‑1 (MQ=255)

                       aaCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCt                                                                   <  1:9689355/93‑1 (MQ=255)

                        aCCTGCCCCAGCGAGTGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCtt                                                                  <  2:1996138/93‑1 (MQ=255)

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

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

                             ccccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTAttt                                                             >  2:6449915/1‑93 (MQ=255)

                              cccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACTTTTATTTc                                                            >  2:8129543/1‑93 (MQ=255)

                                cAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGCCCATCCCTCGTTGCATAAACGTCACCTTTATTTCCt                                                          <  1:1261213/93‑1 (MQ=255)

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

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

                                  gcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGa                                                        <  1:2216557/93‑1 (MQ=255)

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

                                  gcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGa                                                        <  1:1316541/93‑1 (MQ=255)

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

                                     agcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCCTCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAg                                                     >  2:10016796/1‑93 (MQ=255)

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

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

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

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

                                            cttCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTa                                              <  1:5458134/92‑1 (MQ=255)

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

                                             ttCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTaa                                             <  2:3247770/93‑1 (MQ=255)

                                               cAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAc                                           >  2:843763/1‑93 (MQ=255)

                                                   cACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgcg                                       >  1:1998872/1‑91 (MQ=255)

                                                   cACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgag                                       >  2:3490805/1‑93 (MQ=255)

                                                    ccccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgaga                                      <  1:8349707/92‑1 (MQ=255)

                                                     cccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAACGAGaa                                     >  2:2832578/1‑93 (MQ=255)

                                                       cAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGCACt                                   >  1:4935860/1‑93 (MQ=255)

                                                       cAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACt                                   >  2:4151657/1‑93 (MQ=255)

                                                       cAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTACTGATAGGCAGGTAAAATGAGAACt                                   >  1:6919721/1‑93 (MQ=255)

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

                                                        aaCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTc                                  >  1:7856538/1‑93 (MQ=255)

                                                          cGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCgg                                <  2:3392782/93‑1 (MQ=255)

                                                           ggATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGGACTCGGt                               >  1:7271325/1‑93 (MQ=255)

                                                           ggATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGt                               >  1:3684158/1‑93 (MQ=255)

                                                           ggATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGt                               >  1:6272355/1‑93 (MQ=255)

                                                            gATGTAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                              <  2:3751760/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                gAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGcc                          >  1:9769928/1‑93 (MQ=255)

                                                                 aaCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGccc                         >  2:2424209/1‑93 (MQ=255)

                                                                  aCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCt                        >  1:5469903/1‑93 (MQ=255)

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

                                                                   cTATCGTCATCGACTACCCAGCCCATTCCTCCTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCCAGAACGCGGTATGcccag                       >  2:8248528/1‑91 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTg                       <  1:8868890/93‑1 (MQ=255)

                                                                    tATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTgg                      >  1:8113927/1‑93 (MQ=255)

                                                                    tATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTgg                      >  2:7560739/1‑93 (MQ=255)

                                                                     aTCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGc                     >  1:9036215/1‑93 (MQ=255)

                                                                      tCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGcc                    >  1:4665408/1‑93 (MQ=255)

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

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

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

                                                                         tCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAAc                 >  2:1290453/1‑93 (MQ=255)

                                                                         tCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAAc                 <  2:2354771/93‑1 (MQ=255)

                                                                         tCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAAc                 >  2:295777/1‑93 (MQ=255)

                                                                          taTCCACTACCCAGACTATCCCTCGTTGCATAAACGTCCCCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACt                <  1:9200512/92‑1 (MQ=255)

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

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

                                                                            tCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTgg              >  1:990239/1‑93 (MQ=255)

                                                                             cGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGt             <  2:6769581/93‑1 (MQ=255)

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

                                                                               aCTACCCAGCCCATCCCTCGTTGAATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa           <  2:5604614/93‑1 (MQ=255)

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

                                                                               aCTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa           <  1:9364495/93‑1 (MQ=255)

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

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt         <  2:361428/93‑1 (MQ=255)

                                                                                   cccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTc       <  2:7821159/93‑1 (MQ=255)

                                                                                   cccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTc       >  1:5998544/1‑93 (MQ=255)

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

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

                                                                                      aGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGACATCGGTATGCCCTGGCCAACTGGTAAATTCAAt    >  2:5023283/1‑93 (MQ=255)

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

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

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTT  >  NC_000913/4052754‑4052934

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