BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,212	C→G	5.6%	G124R (GGT→CGT) ‡	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,212	0	C	G	5.6%	612.7 / 18.2	178	G124R (GGT→CGT) ‡	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (83/82); new base G (8/2); total (93/85)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

AGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGG  >  NC_000913/3566064‑3566360

agccagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCAcc                                                                                                                                                      >  2:685236/1‑149 (MQ=255)

          cTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                            <  2:932573/149‑1 (MQ=255)

          cTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                            <  3:173939/149‑1 (MQ=255)

             gCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCt                                                                                                                                         >  4:679805/1‑149 (MQ=255)

             gCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCt                                                                                                                                         >  1:627814/1‑149 (MQ=255)

             gCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCt                                                                                                                                         <  4:1454763/149‑1 (MQ=255)

              cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                        >  2:130602/1‑149 (MQ=255)

               gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTgg                                                                                                                                       <  2:261912/149‑1 (MQ=255)

               gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTgg                                                                                                                                       <  4:56486/149‑1 (MQ=255)

               gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTgg                                                                                                                                       >  3:250150/1‑149 (MQ=255)

                  tctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTc                                                                                                                                    >  3:442536/1‑149 (MQ=255)

                    tttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAt                                                                                                                                  <  3:674321/149‑1 (MQ=255)

                        tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                              <  3:679805/149‑1 (MQ=255)

                        tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                              >  4:315521/1‑149 (MQ=255)

                            tACCNTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                          >  3:836472/1‑149 (MQ=255)

                                gTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaa                                                                                                                      >  2:1182101/1‑149 (MQ=255)

                                 ttAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaat                                                                                                                     >  3:388825/1‑149 (MQ=255)

                                    aCGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAg                                                                                                                  >  2:907031/1‑149 (MQ=255)

                                      gATGTTCTGCTTCAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGct                                                                                                                >  1:980698/1‑149 (MQ=255)

                                          ttCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                                            >  2:262763/1‑149 (MQ=255)

                                          ttCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                                            >  1:771822/1‑149 (MQ=255)

                                           tCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGAAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                           >  1:612050/1‑149 (MQ=255)

                                            cTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCga                                                                                                          <  4:423894/149‑1 (MQ=255)

                                             tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttctt                                                                                                                            >  1:1065540/1‑130 (MQ=255)

                                             tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttctt                                                                                                                            <  2:1065540/130‑1 (MQ=255)

                                                ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATCTTCTTCatcaat                                                                                                                     >  2:646434/1‑134 (MQ=255)

                                                 tGAACATACCGTAGTCATAGCGGATGCCGTAACCACGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttctt                                                                                                                            <  4:1241135/126‑1 (MQ=255)

                                                 tGAACATACCGTAGTCATAGCGGATGCCGTAACAACGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttctt                                                                                                                            >  3:1241135/1‑126 (MQ=255)

                                                  gAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTa                                                                                                    >  1:199109/1‑149 (MQ=255)

                                                   aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTaa                                                                                                   <  1:1051433/149‑1 (MQ=255)

                                                    aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAc                                                                                                  <  3:231928/149‑1 (MQ=255)

                                                         ccGTAGTCATAGCGGATGCCGTAACCTCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTCCTTAATCAATCAGCTCTTCGAGATTTAACCCCAt                                                                                             >  4:58235/1‑149 (MQ=255)

                                                          cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAGCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                              <  4:1049562/115‑1 (MQ=255)

                                                          cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                              >  3:1049562/1‑115 (MQ=255)

                                                            tAGTAATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCACAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCATCGAGATTTAACCCCATCGc                                                                                          >  2:1345839/1‑149 (MQ=255)

                                                             aGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCt                                                                                         <  2:731016/149‑1 (MQ=255)

                                                               tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTc                                                                                       <  2:1289250/149‑1 (MQ=255)

                                                                 aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                     >  1:1257670/1‑149 (MQ=255)

                                                                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAg                                                                                    >  1:1104821/1‑149 (MQ=255)

                                                                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAg                                                                                    >  1:798175/1‑149 (MQ=255)

                                                                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAg                                                                                    >  2:488811/1‑149 (MQ=255)

                                                                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAACCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCATCTCTTCGAGATTTAACCCCATCGCTTCCAg                                                                                    >  1:958972/1‑149 (MQ=255)

                                                                   aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGt                                                                                   >  2:1348868/1‑149 (MQ=255)

                                                                    gCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                  >  1:208693/1‑149 (MQ=255)

                                                                      ggATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgcg                                                                                >  3:310229/1‑149 (MQ=255)

                                                                          gCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCt                                                                            >  4:341433/1‑149 (MQ=255)

                                                                           ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTg                                                                           >  1:1004309/1‑149 (MQ=255)

                                                                                 ccGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATc                                                                     <  2:1091743/149‑1 (MQ=255)

                                                                                    cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                                  <  3:566570/149‑1 (MQ=255)

                                                                                    cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                                  <  1:380078/149‑1 (MQ=255)

                                                                                    cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                                  <  3:341433/149‑1 (MQ=255)

                                                                                    cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                                  <  1:1359284/149‑1 (MQ=255)

                                                                                      ccccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatca                                                                                                                   >  4:130782/1‑98 (MQ=255)

                                                                                      ccccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGt                                                                <  4:584261/149‑1 (MQ=255)

                                                                                       cccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTa                                                               >  4:753050/1‑149 (MQ=255)

                                                                                        ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaa                                                              >  4:510674/1‑149 (MQ=255)

                                                                                        ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATAGCTTCCAGTGCGCCCTGTACATCTTCGTaa                                                              >  3:871756/1‑149 (MQ=255)

                                                                                         cGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaaa                                                             <  2:1222710/149‑1 (MQ=255)

                                                                                          ggCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACTGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAt                                                            <  2:783786/149‑1 (MQ=255)

                                                                                            cAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTc                                                          <  2:1378341/149‑1 (MQ=255)

                                                                                             aaCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTcc                                                         >  4:1372409/1‑149 (MQ=255)

                                                                                               cccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTa                                                       <  3:1323030/149‑1 (MQ=255)

                                                                                                ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                      <  1:528278/149‑1 (MQ=255)

                                                                                                ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                      <  2:1472644/149‑1 (MQ=255)

                                                                                                ccTAACGTCGCCAGAGAATCAACGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTTAATTCCTAg                                                      <  1:1270156/149‑1 (MQ=255)

                                                                                                  tAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTGAATTCCTa                                                       >  3:946560/1‑146 (MQ=255)

                                                                                                      gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaa                                                >  1:1289054/1‑149 (MQ=255)

                                                                                                      gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaa                                                >  4:1284752/1‑149 (MQ=255)

                                                                                                          ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                    <  1:150624/125‑1 (MQ=255)

                                                                                                           cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATc                                                                                            >  1:1379182/1‑100 (MQ=255)

                                                                                                           cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCCGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATgg                                           <  2:116190/149‑1 (MQ=255)

                                                                                                             gagaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCg                                         >  1:559358/1‑149 (MQ=255)

                                                                                                               gaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAATCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGtt                                       <  1:685236/149‑1 (MQ=255)

                                                                                                                  tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGga                                    >  4:199300/1‑149 (MQ=255)

                                                                                                                  tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGga                                    >  1:1292658/1‑149 (MQ=255)

                                                                                                                    aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaaa                                                             <  4:166853/122‑1 (MQ=255)

                                                                                                                    aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAACTTCCTCGCGACAACATGGCGTTGgaga                                  <  4:388825/149‑1 (MQ=255)

                                                                                                                     aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGagag                                 >  4:950715/1‑149 (MQ=255)

                                                                                                                      ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGACCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGc                                <  4:871756/149‑1 (MQ=255)

                                                                                                                            aagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAc                                                                                                  <  2:1082133/77‑1 (MQ=255)

                                                                                                                                cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCaa                      <  2:199109/149‑1 (MQ=255)

                                                                                                                                cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCaa                      <  1:1182101/149‑1 (MQ=255)

                                                                                                                                cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCaa                      <  4:263813/149‑1 (MQ=255)

                                                                                                                                 gccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAAt                     <  2:597499/149‑1 (MQ=255)

                                                                                                                                   cAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCa                   >  2:67937/1‑149 (MQ=255)

                                                                                                                                    aaatgacccAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaa                  >  2:103670/7‑149 (MQ=38)

                                                                                                                                       cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaaa               <  3:561036/149‑1 (MQ=255)

                                                                                                                                       cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaaa               <  3:472535/149‑1 (MQ=255)

                                                                                                                                       cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaaa               <  2:958972/149‑1 (MQ=255)

                                                                                                                                           ccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc           >  2:103536/1‑149 (MQ=255)

                                                                                                                                           ccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc           >  1:866042/1‑149 (MQ=255)

                                                                                                                                           ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc           >  2:402066/1‑149 (MQ=255)

                                                                                                                                           ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc           <  3:1286712/149‑1 (MQ=255)

                                                                                                                                            cAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa          >  3:710431/1‑149 (MQ=255)

                                                                                                                                            cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa          <  1:801716/149‑1 (MQ=255)

                                                                                                                                             aGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt         <  3:1256943/149‑1 (MQ=255)

                                                                                                                                                ccACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGa      >  3:178515/1‑149 (MQ=255)

                                                                                                                                                  aCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                  <  2:174568/119‑1 (MQ=255)

                                                                                                                                                  aCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATCTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACa    >  2:971641/1‑149 (MQ=255)

                                                                                                                                                   ccGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAg   >  3:1422570/1‑149 (MQ=255)

AGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGG  >  NC_000913/3566064‑3566360

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