BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,092,774:1	+A	11.5%	coding (897/951 nt)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,092,774	1	.	A	11.5%	465.7 / 67.3	131	coding (897/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand): ref base . (62/54); new base A (8/7); total (70/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 132 total aligned reads displayed.

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092631‑3092920

  gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCa                                                                                                                                             >  4:202805/1‑149 (MQ=255)

    aaaTCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCAcc                                                                                                                                           <  1:424587/149‑1 (MQ=255)

         gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                      <  3:747234/149‑2 (MQ=255)

         gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                      <  2:273980/149‑2 (MQ=255)

         gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                      <  1:734891/149‑2 (MQ=255)

              tCAGATCGTGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                 <  3:100381/149‑1 (MQ=255)

              tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                 <  2:151672/149‑1 (MQ=255)

              tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                 <  2:792023/149‑1 (MQ=255)

               cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                <  4:73748/149‑1 (MQ=255)

                  aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGa                                                                                                                             >  1:589236/1‑149 (MQ=255)

                   tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGAt                                                                                                                            >  1:854556/1‑149 (MQ=255)

                   tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGAt                                                                                                                            >  3:473336/1‑149 (MQ=255)

                     gggCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGc                                                                                                                          <  3:879355/149‑1 (MQ=255)

                      ggCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGTCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGcc                                                                                                                         >  1:760635/1‑149 (MQ=255)

                            cGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGaa                                                                                                                   <  1:631310/149‑1 (MQ=255)

                                 aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCgg                                                                                                                                         <  4:831776/122‑1 (MQ=255)

                                 aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg                                                                                                              >  1:160055/1‑149 (MQ=255)

                                  aagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGt                                                                                                             >  1:311207/1‑149 (MQ=255)

                                   agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                            <  3:44969/149‑1 (MQ=255)

                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCGCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGAATGCCATCGAAGCACGTTTAc                                                                                                         <  4:343168/149‑1 (MQ=255)

                                          gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                     <  2:742291/149‑1 (MQ=255)

                                             cTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                 <  1:862010/118‑1 (MQ=255)

                                             cTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                 >  2:862010/1‑118 (MQ=255)

                                              tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                 >  3:531467/1‑149 (MQ=255)

                                               gATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGt                                                                                                <  4:113552/149‑1 (MQ=255)

                                               gATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGt                                                                                                >  2:405476/1‑149 (MQ=255)

                                               gATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGt                                                                                                >  1:907905/1‑149 (MQ=255)

                                                 tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAGGCAGAGTTTCCGGTG‑TCGATCACGGGAATGTTAATGGATGCCATCGAAGCACGGTTACAGCAGCAGTaa                                                                                              >  1:205386/1‑149 (MQ=255)

                                                  ttatgttgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAc                                                                                             <  2:723463/146‑1 (MQ=255)

                                                     tgttgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCa                                                                                          <  1:728803/149‑1 (MQ=255)

                                                      gttgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCAc                                                                                         >  4:115073/1‑149 (MQ=255)

                                                              ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCga                                                                                 <  4:102964/149‑1 (MQ=255)

                                                               gATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCgag                                                                                >  2:26772/1‑149 (MQ=255)

                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCAcc                                                                                        <  2:862758/139‑1 (MQ=255)

                                                                  atcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAg                                                                             <  2:627894/149‑1 (MQ=255)

                                                                   tcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAgg                                                                            >  3:219613/1‑149 (MQ=255)

                                                                   tcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAgg                                                                            >  3:440581/1‑149 (MQ=255)

                                                                     atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                       >  3:310059/1‑120 (MQ=255)

                                                                     atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                       <  4:310059/120‑1 (MQ=255)

                                                                     atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                          >  3:904222/1‑149 (MQ=255)

                                                                        ggCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                       <  4:28266/149‑1 (MQ=255)

                                                                          cGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                  >  3:339541/1‑120 (MQ=255)

                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  3:287322/149‑1 (MQ=255)

                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  2:311207/149‑1 (MQ=255)

                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  2:830390/149‑1 (MQ=255)

                                                                             ccGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                  <  3:124326/149‑1 (MQ=255)

                                                                               gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                >  1:15273/1‑149 (MQ=255)

                                                                                tCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                               <  4:440581/149‑1 (MQ=255)

                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                          >  2:435108/1‑149 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                          >  4:231939/1‑148 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                          <  3:231939/148‑1 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         >  2:389253/1‑149 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         >  4:319052/1‑149 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         >  1:769780/1‑149 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         >  2:295158/1‑149 (MQ=255)

                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         <  2:371565/149‑1 (MQ=255)

                                                                                          aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAAt                                                                                                                                      >  1:860232/1‑68 (MQ=255)

                                                                                           ttAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACa                                                    >  4:629821/1‑149 (MQ=255)

                                                                                           ttAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACa                                                    >  4:823721/1‑149 (MQ=255)

                                                                                           ttAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACa                                                    >  4:850373/1‑149 (MQ=255)

                                                                                                gTCACCAGCCCAACCTGTATTCGTGAGATTGACGCAGAGTTTCCGGTG‑TCGAACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGcc                                               >  4:159436/1‑149 (MQ=255)

                                                                                                   aCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                          >  4:258549/1‑119 (MQ=255)

                                                                                                   aCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCctt                                            <  3:115073/149‑1 (MQ=255)

                                                                                                   aCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCctt                                            >  2:297516/1‑149 (MQ=255)

                                                                                                         ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                      >  4:441657/1‑149 (MQ=255)

                                                                                                         ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                      <  4:473336/149‑1 (MQ=255)

                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCa                                     >  3:342166/1‑149 (MQ=255)

                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCa                                     >  2:343115/1‑149 (MQ=255)

                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                 >  1:664712/1‑149 (MQ=255)

                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                 >  3:19332/1‑149 (MQ=255)

                                                                                                                    tCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTgg                           <  3:179625/149‑1 (MQ=255)

                                                                                                                      gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                         >  3:584121/1‑149 (MQ=255)

                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTGACCCACCGTAGAGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcg                        >  2:667975/1‑149 (MQ=255)

                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                  <  4:368516/107‑1 (MQ=255)

                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                  >  3:368516/1‑107 (MQ=255)

                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcg                        >  3:280596/1‑149 (MQ=255)

                                                                                                                        gagaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcgc                       <  2:603466/149‑1 (MQ=255)

                                                                                                                         agaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCt                      >  1:425765/1‑149 (MQ=255)

                                                                                                                           aTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGt                    >  3:76482/1‑149 (MQ=255)

                                                                                                                           aTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGt                    >  1:886698/1‑149 (MQ=255)

                                                                                                                             tGAAGCAGAGTTTCCGGTG‑TCGAGCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAACCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTGTCCACGCATACTGGGCGCTGTTg                  <  4:425470/149‑1 (MQ=255)

                                                                                                                                 gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttt              >  4:238495/1‑149 (MQ=255)

                                                                                                                                   agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCTTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            <  1:578552/149‑1 (MQ=255)

                                                                                                                                   agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            <  3:240147/149‑1 (MQ=255)

                                                                                                                                   agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATACCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            <  2:464550/149‑1 (MQ=255)

                                                                                                                                    gagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt           <  1:229584/149‑1 (MQ=255)

                                                                                                                                     agTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg          >  2:844057/1‑149 (MQ=255)

                                                                                                                                      cattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttt             <  1:149718/143‑1 (MQ=255)

                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt           >  4:32078/1‑146 (MQ=255)

                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  3:69961/1‑149 (MQ=255)

                                                                                                                                       attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCAt                                >  2:835531/2‑125 (MQ=255)

                                                                                                                                       attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtt               <  2:26427/141‑1 (MQ=255)

                                                                                                                                       attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtt               >  1:26427/2‑142 (MQ=255)

                                                                                                                                       attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  1:530802/2‑149 (MQ=255)

                                                                                                                                       attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAACACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  3:565842/2‑149 (MQ=255)

                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc       >  3:832830/1‑149 (MQ=255)

                                                                                                                                          ccGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa     <  4:342166/149‑1 (MQ=255)

                                                                                                                                          ccGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa     <  1:909942/149‑1 (MQ=255)

                                                                                                                                          ccGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa     <  1:26772/149‑1 (MQ=255)

                                                                                                                                             gtg‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGa  >  4:280230/1‑149 (MQ=255)

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092631‑3092920

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