BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,937	T→C	84.4%	K96E (AAG→GAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,937	0	T	C	84.4%	256.7 / 33.9	102	K96E (AAG→GAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (12/4); new base C (45/41); total (57/45)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 103 total aligned reads displayed.

CACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCT  >  NC_000913/4116792‑4117085

cacttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGcc                                                                                                                                                   >  1:232198/1‑149 (MQ=255)

 acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCg                                                                                                                                                  >  2:186714/1‑149 (MQ=255)

   ttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGt                                                                                                                                                >  1:138507/1‑149 (MQ=255)

       cttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCt                                                                                                                                            <  4:120797/149‑1 (MQ=255)

       cttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCt                                                                                                                                            <  4:502585/149‑1 (MQ=255)

        tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCtt                                                                                                                                           >  4:782209/1‑149 (MQ=255)

        tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCtt                                                                                                                                           >  4:582517/1‑149 (MQ=255)

           ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCttttt                                                                                                                                        <  3:606208/149‑1 (MQ=255)

            gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGTCTCGCCGGTTTCTTTTTc                                                                                                                                       <  4:465072/149‑1 (MQ=255)

            gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTc                                                                                                                                       <  4:157391/149‑1 (MQ=255)

            gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTc                                                                                                                                       <  4:281751/149‑1 (MQ=255)

             tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTcc                                                                                                                                      <  2:191290/149‑1 (MQ=255)

                      aGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATa                                                                                                                             <  4:239960/149‑1 (MQ=255)

                        tACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGt                                                                                                                           >  4:642637/1‑149 (MQ=255)

                          cGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTgg                                                                                                                         <  4:63425/149‑1 (MQ=255)

                             gTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGttt                                                                                                                      >  2:659599/1‑149 (MQ=255)

                             gTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGttt                                                                                                                      <  2:91899/149‑1 (MQ=255)

                               cAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCa                                                                                                                    >  4:64466/1‑149 (MQ=255)

                                aaacACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGGCACGTCTTAAAGCCTAGAAGATTTCTGCGGTACGACGGCACTGCCAGACACTGGCGTTATACATAGGCGCGCCGGTTTCTTTTTCACAGACAATAGTTGTTTCAc                                                                                                                   <  2:562706/146‑1 (MQ=255)

                                aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                                   >  3:447909/1‑149 (MQ=255)

                                aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCACGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                                   >  3:753174/1‑149 (MQ=255)

                                 aTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACg                                                                                                                  >  1:737070/1‑149 (MQ=255)

                                 aTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCa                                                                                                                    >  4:650299/1‑147 (MQ=255)

                                 aTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCa                                                                                                                    <  3:650299/147‑1 (MQ=255)

                                 aTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACg                                                                                                                  >  2:69992/1‑149 (MQ=255)

                                     ccAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                                                                                                              <  1:607321/149‑1 (MQ=255)

                                      cAGACCGGTATTGCTGCGGATATAATCTTCTTAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                                                                                                             <  3:723593/149‑1 (MQ=255)

                                      cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                                                                                                             >  1:62906/1‑149 (MQ=255)

                                             gTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAt                                                                                                      >  1:713443/1‑149 (MQ=255)

                                             gTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAt                                                                                                      >  3:278213/1‑149 (MQ=255)

                                               aTTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                                                                                                    <  2:375184/149‑1 (MQ=255)

                                                    tgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa                                                                                               >  4:85059/1‑149 (MQ=255)

                                                       ggATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCt                                                                                            <  1:608431/149‑1 (MQ=255)

                                                            taATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt                                                                                       >  2:76479/1‑149 (MQ=255)

                                                             aaTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                                      >  2:829201/1‑149 (MQ=255)

                                                               tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                    >  2:436907/1‑149 (MQ=255)

                                                               tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                    >  4:607674/1‑149 (MQ=255)

                                                               tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                    >  4:670970/1‑149 (MQ=255)

                                                                 ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                  <  3:185547/149‑1 (MQ=255)

                                                                 ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                  <  3:159279/149‑1 (MQ=255)

                                                                 ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                  >  1:647245/1‑149 (MQ=255)

                                                                  tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                 <  3:256081/149‑1 (MQ=255)

                                                                  tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                 <  1:268257/149‑1 (MQ=255)

                                                                    tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCgg                                                                               >  4:599105/1‑149 (MQ=255)

                                                                    tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCgg                                                                               >  2:240716/1‑149 (MQ=255)

                                                                      aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                             <  4:682157/149‑1 (MQ=255)

                                                                      aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                             <  2:27611/149‑1 (MQ=255)

                                                                       aCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAAc                                                                            >  4:292265/1‑149 (MQ=255)

                                                                          gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGa                                                                         <  3:512293/149‑1 (MQ=255)

                                                                           tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                        >  3:752614/1‑149 (MQ=255)

                                                                               gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCaa                                                                                        >  1:563366/1‑129 (MQ=255)

                                                                               gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCaa                                                                                        <  2:563366/129‑1 (MQ=255)

                                                                                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCa                                                                                                                                    <  3:517866/84‑1 (MQ=255)

                                                                                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCa                                                                                                                                    >  4:517866/1‑84 (MQ=255)

                                                                                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGa                                                                                   <  3:116238/133‑1 (MQ=255)

                                                                                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGa                                                                                   >  4:116238/1‑133 (MQ=255)

                                                                                   tAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                                >  2:312828/1‑149 (MQ=255)

                                                                                     aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                              >  3:306294/1‑149 (MQ=255)

                                                                                     aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                              <  2:138507/149‑1 (MQ=255)

                                                                                        gCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGcc                                                           <  4:54591/149‑1 (MQ=255)

                                                                                         cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                          <  3:529301/149‑1 (MQ=255)

                                                                                          tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAg                                                         >  3:319943/1‑149 (MQ=255)

                                                                                            gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                                                                                                    <  4:661319/104‑1 (MQ=255)

                                                                                            gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                                                                                                    >  3:661319/1‑104 (MQ=255)

                                                                                            gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                       >  2:91701/1‑149 (MQ=255)

                                                                                            gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                       >  4:439005/1‑149 (MQ=255)

                                                                                                aTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTc                                                   <  4:137311/149‑1 (MQ=255)

                                                                                                aTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTc                                                   >  2:551734/1‑149 (MQ=255)

                                                                                                         gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                          <  1:712325/149‑1 (MQ=255)

                                                                                                         gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                          >  1:679059/1‑149 (MQ=255)

                                                                                                          tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                         >  3:477526/1‑149 (MQ=255)

                                                                                                             gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                      >  1:621147/1‑149 (MQ=255)

                                                                                                              acgGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                     >  2:603054/1‑149 (MQ=255)

                                                                                                               cgGCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCt                                    <  1:603054/149‑1 (MQ=255)

                                                                                                                 gCACTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                  <  1:186714/149‑1 (MQ=255)

                                                                                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                <  2:852356/149‑1 (MQ=255)

                                                                                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                >  1:303498/1‑149 (MQ=255)

                                                                                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTCGCTGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                <  3:71308/149‑1 (MQ=255)

                                                                                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                <  1:46124/149‑1 (MQ=255)

                                                                                                                     tGCCAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGt                              <  1:560842/149‑1 (MQ=255)

                                                                                                                      gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                             >  3:110760/1‑149 (MQ=255)

                                                                                                                       ccAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                            >  4:243871/1‑149 (MQ=255)

                                                                                                                       ccAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                            >  2:750465/1‑149 (MQ=255)

                                                                                                                        cAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                 >  2:869167/1‑95 (MQ=255)

                                                                                                                        cAGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                 <  1:869167/95‑1 (MQ=255)

                                                                                                                         aGACAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                          <  2:98906/149‑1 (MQ=255)

                                                                                                                           aCAATGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCa                        <  2:393709/149‑1 (MQ=255)

                                                                                                                             aaTGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                      >  1:409255/1‑149 (MQ=255)

                                                                                                                              aTGGCGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAt                     <  1:632994/149‑1 (MQ=255)

                                                                                                                                 gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc                  >  1:93444/1‑149 (MQ=255)

                                                                                                                                  cGTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                 <  1:71791/149‑1 (MQ=255)

                                                                                                                                   gTTATAGATAGGCTCGCCGGTTTCTTTTTCCCAAACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                >  3:548813/1‑149 (MQ=255)

                                                                                                                                    ttATAGATAGGCTCGCCGGTTTCTTTTTCCCAGACACTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTTGCCCAGATTTCCAt               <  4:510196/149‑1 (MQ=255)

                                                                                                                                        agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg           >  2:838339/1‑149 (MQ=255)

                                                                                                                                        agatagGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGCGGCCCAGATTTCCATTggg           >  1:724602/1‑149 (MQ=255)

                                                                                                                                         gatagGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGCACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt          >  4:647129/1‑149 (MQ=255)

                                                                                                                                          atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc         >  2:581809/1‑149 (MQ=255)

                                                                                                                                            agGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt       >  2:436633/1‑149 (MQ=255)

                                                                                                                                            agGCTCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt       >  3:458317/1‑149 (MQ=255)

                                                                                                                                                 tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt  <  3:495088/149‑1 (MQ=255)

CACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCT  >  NC_000913/4116792‑4117085

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