BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,735,666	T→C	S181P (TCC→CCC)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,735,666	0	T	C	99.7%	1250.4 / ‑4.9	368	S181P (TCC→CCC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (202/165); total (203/165)
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.76e-01

Only 100 of 408 total aligned reads displayed.

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAC  >  NC_000913/3735576‑3735755

aTTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACAccc                                                                                         <  2:6176724/93‑1 (MQ=255)

aTTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACAccc                                                                                         >  2:5478441/1‑93 (MQ=255)

       tAGGGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCtt                                                                                  <  1:2440834/93‑1 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                >  1:151557/1‑93 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                <  1:3687552/93‑1 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                >  2:4570785/1‑93 (MQ=255)

           gaagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGa                                                                              <  2:4084283/93‑1 (MQ=255)

           gaagaaCGGTCAAATAATGATACAACCCGCTACCACCAGGAACGGGTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGa                                                                              >  2:7680614/1‑93 (MQ=255)

            acgaaCGGTCAAATAATGTTCCATCCCGCTGCCCCCAGAAACTGTTTTCTCCTGCTGGAACGGGCAGAATCGGACACACCCGCCCCTTTCgtc                                                                             <  2:3554740/91‑3 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGAc                                                                             >  1:8081997/1‑93 (MQ=255)

             agaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACCCACCCGCCCCTTTCGACt                                                                            <  2:7227562/93‑1 (MQ=255)

             agaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACt                                                                            >  2:5511584/1‑93 (MQ=255)

             agaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACt                                                                            >  1:2272969/1‑93 (MQ=255)

              gaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACACGCCCCTTTCGACTg                                                                           >  2:5859062/1‑93 (MQ=255)

                 cGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCa                                                                        >  2:6090642/1‑93 (MQ=255)

                  ggTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCAt                                                                       >  1:3851108/1‑93 (MQ=255)

                  cgTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCAt                                                                       <  1:1885378/92‑1 (MQ=255)

                   gTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATa                                                                      >  1:5712247/1‑93 (MQ=255)

                   gTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATa                                                                      >  1:167935/1‑93 (MQ=255)

                   gTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATa                                                                      >  1:4567335/1‑93 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     <  2:8497833/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     >  1:5573652/1‑93 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     <  1:7714990/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     >  1:3355672/1‑93 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     <  1:3249514/93‑1 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                    >  1:3117217/1‑93 (MQ=255)

                      aaaTAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACg                                                                   <  2:1612195/93‑1 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                  >  1:1045684/1‑93 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                  >  2:4479614/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGTTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                <  2:8330252/93‑1 (MQ=255)

                             gTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGt                                                            >  2:1321446/1‑93 (MQ=255)

                              ttACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGtt                                                           >  1:709550/1‑93 (MQ=255)

                              ttACAACCCGCGGCCACCAGCAACGGTTTACTCCTGATGGAACGGGCAGAAACTGACACAACCGCCCCTTTCGACTGGCATAACGCCACGGtt                                                           >  1:4429407/1‑93 (MQ=255)

                                aCAACCCGCTGCCACCAGCAACGGTTTCCGCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTa                                                         <  1:5235014/93‑1 (MQ=255)

                                 cAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTAc                                                        <  1:4098715/93‑1 (MQ=255)

                                  caCCCGCTGCTACCACAAACGTTGAACCCCGGCAGTACCGGGCACAAACTGACACACCCGCCCCTTTCGCCTGGCATATCGCCACGGTTTACt                                                       <  2:7919986/92‑1 (MQ=255)

                                  aaCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACt                                                       >  2:3906625/1‑93 (MQ=255)

                                   aCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACtt                                                      <  2:510255/93‑1 (MQ=255)

                                    cccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACttt                                                     <  2:3715929/93‑1 (MQ=255)

                                     ccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtg                                                    >  1:8514919/1‑93 (MQ=255)

                                     ccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtg                                                    <  2:7329230/93‑1 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                   <  1:8682056/93‑1 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                   >  1:1206421/1‑93 (MQ=255)

                                       gCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgtg                                                  >  1:1112949/1‑93 (MQ=255)

                                       ccTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgtg                                                  <  1:6577231/92‑1 (MQ=255)

                                        cTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGc                                                 <  1:6496977/93‑1 (MQ=255)

                                        cTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGc                                                 >  1:2855747/1‑93 (MQ=255)

                                         tGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCt                                                <  2:5759620/93‑1 (MQ=255)

                                            caccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACTCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGAc                                             <  1:6758608/93‑1 (MQ=255)

                                             cccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACa                                            <  2:2397026/92‑1 (MQ=255)

                                             accaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCTACTGGCATAACGCCACGGTTTACTTTGTGCTGACa                                            >  2:220322/1‑93 (MQ=255)

                                              ccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAg                                           >  1:6003957/1‑93 (MQ=255)

                                               tagtaaCGTTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCTCTTTCGACTGGCATAACGCCACGGTTTTCTCTGTGcagacaga                                          >  2:7825162/5‑93 (MQ=255)

                                               caGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGa                                          >  1:5825109/1‑93 (MQ=255)

                                               caGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGa                                          <  2:4149489/93‑1 (MQ=255)

                                                gggaaCGGTTTGCTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         <  2:8664145/90‑1 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         <  2:1445025/93‑1 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         >  2:2049561/1‑93 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         >  1:2760041/1‑93 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAtg                                        >  1:1990717/1‑92 (MQ=255)

                                                    aCGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGtt                                     <  2:3061108/93‑1 (MQ=255)

                                                     cGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGttt                                    <  1:6092594/93‑1 (MQ=255)

                                                      ggTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTc                                   >  1:135165/1‑93 (MQ=255)

                                                       gTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCg                                  >  2:1811667/1‑93 (MQ=255)

                                                       gTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCg                                  <  2:6775879/93‑1 (MQ=255)

                                                             tCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                            >  2:6581125/1‑93 (MQ=255)

                                                             tCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                            >  2:1235629/1‑93 (MQ=255)

                                                             tCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                            >  1:348129/1‑93 (MQ=255)

                                                              ccTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACgg                           <  1:7277790/93‑1 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                          >  1:4744083/1‑93 (MQ=255)

                                                                tgctgGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCg                         >  2:8042502/1‑93 (MQ=255)

                                                                tgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCg                         >  2:2547434/1‑93 (MQ=255)

                                                                  ctgGGACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGAt                       <  1:6713542/93‑1 (MQ=255)

                                                                   tgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATc                      <  2:4367959/93‑1 (MQ=255)

                                                                   tgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATc                      >  2:4912793/1‑93 (MQ=255)

                                                                    ggAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATcc                     >  2:3316185/1‑93 (MQ=255)

                                                                     gagcGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATccc                    <  1:835402/90‑1 (MQ=255)

                                                                      aaCGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCa                   <  1:4899953/93‑1 (MQ=255)

                                                                       aCGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAg                  <  1:752267/93‑1 (MQ=255)

                                                                        cGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGt                 <  1:4303111/93‑1 (MQ=255)

                                                                        cGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGt                 >  1:5104601/1‑93 (MQ=255)

                                                                        agggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGt                 <  2:6852179/92‑1 (MQ=255)

                                                                         gggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTa                >  1:2561619/1‑93 (MQ=255)

                                                                         gggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTa                >  1:4995076/1‑93 (MQ=255)

                                                                         gggCAGAAACGCACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTCCTTTGTGCTGACAGATAGTTTCGAAAACGGCGATCCCAGTa                <  2:6273607/93‑1 (MQ=255)

                                                                          ggCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTaa               <  1:2612703/93‑1 (MQ=255)

                                                                           gCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAAt              <  2:8757120/93‑1 (MQ=255)

                                                                            cAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATg             >  1:6360197/1‑93 (MQ=255)

                                                                             aGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGa            <  2:6946677/93‑1 (MQ=255)

                                                                             aGAAACTCACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAGCGGTTGGAAAACGGCGATCCCAGTAATgc            >  2:2759007/1‑92 (MQ=255)

                                                                              gAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGAc           >  2:710330/1‑93 (MQ=255)

                                                                                aaCTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCa         <  2:5193728/93‑1 (MQ=255)

                                                                                aaCTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCa         >  1:2888703/1‑93 (MQ=255)

                                                                                 aCTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCag        >  2:2921119/1‑93 (MQ=255)

                                                                                  cTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACATTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCaga       >  2:7446946/1‑93 (MQ=255)

                                                                                   tGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCagag      >  1:5368472/1‑93 (MQ=255)

                                                                                    gACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGt     <  2:1712387/93‑1 (MQ=255)

                                                                                     acacacCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGtt    <  1:6912778/93‑1 (MQ=255)

                                                                                      cccctCCGCCCCTCTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTa   <  2:5914333/89‑1 (MQ=255)

                                                                                       acacCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAc  <  2:3826493/93‑1 (MQ=255)

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAC  >  NC_000913/3735576‑3735755

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