BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,703,598	C→T	7.8%	R237R (CGG→CGA)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,598	0	C	T	7.8%	484.1 / 20.7	152	R237R (CGG→CGA)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (73/67); new base T (6/6); total (79/73)
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.36e-01

Only 100 of 165 total aligned reads displayed.

GTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAA  >  NC_000913/1703450‑1703746

    cTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTg                                                                                                                                                  >  3:1126460/1‑149 (MQ=255)

     tCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGc                                                                                                                                                 >  1:1437335/1‑149 (MQ=255)

       aCCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcat                                                                                                                                               <  4:567596/149‑3 (MQ=255)

       aCCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTg                                                                                                                                               <  2:506593/149‑1 (MQ=255)

         cTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATAc                                                                                                                                                      <  3:669064/140‑1 (MQ=255)

          ttCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatgga                                                                                                                                            <  1:319394/149‑6 (MQ=255)

                  gggTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCt                                                                                                                                    <  1:1562728/149‑1 (MQ=255)

                    gTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGc                                                                                                                                  >  2:1879951/1‑149 (MQ=255)

                    gTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGc                                                                                                                                  <  3:1660970/149‑1 (MQ=255)

                      aaGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCa                                                                                                                                <  1:158047/149‑1 (MQ=255)

                            gCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATgg                                                                                                                          <  1:230358/149‑1 (MQ=255)

                              cATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTa                                                                                                                        >  1:1540157/1‑149 (MQ=255)

                                   cGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTa                                                                                                                        <  1:777245/144‑1 (MQ=255)

                                    gCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACTTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                                                                              <  2:1402120/121‑1 (MQ=255)

                                    gCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTG‑CACGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTa                                                                                                                 <  2:1365344/149‑1 (MQ=255)

                                        atCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTaaaa                                                                                                              <  4:478955/149‑1 (MQ=255)

                                            tAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGAc                                                                                                          >  3:1665972/1‑149 (MQ=255)

                                                 ccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgcc                                                                                                     >  2:1437787/1‑149 (MQ=255)

                                                   aGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccgc                                                                                                   >  1:1207960/1‑149 (MQ=255)

                                                     tCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtg                                                                                                 <  1:713188/149‑1 (MQ=255)

                                                     tCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtg                                                                                                 >  1:1308159/1‑149 (MQ=255)

                                                       ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATAcc                                                                                                                                                     <  2:763555/95‑1 (MQ=255)

                                                        tCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTgg                                                                                              <  2:590867/149‑1 (MQ=255)

                                                         ccTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGTCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGa                                                                                             <  4:1859616/149‑1 (MQ=255)

                                                           tGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATa                                                                                           <  4:1665972/149‑1 (MQ=255)

                                                            gCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAg                                                                                          >  3:932/1‑149 (MQ=255)

                                                              ggTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGa                                                                                        <  2:1207960/149‑1 (MQ=255)

                                                              ggTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGa                                                                                        <  4:1844319/149‑1 (MQ=255)

                                                               gTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                                                                        >  3:1690543/1‑132 (MQ=255)

                                                               gTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                                                                        <  4:1690543/132‑1 (MQ=255)

                                                                tAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                           >  1:1245625/1‑128 (MQ=255)

                                                                tAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                           <  2:1245625/128‑1 (MQ=255)

                                                                     gCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                           <  2:1036482/123‑1 (MQ=255)

                                                                     gCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                           >  1:1036482/1‑123 (MQ=255)

                                                                     gCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 >  3:1075615/1‑149 (MQ=255)

                                                                       tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                                                                   >  2:1731200/1‑145 (MQ=255)

                                                                       tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                                                                   <  1:1731200/145‑1 (MQ=255)

                                                                       tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                               >  3:1418961/1‑149 (MQ=255)

                                                                       tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGAACCATGCACGATATTCCGTGCTGACTAAGCGGCTGCCAGCATGGTACCGTGGAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                               >  1:1346374/1‑149 (MQ=255)

                                                                             tttttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 <  3:1466397/141‑1 (MQ=255)

                                                                             tttttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 >  4:1466397/1‑141 (MQ=255)

                                                                                    gACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg                                                                  >  3:1037989/1‑149 (MQ=255)

                                                                                    gACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg                                                                  >  4:416414/1‑149 (MQ=255)

                                                                                            tcCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTg                                                          >  1:211613/1‑149 (MQ=255)

                                                                                                gtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCt                                                      >  1:318145/1‑149 (MQ=255)

                                                                                                gtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATAGTACAGTGTAAAATGACTCGCCTCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCt                                                      >  3:1376945/1‑149 (MQ=255)

                                                                                                  ggATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTcc                                                    <  1:492531/149‑1 (MQ=255)

                                                                                                   gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCg                                                   >  4:718241/1‑149 (MQ=255)

                                                                                                    aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                                  >  1:54253/1‑149 (MQ=255)

                                                                                                    aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                                  >  3:291004/1‑149 (MQ=255)

                                                                                                    aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                                  >  3:1186574/1‑149 (MQ=255)

                                                                                                         ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATa                                                                                           <  2:376734/103‑1 (MQ=255)

                                                                                                         ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATa                                                                                           >  1:376734/1‑103 (MQ=255)

                                                                                                         ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                             >  2:1334190/1‑149 (MQ=255)

                                                                                                         ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                             >  3:94598/1‑149 (MQ=255)

                                                                                                         ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                             <  1:737058/149‑1 (MQ=255)

                                                                                                              tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                               >  3:493747/1‑110 (MQ=255)

                                                                                                              tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                               <  2:463999/110‑1 (MQ=255)

                                                                                                               atTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                       <  4:291004/149‑1 (MQ=255)

                                                                                                                tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTc                                                               >  1:1031626/1‑124 (MQ=255)

                                                                                                                tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                      >  3:1208643/1‑149 (MQ=255)

                                                                                                                tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                      >  4:1377517/1‑149 (MQ=255)

                                                                                                                tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                      >  3:1564005/1‑149 (MQ=255)

                                                                                                                 ttCATATACCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGGGc                                     <  1:1205067/149‑1 (MQ=255)

                                                                                                                 ttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGc                                     >  4:1240842/1‑149 (MQ=255)

                                                                                                                    acaTGACTGCCGCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCaa                                  >  2:222345/1‑149 (MQ=255)

                                                                                                                    acaTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATa                                                                                           >  4:726162/1‑92 (MQ=255)

                                                                                                                     caTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCaaa                                 <  2:318145/149‑1 (MQ=255)

                                                                                                                      aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAAc                                          <  2:166729/139‑1 (MQ=255)

                                                                                                                        aaCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 <  3:444554/98‑1 (MQ=255)

                                                                                                                         aCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGcc                                    <  3:358024/142‑1 (MQ=255)

                                                                                                                         aCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTa                             <  2:287588/149‑1 (MQ=255)

                                                                                                                           tGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTAcc                           <  2:1039999/149‑1 (MQ=255)

                                                                                                                                cGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCAt                      >  2:1416513/1‑149 (MQ=255)

                                                                                                                                  ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 <  3:475067/88‑1 (MQ=255)

                                                                                                                                  ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                                 >  4:475067/1‑88 (MQ=255)

                                                                                                                                    gagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGTATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTg                  >  4:1100649/7‑149 (MQ=255)

                                                                                                                                      gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAc                >  2:211395/5‑149 (MQ=255)

                                                                                                                                       cATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc               <  4:1037989/149‑1 (MQ=255)

                                                                                                                                        aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                       <  2:1578508/124‑1 (MQ=255)

                                                                                                                                        aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCg              >  3:828860/1‑149 (MQ=255)

                                                                                                                                         tGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtg                                        <  3:1577046/122‑1 (MQ=255)

                                                                                                                                         tGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtg                                        >  4:1577046/1‑122 (MQ=255)

                                                                                                                                         agCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCgg             >  2:1854349/2‑149 (MQ=255)

                                                                                                                                          gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGt            <  1:1879951/149‑1 (MQ=255)

                                                                                                                                          gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGt            >  2:995498/1‑149 (MQ=255)

                                                                                                                                           cACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAAGAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTa           <  3:1100649/149‑1 (MQ=255)

                                                                                                                                            aCGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaa          >  4:334394/1‑149 (MQ=255)

                                                                                                                                             cggtataTCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaaa         <  1:211395/146‑1 (MQ=255)

                                                                                                                                              ggtataTCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc        <  1:1854349/147‑1 (MQ=255)

                                                                                                                                              gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc        <  4:1087580/149‑1 (MQ=255)

                                                                                                                                              gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc        <  3:1852041/149‑1 (MQ=255)

                                                                                                                                              gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc        <  4:932/149‑1 (MQ=255)

                                                                                                                                               atatACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAcc       <  3:718241/149‑1 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaaa         >  4:1547097/1‑142 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaaa         >  3:1164116/1‑142 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaaa         <  3:1547097/142‑1 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaaa         <  4:1164116/142‑1 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc        <  4:1353269/143‑1 (MQ=255)

                                                                                                                                                    ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCaaaaa  <  4:1392651/149‑1 (MQ=255)

GTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAA  >  NC_000913/1703450‑1703746

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