BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	866,644	A→G	N285S (AAT→AGT)	ybiK →	putative asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	866,644	0	A	G	99.2%	742.5 / ‑3.1	241	N285S (AAT→AGT)	ybiK	putative asparaginase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (133/106); total (134/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 253 total aligned reads displayed.

GTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866554‑866733

gTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTg                                                                                         >  2:5713233/1‑93 (MQ=255)

  ctcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTc                                                                                       <  1:701474/93‑1 (MQ=255)

   tcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcg                                                                                      >  2:3602681/1‑93 (MQ=255)

   tcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcg                                                                                      <  1:5438235/93‑1 (MQ=255)

    cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     >  2:1189370/1‑93 (MQ=255)

    cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     <  2:6660099/93‑1 (MQ=255)

    cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTcgc                                                                                     >  2:6335173/1‑93 (MQ=255)

     gcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgcg                                                                                    >  1:3691872/1‑93 (MQ=255)

       ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCt                                                                                  >  2:2253363/1‑93 (MQ=255)

           gCCTGCGAGCGGGTAGTAATGGAAAACCTCCCTGCGCTTGGCGGTAGCGGTGTCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCg                                                                              <  1:2854304/93‑1 (MQ=255)

           gCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCg                                                                              <  2:7721843/93‑1 (MQ=255)

           gCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCg                                                                              >  2:1633246/1‑93 (MQ=255)

             cTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGtt                                                                            <  2:6886306/93‑1 (MQ=255)

              tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           <  1:956125/93‑1 (MQ=255)

              tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           >  1:5254332/1‑93 (MQ=255)

              tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           >  2:4864609/1‑93 (MQ=255)

              tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           >  1:6400112/1‑93 (MQ=255)

               gCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTa                                                                          <  2:3382543/93‑1 (MQ=255)

                  aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       <  1:1047374/93‑1 (MQ=255)

                  aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       >  2:1340502/1‑93 (MQ=255)

                    cGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACAcc                                                                     <  1:3080753/93‑1 (MQ=255)

                    cGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACAcc                                                                     >  1:9139/1‑93 (MQ=255)

                      ggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGa                                                                   <  2:1427875/93‑1 (MQ=255)

                       gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGATGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  <  1:6622751/93‑1 (MQ=255)

                        tagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAg                                                                 >  1:6204663/1‑93 (MQ=255)

                          gtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGa                                                               >  2:5086198/1‑93 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCTCGCTACCGTTTAACACCGAAGGaa                                                              >  2:7204696/1‑93 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              >  1:4030772/1‑93 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              >  1:547443/1‑93 (MQ=255)

                            aaTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAAt                                                             >  2:3341131/1‑93 (MQ=255)

                            aaTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAAt                                                             <  1:2005113/93‑1 (MQ=255)

                             aTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATg                                                            <  2:7580712/93‑1 (MQ=255)

                              tGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGt                                                           <  1:7680343/93‑1 (MQ=255)

                              tGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGGGTCGCGCTACCGTTTAACACCGAAGGAATGt                                                           >  2:4111372/1‑93 (MQ=255)

                               ggAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTa                                                          <  1:6569103/93‑1 (MQ=255)

                                gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         >  1:5530882/1‑93 (MQ=255)

                                gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         <  1:4656432/93‑1 (MQ=255)

                                gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         <  2:5187901/93‑1 (MQ=255)

                                 aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        <  1:1991532/93‑1 (MQ=255)

                                  aaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcg                                                       >  1:2687796/1‑93 (MQ=255)

                                   aaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcgc                                                      <  2:7115428/93‑1 (MQ=255)

                                   aaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcgc                                                      >  2:7123439/1‑93 (MQ=255)

                                     aCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCgcgc                                                    >  1:3124296/1‑93 (MQ=255)

                                     aCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCgcgc                                                    <  2:7472731/93‑1 (MQ=255)

                                      cTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGcc                                                   <  1:5931870/93‑1 (MQ=255)

                                        cccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTg                                                 <  1:2431810/93‑1 (MQ=255)

                                        cccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTg                                                 <  1:1388026/93‑1 (MQ=255)

                                         ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                >  2:4209949/1‑93 (MQ=255)

                                         ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                >  2:3562294/1‑93 (MQ=255)

                                          cTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTggg                                               >  1:4447426/1‑93 (MQ=255)

                                           tGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgggg                                              >  1:1917874/1‑93 (MQ=255)

                                            gcgcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGc                                             >  2:2964958/1‑93 (MQ=255)

                                              gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTa                                           >  1:1483728/1‑93 (MQ=255)

                                              gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTa                                           >  1:4064021/1‑93 (MQ=255)

                                               cTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTAc                                          <  2:3938262/93‑1 (MQ=255)

                                                ttGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACg                                         >  2:7366783/1‑93 (MQ=255)

                                                      gTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCg                                   >  1:2030213/1‑93 (MQ=255)

                                                       tAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGa                                  <  1:3197213/93‑1 (MQ=255)

                                                       tAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGCATCGCGACTGGGGCTCCGCCGGGGa                                  >  1:2223447/1‑93 (MQ=255)

                                                         gCGGTGTCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  1:2204283/93‑1 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                >  2:4715405/1‑93 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  2:522393/93‑1 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  1:7845324/93‑1 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                >  1:7248332/1‑93 (MQ=255)

                                                          cGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATAc                               >  2:3073147/1‑93 (MQ=255)

                                                           ggTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACg                              >  2:2597712/1‑93 (MQ=255)

                                                           ggTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACg                              >  1:909479/1‑93 (MQ=255)

                                                             tGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGcc                            >  2:305482/1‑93 (MQ=255)

                                                              ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           >  2:2367108/1‑93 (MQ=255)

                                                              ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           >  1:2479487/1‑93 (MQ=255)

                                                              ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           >  1:6163677/1‑93 (MQ=255)

                                                               gCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCaa                          <  2:2715506/93‑1 (MQ=255)

                                                                cTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAac                         <  1:534009/93‑1 (MQ=255)

                                                                 ttAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacc                        <  1:4433026/93‑1 (MQ=255)

                                                                 ttAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacc                        >  1:7412121/1‑93 (MQ=255)

                                                                  tAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacca                       >  1:4822290/1‑93 (MQ=255)

                                                                  tAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacca                       >  2:245737/1‑93 (MQ=255)

                                                                   aaTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccac                      <  2:1914734/93‑1 (MQ=255)

                                                                    aTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccacc                     >  1:6916248/1‑93 (MQ=255)

                                                                     tCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCg                    >  1:3870024/1‑93 (MQ=255)

                                                                       gCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGt                  <  2:4026275/93‑1 (MQ=255)

                                                                       gCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGt                  >  2:5869479/1‑93 (MQ=255)

                                                                         tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                >  1:3474516/1‑93 (MQ=255)

                                                                         tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                >  2:62248/1‑93 (MQ=255)

                                                                         tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                >  1:1646538/1‑93 (MQ=255)

                                                                         tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACCCCAACCACCGGTAt                >  2:6722015/1‑93 (MQ=255)

                                                                          aTCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATc               >  2:6113183/1‑93 (MQ=255)

                                                                          aTCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATc               <  2:6139949/93‑1 (MQ=255)

                                                                          aTCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATc               >  1:5157941/1‑93 (MQ=255)

                                                                           tCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCt              >  1:1597481/1‑93 (MQ=255)

                                                                             gACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTAc            <  2:1668370/93‑1 (MQ=255)

                                                                             gACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTAc            >  1:6503284/1‑93 (MQ=255)

                                                                              aCCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTAcc           <  1:5633692/93‑1 (MQ=255)

                                                                                cATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGt         <  2:4906746/93‑1 (MQ=255)

                                                                                cATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGt         <  1:1171341/93‑1 (MQ=255)

                                                                                  tGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGa       >  1:57985/1‑93 (MQ=255)

                                                                                     aGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaa    <  2:551841/93‑1 (MQ=255)

                                                                                      gggAGTGTCGCGCTACCGTTTAACAGCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCAATACACTAACCACCGGTATCTACCGTGaaaaa   <  1:2633875/93‑1 (MQ=255)

                                                                                      gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   <  2:6879363/93‑1 (MQ=255)

                                                                                      gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   <  1:7806925/93‑1 (MQ=255)

GTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866554‑866733

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