BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	199,879	A→G	N651S (AAC→AGC)	yaeT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	199,879	0	A	G	99.1%	1041.4 / ‑3.2	326	N651S (AAC→AGC)	yaeT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (146/177); total (147/179)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 346 total aligned reads displayed.

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCt                                                                                         <  1:4914606/93‑1 (MQ=255)

tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCt                                                                                         <  2:9168633/93‑1 (MQ=255)

 cgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCtt                                                                                        >  2:5307053/1‑93 (MQ=255)

  gatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTc                                                                                       <  1:6619915/93‑1 (MQ=255)

  gatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTc                                                                                       >  1:4447162/1‑93 (MQ=255)

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      >  1:5627709/1‑93 (MQ=255)

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      >  1:5196096/1‑93 (MQ=255)

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      <  1:8707662/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  2:6806723/93‑1 (MQ=255)

      aCGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATg                                                                                   >  1:1641578/1‑93 (MQ=255)

        gATCACAAATGTGTTGTTCTGTGGCGTACCCGCTGGGGGTATGGTGGTGGTGTATGCGGCAAAGAGATGACGATCTACTAGAGCTTCTATGcc                                                                                 >  2:5566122/1‑93 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               >  2:3682876/1‑93 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               <  1:3608541/93‑1 (MQ=255)

            acaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtg                                                                             <  1:2038987/93‑1 (MQ=255)

              aaaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtggt                                                                           <  1:7065195/93‑1 (MQ=255)

               aaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGtt                                                                          >  2:5214285/1‑93 (MQ=255)

                  gggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCa                                                                       <  1:9248250/93‑1 (MQ=255)

                  gggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCa                                                                       <  1:2397189/93‑1 (MQ=255)

                   ggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCGAg                                                                      >  2:6207138/1‑93 (MQ=255)

                   ggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAg                                                                      <  2:8662598/93‑1 (MQ=255)

                    gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGc                                                                     <  2:7013338/93‑1 (MQ=255)

                    gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGc                                                                     >  1:14210/1‑93 (MQ=255)

                      tgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCAc                                                                   <  1:9964485/93‑1 (MQ=255)

                      tgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCAc                                                                   <  1:10015294/93‑1 (MQ=255)

                       gttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATACCGGTGGTTCCAGCAcc                                                                  >  2:1803386/1‑93 (MQ=255)

                        ttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcg                                                                 <  1:5971886/93‑1 (MQ=255)

                        ttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcg                                                                 <  2:9286245/93‑1 (MQ=255)

                         tCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgt                                                                <  2:7619657/93‑1 (MQ=255)

                            ggggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcg                                                             <  1:3168886/93‑1 (MQ=255)

                            ggggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcg                                                             <  2:8080564/93‑1 (MQ=255)

                             gggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcgt                                                            >  2:9891396/1‑93 (MQ=255)

                             gggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcgt                                                            >  1:3386717/1‑93 (MQ=255)

                               gCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTgg                                                          <  1:6286981/93‑1 (MQ=255)

                                  tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCGAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCtt                                                       >  1:2123852/1‑93 (MQ=255)

                                  tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCtt                                                       >  2:5832358/1‑93 (MQ=255)

                                  tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCtt                                                       <  2:2117070/93‑1 (MQ=255)

                                    cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTcc                                                     <  1:1254077/93‑1 (MQ=255)

                                    cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTcc                                                     <  1:8864634/93‑1 (MQ=255)

                                    cccGCTGGGGTTATGCTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTcc                                                     <  2:5502238/93‑1 (MQ=255)

                                     ccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCa                                                    <  1:7792098/93‑1 (MQ=255)

                                      cGCTGGGGTTATGGTGATGGTTTAGGCGGCACAGAGATGCCGTTCTACGCGAGCTTCTATGCCGGTGGTTCCAGCACCGCGCGTGGCTACCAg                                                   >  2:201959/1‑93 (MQ=255)

                                       gCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGt                                                  >  1:3720200/1‑93 (MQ=255)

                                       gCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGt                                                  >  2:7320891/1‑93 (MQ=255)

                                        cTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTc                                                 <  2:6501322/93‑1 (MQ=255)

                                        cTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTc                                                 >  1:6081108/1‑93 (MQ=255)

                                          ggggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCa                                               <  2:5357065/93‑1 (MQ=255)

                                           gggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCaa                                              <  1:6412551/93‑1 (MQ=255)

                                            ggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAAt                                             <  2:7929668/93‑1 (MQ=255)

                                              ttATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATAc                                           >  1:5433198/1‑93 (MQ=255)

                                              ttATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCCTCCAGTCCCATAc                                           >  1:7586127/1‑93 (MQ=255)

                                               tagggTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATAcc                                          <  2:4888930/90‑1 (MQ=255)

                                               tATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATAcc                                          <  1:7887813/93‑1 (MQ=255)

                                                   gTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTg                                      <  2:2939765/93‑1 (MQ=255)

                                                    tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTgg                                     >  2:1309768/1‑93 (MQ=255)

                                                     gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGt                                    >  2:882408/1‑93 (MQ=255)

                                                     gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGt                                    <  2:3783388/93‑1 (MQ=255)

                                                      aTGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTc                                   <  2:5421931/93‑1 (MQ=255)

                                                       tGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTcc                                  >  2:4186346/1‑93 (MQ=255)

                                                          tttAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGaa                               >  1:9172639/1‑93 (MQ=255)

                                                             aGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGc                            >  1:2350032/1‑93 (MQ=255)

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCa                           <  1:6783654/93‑1 (MQ=255)

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCa                           <  1:7151317/93‑1 (MQ=255)

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCa                           <  1:6941967/93‑1 (MQ=255)

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCa                           <  2:189894/93‑1 (MQ=255)

                                                               gcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAg                          >  2:2408368/1‑93 (MQ=255)

                                                                 ggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                        <  2:4353129/93‑1 (MQ=255)

                                                                 ggcAAAGACATGCCGTTCTTCGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                        <  2:9487032/93‑1 (MQ=255)

                                                                  gcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGttt                       >  2:4720021/1‑93 (MQ=255)

                                                                  gcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGttt                       >  2:1214710/1‑93 (MQ=255)

                                                                   cAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTa                      <  2:8422177/93‑1 (MQ=255)

                                                                    caaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     <  1:4503993/92‑1 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTat                     >  1:2728980/1‑92 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     <  1:4646462/93‑1 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     >  1:1890155/1‑93 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     <  2:9474513/93‑1 (MQ=255)

                                                                     aaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACt                    <  1:5077156/93‑1 (MQ=255)

                                                                     aaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACt                    >  2:7647597/1‑93 (MQ=255)

                                                                      agagATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACtt                   <  2:3236538/93‑1 (MQ=255)

                                                                      agagATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACtt                   <  1:8121361/93‑1 (MQ=255)

                                                                      agagATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACtt                   >  2:9693157/1‑93 (MQ=255)

                                                                       gagaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTc                  <  2:5931247/93‑1 (MQ=255)

                                                                        agaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTGcc                 <  2:6722999/93‑1 (MQ=255)

                                                                         gaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTccc                >  1:9371979/1‑93 (MQ=255)

                                                                         gaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTccc                >  1:3894136/1‑93 (MQ=255)

                                                                         gaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTccc                >  1:262416/1‑93 (MQ=255)

                                                                         gaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTccc                <  2:2532696/93‑1 (MQ=255)

                                                                          aTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTcccc               >  1:4085817/1‑92 (MQ=255)

                                                                          aTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCg               >  1:8369754/1‑93 (MQ=255)

                                                                            gCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCa             <  2:8235299/93‑1 (MQ=255)

                                                                            gCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCa             <  1:5528541/93‑1 (MQ=255)

                                                                              cGTTCTACTAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATc           <  2:2764006/93‑1 (MQ=255)

                                                                               gTTCTACGAGAGCTTGTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          <  2:8732849/93‑1 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATTCCGGTGGTTCCAGCTCCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          <  1:8449503/93‑1 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          <  1:108106/93‑1 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          >  1:7267791/1‑93 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          >  1:8194720/1‑93 (MQ=255)

                                                                                    acaagagCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCa     <  2:4059327/90‑1 (MQ=255)

                                                                                     cGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAg    <  2:4473592/93‑1 (MQ=255)

                                                                                       agagCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTa  >  2:6622799/1‑93 (MQ=255)

                                                                                       agagCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTa  <  1:9747586/93‑1 (MQ=255)

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

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