BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,933	G→T	50.8%	P97H (CCT→CAT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,933	0	G	T	50.8%	11.9 / 235.1	128	P97H (CCT→CAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (31/32); new base T (32/33); total (63/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 139 total aligned reads displayed.

GGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGT  >  NC_000913/4116787‑4117080

ggATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAgg                                                                                                                                                    >  1:594374/1‑148 (MQ=255)

ggATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGc                                                                                                                                                   <  3:223168/149‑1 (MQ=255)

       cttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTACATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCg                                                                                                                                             >  2:414500/1‑148 (MQ=255)

       cttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCgg                                                                                                                                            >  3:673107/1‑149 (MQ=255)

        ttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCgg                                                                                                                                            >  1:367571/1‑148 (MQ=255)

           acttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTc                                                                                                                                        >  4:340184/1‑149 (MQ=255)

             tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCtt                                                                                                                                      >  4:100960/1‑149 (MQ=255)

             tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCtt                                                                                                                                      <  4:180646/149‑1 (MQ=255)

              ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCtt                                                                                                                                      >  2:119136/1‑148 (MQ=255)

                ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCttttt                                                                                                                                   <  3:721775/149‑1 (MQ=255)

                  tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTcc                                                                                                                                 <  3:490037/149‑1 (MQ=255)

                  tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTc                                                                                                                                  <  2:594374/148‑1 (MQ=255)

                         aaaGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACa                                                                                                                           >  1:623770/1‑148 (MQ=255)

                           aGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATa                                                                                                                        <  4:359995/149‑1 (MQ=255)

                             tACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAg                                                                                                                       >  2:924570/1‑148 (MQ=255)

                             tACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGCTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTGTCTTTTTCCCAGACAATAGt                                                                                                                      >  3:529768/1‑149 (MQ=255)

                                gggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCACACAAGAGTgg                                                                                                                    <  2:541777/148‑1 (MQ=255)

                                 ggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                   >  1:371008/1‑148 (MQ=255)

                                     aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCa                                                                                                               <  2:906339/148‑1 (MQ=255)

                                       tCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACg                                                                                                             >  1:518501/1‑148 (MQ=255)

                                         aCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCt                                                                                                           >  2:775088/1‑148 (MQ=255)

                                          ccAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                                                                                                         >  3:384803/1‑149 (MQ=255)

                                          ccAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCt                                                                                                           <  3:397401/147‑1 (MQ=255)

                                           cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                                                                                                        <  4:816089/149‑1 (MQ=255)

                                           cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                                                                                                         >  1:185617/1‑148 (MQ=255)

                                             gACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGtt                                                                                                       <  2:823094/148‑1 (MQ=255)

                                               ccGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGTCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                                                                                                    <  3:813500/149‑1 (MQ=255)

                                                cGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                                                                                                    <  1:119136/148‑1 (MQ=255)

                                                 ggTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                                                                                                   >  1:638792/1‑148 (MQ=255)

                                                 ggTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                                                                                                   <  2:518500/148‑1 (MQ=255)

                                                    aTTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                                                                                                <  1:150754/148‑1 (MQ=255)

                                                    aTTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                                                                                               <  3:489181/149‑1 (MQ=255)

                                                        ctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAt                                                                                           >  3:182910/1‑149 (MQ=255)

                                                         tgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAt                                                                                           >  1:759543/1‑148 (MQ=255)

                                                           cGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc                                                                                        <  3:528889/149‑1 (MQ=255)

                                                                ataATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCaa                                                                                   <  3:218753/149‑1 (MQ=255)

                                                                   aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                                >  3:280849/1‑149 (MQ=255)

                                                                   aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                                <  3:340184/149‑1 (MQ=255)

                                                                    tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCTGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                                >  2:810621/1‑148 (MQ=255)

                                                                    tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                               >  4:468410/1‑149 (MQ=255)

                                                                    tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                               >  3:65352/1‑149 (MQ=255)

                                                                      ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                             >  3:52537/1‑149 (MQ=255)

                                                                      ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                             <  4:512755/149‑1 (MQ=255)

                                                                      ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                             <  3:386095/149‑1 (MQ=255)

                                                                      ctctAATCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGtt                                                                                                                                          <  4:10727/87‑1 (MQ=255)

                                                                       tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                            <  4:57947/149‑1 (MQ=255)

                                                                       tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                            <  4:52537/149‑1 (MQ=255)

                                                                         tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACg                                                                           >  1:401824/1‑148 (MQ=255)

                                                                          aaaCCGTCACGTTTTAAATGCTCGCAGCTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGa                                                                         >  4:845706/1‑149 (MQ=255)

                                                                            aCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTGCCGGACACCAGTGGTTTCACGCTGGTTCGTAATACCGATAACTGCAATTTGATCGGaa                                                                        <  2:235572/148‑1 (MQ=255)

                                                                             ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                      >  3:682723/1‑149 (MQ=255)

                                                                                tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                   <  3:100960/149‑1 (MQ=255)

                                                                                 cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGAAAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                   >  1:69406/1‑148 (MQ=255)

                                                                                   cGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                 >  2:305710/1‑148 (MQ=255)

                                                                                    gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                               >  3:204740/1‑149 (MQ=255)

                                                                                     ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                              >  3:614390/1‑149 (MQ=255)

                                                                                      tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                              >  2:931224/1‑148 (MQ=255)

                                                                                      tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCgg                                                              >  2:77769/1‑148 (MQ=255)

                                                                                       ttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGc                                                             <  2:623770/148‑1 (MQ=255)

                                                                                          aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCttt                                                          <  2:47840/148‑1 (MQ=255)

                                                                                          aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                         <  4:384803/149‑1 (MQ=255)

                                                                                           aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                             <  2:455886/128‑1 (MQ=255)

                                                                                           aTGCTCGCAGAGTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                             >  1:455886/1‑128 (MQ=255)

                                                                                            tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGt                                                                                                                      <  2:325839/86‑1 (MQ=255)

                                                                                            tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                        >  1:727850/1‑148 (MQ=255)

                                                                                             gCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCGGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGcc                                                      <  3:862910/149‑1 (MQ=255)

                                                                                               tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAg                                                    >  4:436657/1‑149 (MQ=255)

                                                                                                        tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                           <  4:65735/149‑1 (MQ=255)

                                                                                                            cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATTGGCTTTCGCCAGCACTTCTACCAg                                        >  1:365646/1‑148 (MQ=255)

                                                                                                            cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAg                                        >  1:812308/1‑148 (MQ=255)

                                                                                                              gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTGGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCg                                      >  2:102808/1‑148 (MQ=255)

                                                                                                                      gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                              <  1:151677/148‑1 (MQ=255)

                                                                                                                        aCTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                           >  3:737537/1‑149 (MQ=255)

                                                                                                                        aCTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                            <  2:69406/148‑1 (MQ=255)

                                                                                                                        aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTgg                           >  4:388006/1‑149 (MQ=255)

                                                                                                                        aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCGTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGaa                                                                             >  4:552849/1‑98 (MQ=255)

                                                                                                                        aCTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTgg                           <  4:614390/149‑1 (MQ=255)

                                                                                                                              aGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                      <  2:179588/148‑1 (MQ=255)

                                                                                                                                aCAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                          <  1:256910/142‑1 (MQ=255)

                                                                                                                                aCAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                          >  2:256910/1‑142 (MQ=255)

                                                                                                                                aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                                                                                                <  4:317993/56‑1 (MQ=255)

                                                                                                                                aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                                                                                                >  3:317993/1‑56 (MQ=255)

                                                                                                                                aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                    <  1:305710/148‑1 (MQ=255)

                                                                                                                                  aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCTCAGACAATAGTGGTTTCACGCTGGTTTGTAATTCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCCTGGAGCTTTGGGTGGCCCAGa                 <  4:800561/149‑1 (MQ=255)

                                                                                                                                       cGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc             <  1:766134/148‑1 (MQ=255)

                                                                                                                                        gTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc            <  2:759543/148‑1 (MQ=255)

                                                                                                                                        gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                         >  3:881476/1‑119 (MQ=255)

                                                                                                                                        gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                         <  4:881476/119‑1 (MQ=255)

                                                                                                                                         ttATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt          >  3:426284/1‑149 (MQ=255)

                                                                                                                                          tataGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                 <  2:704068/141‑1 (MQ=255)

                                                                                                                                          tataGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt         <  3:569528/149‑1 (MQ=255)

                                                                                                                                           atagatatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt         <  1:226722/148‑1 (MQ=255)

                                                                                                                                            tagatatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg       <  4:885837/149‑1 (MQ=255)

                                                                                                                                            tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg       <  4:65352/149‑1 (MQ=255)

                                                                                                                                               atatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc                                                                                        <  4:2695/65‑1 (MQ=255)

                                                                                                                                               atatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc                                                                                        >  3:2695/1‑65 (MQ=255)

                                                                                                                                               atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                                <  3:939019/73‑1 (MQ=255)

                                                                                                                                                tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                               >  4:451977/1‑105 (MQ=255)

                                                                                                                                                tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                               <  3:451977/105‑1 (MQ=255)

                                                                                                                                                  ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt  >  1:238244/1‑148 (MQ=255)

GGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGT  >  NC_000913/4116787‑4117080

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: `‑`