BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,005	T→G	94.7%	D73A (GAT→GCT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,005	0	T	G	94.7%	722.5 / 19.0	204	D73A (GAT→GCT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (3/8); new base G (93/100); total (96/108)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

CTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATC  >  NC_000913/4116859‑4117150

ctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCg                                                                                                                                                 <  1:67696/149‑1 (MQ=255)

    aCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAAc                                                                                                                                             <  2:1319890/149‑1 (MQ=255)

        tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGat                                                                                                                                         >  2:986964/1‑149 (MQ=255)

        tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGat                                                                                                                                         >  3:902711/1‑149 (MQ=255)

          aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGatat                                                                                                                                       >  2:1631652/1‑149 (MQ=255)

            gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCg                                                                                                                                     <  3:714120/149‑1 (MQ=255)

              tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGc                                                                                                                                   <  3:1396097/149‑1 (MQ=255)

                  aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTc                                                                                                                               <  2:894546/149‑1 (MQ=255)

                       tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAg                                                                                                                          >  4:153738/1‑149 (MQ=255)

                       tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAg                                                                                                                          >  3:457642/1‑149 (MQ=255)

                       tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAg                                                                                                                          >  1:899084/1‑149 (MQ=255)

                         gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTGCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGAGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                        >  2:1250901/1‑149 (MQ=255)

                               ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTa                                                                                                                  <  1:1377675/149‑1 (MQ=255)

                               ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTa                                                                                                                  <  2:1668751/149‑1 (MQ=255)

                                      gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCCCTTCTACCAGCGt                                                                                                           >  1:70821/1‑149 (MQ=255)

                                      gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                                                                                           <  1:582865/149‑1 (MQ=255)

                                         cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                                                                                        <  1:1063547/149‑1 (MQ=255)

                                          gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                  <  1:389011/106‑1 (MQ=255)

                                          gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                                                                                       >  1:1678678/1‑149 (MQ=255)

                                           acgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                                                                                      >  1:881733/1‑149 (MQ=255)

                                             ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                                                                                                    >  3:854388/1‑149 (MQ=255)

                                              gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                                                                                   <  4:744626/149‑1 (MQ=255)

                                               cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                                                                                                  <  4:1121696/149‑1 (MQ=255)

                                                aatGCCAGACACTGCCGGTATAGAGAGGCTTGCCGGTTTCTTTTACCCAGACAATAGTGCTTTCGCGCTGGTTCGTAATACCGAGAGCTGCAATTTGAGCGGAACGGATATCAGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  2:1069431/147‑1 (MQ=255)

                                                aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  1:156343/149‑1 (MQ=255)

                                                aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  4:1600526/149‑1 (MQ=255)

                                                aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  2:971438/149‑1 (MQ=255)

                                                aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  4:912891/149‑1 (MQ=255)

                                                aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                 <  3:497842/149‑1 (MQ=255)

                                                 cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                                <  4:208775/149‑1 (MQ=255)

                                                   gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                                                                                              >  4:274064/1‑149 (MQ=255)

                                                    ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                             <  1:1515259/149‑1 (MQ=255)

                                                     cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTc                                                                                                                    <  1:786114/125‑1 (MQ=255)

                                                      aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                           >  4:441517/1‑149 (MQ=255)

                                                       gACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                                                                                          >  2:100775/1‑149 (MQ=255)

                                                          aaTGGCGTTATAGATAGGCTTGCCGGTTTTTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                                       >  3:1477794/1‑149 (MQ=255)

                                                                gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                 <  1:1151816/149‑1 (MQ=255)

                                                                 ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATTGGCTCTCGc                                                                                                                             >  2:1729072/1‑104 (MQ=255)

                                                                 ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                                                <  3:274064/149‑1 (MQ=255)

                                                                 ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                                                >  2:230878/1‑149 (MQ=255)

                                                                    tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                             <  1:840750/149‑1 (MQ=255)

                                                                     agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                                                                            >  1:1613446/1‑149 (MQ=255)

                                                                        tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCg                                                                         >  2:294226/1‑149 (MQ=255)

                                                                               gCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                                                  <  1:283221/149‑1 (MQ=255)

                                                                                ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAGTTTGAGAGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                 <  1:1250901/149‑1 (MQ=255)

                                                                                ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                 <  4:1319348/149‑1 (MQ=255)

                                                                                ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                 <  4:1581187/149‑1 (MQ=255)

                                                                                ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACAGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                 <  2:306788/149‑1 (MQ=255)

                                                                                   gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCa                                                              >  2:1083840/1‑149 (MQ=255)

                                                                                        tttttCCCAGTCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         <  4:1477794/149‑1 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         >  4:790842/1‑149 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                                       >  1:1071406/1‑119 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                               >  3:1190186/1‑127 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                               <  4:1190186/127‑1 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         >  3:174652/1‑149 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         >  4:313400/1‑149 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         >  4:294403/1‑149 (MQ=255)

                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                         >  1:74790/1‑149 (MQ=255)

                                                                                           ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                      <  4:329792/149‑1 (MQ=255)

                                                                                            tCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGttt                                                     >  1:641279/1‑149 (MQ=255)

                                                                                               cataCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                  <  4:848595/146‑1 (MQ=255)

                                                                                               cagtcagtAGTGGTTTCACGCTTGTTCGTATTACCGATTGCTGCAATTTGTGTTGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGTTTACCGTTGGGTCGTGTTCTATACAACCTGGTTTTgg                                                  <  2:845062/149‑1 (MQ=255)

                                                                                               cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                                <  4:1052517/103‑1 (MQ=255)

                                                                                                aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                 <  3:807669/149‑1 (MQ=255)

                                                                                                aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                 <  2:667232/149‑1 (MQ=255)

                                                                                                 gACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGt                                                >  3:1007536/1‑149 (MQ=255)

                                                                                                    aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                                                                                              <  3:1038471/100‑1 (MQ=255)

                                                                                                    aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                                             >  2:208411/1‑149 (MQ=255)

                                                                                                    aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                                             <  4:1007536/149‑1 (MQ=255)

                                                                                                     aTAGTGGTTTCGCGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                            >  3:1513719/1‑149 (MQ=255)

                                                                                                     aTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                            <  1:237973/149‑1 (MQ=255)

                                                                                                      tAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                  <  4:368011/142‑1 (MQ=255)

                                                                                                      tAGTGGCTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                  >  3:368011/1‑142 (MQ=255)

                                                                                                         tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                        <  2:522627/149‑1 (MQ=255)

                                                                                                         tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                        >  1:949895/1‑149 (MQ=255)

                                                                                                           gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                 >  2:620589/1‑74 (MQ=255)

                                                                                                           gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                      <  3:59007/149‑1 (MQ=255)

                                                                                                               cACGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAt                                  <  3:1549/149‑1 (MQ=255)

                                                                                                                aCGCTGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                                 >  3:209506/1‑149 (MQ=255)

                                                                                                                     ggTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCgcgc                            <  2:74790/149‑1 (MQ=255)

                                                                                                                        tCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                         >  3:293553/1‑149 (MQ=255)

                                                                                                                          gTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                           <  4:1095883/97‑1 (MQ=255)

                                                                                                                          gTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                           >  3:1095883/1‑97 (MQ=255)

                                                                                                                             aTACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                    >  1:1188994/1‑149 (MQ=255)

                                                                                                                             aTACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                    >  2:1500371/1‑149 (MQ=255)

                                                                                                                                ccGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAtt                                           >  4:504016/1‑123 (MQ=255)

                                                                                                                                ccGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt                 >  4:1713107/1‑149 (MQ=255)

                                                                                                                                 cGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa                >  3:926639/1‑149 (MQ=255)

                                                                                                                                  gATAGCTGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGAAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAGCCTGGTTTTGGGTAGATTTGCTCACATTCGCGCTGCGACACGCTaa               <  4:425651/149‑1 (MQ=255)

                                                                                                                                        tGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                    >  1:851913/1‑106 (MQ=255)

                                                                                                                                        tGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                    <  2:851913/106‑1 (MQ=255)

                                                                                                                                        tGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                         >  3:321179/1‑133 (MQ=255)

                                                                                                                                        tGCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat         >  4:1388725/1‑149 (MQ=255)

                                                                                                                                         gCAATTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt        <  3:7667/149‑1 (MQ=255)

                                                                                                                                          cAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                        >  1:300412/1‑116 (MQ=255)

                                                                                                                                           aaTTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTgg      >  3:674134/1‑149 (MQ=255)

                                                                                                                                           aaTTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTgg      <  4:651855/149‑1 (MQ=255)

                                                                                                                                             tttGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                      <  3:541435/131‑1 (MQ=255)

                                                                                                                                             tttGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                      >  4:541435/1‑131 (MQ=255)

                                                                                                                                             tttGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa    >  3:601825/1‑149 (MQ=255)

CTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATC  >  NC_000913/4116859‑4117150

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