BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,516	G→T	72.6%	T236K (ACG→AAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,516	0	G	T	72.6%	217.5 / 131.9	128	T236K (ACG→AAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (19/16); new base T (52/41); total (71/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 130 total aligned reads displayed.

TTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAG  >  NC_000913/4116369‑4116662

ttttCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                                                                                                                   >  3:427772/1‑149 (MQ=255)

     tGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTgccgc                                                                                                                                              >  3:359508/1‑149 (MQ=255)

          tCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTg                                                                                                                                         >  4:330913/1‑149 (MQ=255)

          tCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTg                                                                                                                                         >  1:288549/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgcc                                                                                                                                       >  1:499367/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgcc                                                                                                                                       >  2:520655/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgcc                                                                                                                                       >  2:244069/1‑149 (MQ=255)

                cTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCa                                                                                                                                   >  1:400985/1‑149 (MQ=255)

                  ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                                                                                                 <  3:302443/149‑1 (MQ=255)

                      cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                             <  4:591046/149‑1 (MQ=255)

                      cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                             <  2:7298/149‑1 (MQ=255)

                      cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                             <  3:330913/149‑1 (MQ=255)

                        ccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAGCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                           <  4:655958/149‑1 (MQ=255)

                          aGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCt                                                                                                                         <  2:400985/149‑1 (MQ=255)

                             gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                           >  2:128452/1‑144 (MQ=255)

                             gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                           <  1:128452/144‑1 (MQ=255)

                             gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                      >  3:781418/1‑149 (MQ=255)

                             gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                      <  1:520655/149‑1 (MQ=255)

                                catcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGt                                                                                                                   >  4:228393/1‑149 (MQ=255)

                                    aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATAc                                                                                                               >  2:63940/1‑149 (MQ=255)

                                       aTAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAAt                                                                                                                                 <  3:515931/128‑1 (MQ=255)

                                       aTAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAAt                                                                                                                                 >  4:515931/1‑128 (MQ=255)

                                          aaGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCt                                                                                                                         >  3:600846/1‑133 (MQ=255)

                                          aaGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGga                                                                                                         >  4:396454/1‑149 (MQ=255)

                                           agcagcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                        <  1:554811/149‑1 (MQ=255)

                                           agcagcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                        <  3:396454/149‑1 (MQ=255)

                                               gcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                                                                                    >  3:728498/1‑149 (MQ=255)

                                                  aGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                                                                                                 <  2:39291/149‑1 (MQ=255)

                                                  aGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                                                                                                 <  4:68263/149‑1 (MQ=255)

                                                       cATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                                                                                            >  1:516497/1‑149 (MQ=255)

                                                       cATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                                                                                            >  1:460506/1‑149 (MQ=255)

                                                       cATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                                                                                            >  2:99283/1‑149 (MQ=255)

                                                             tgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                      <  2:218764/117‑1 (MQ=255)

                                                              gtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                     <  3:770576/149‑1 (MQ=255)

                                                              gtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                     <  4:750309/149‑1 (MQ=255)

                                                                  ttcGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt                                                                                 >  3:706931/1‑149 (MQ=255)

                                                                     gCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGtt                                                                                                                              >  1:783787/1‑101 (MQ=255)

                                                                           ccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtctc                                                                              >  4:342989/1‑143 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                       >  2:258515/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                       >  1:112604/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                       >  3:707347/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                       >  3:454628/1‑149 (MQ=255)

                                                                                tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat                                                                   >  3:37606/1‑149 (MQ=255)

                                                                                 ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATc                                                                  >  1:422823/1‑149 (MQ=255)

                                                                                   cACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCa                                                                <  1:258515/149‑1 (MQ=255)

                                                                                    aCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                               <  3:489354/149‑1 (MQ=255)

                                                                                    aCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                               <  2:511225/149‑1 (MQ=255)

                                                                                      gCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGAAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCa                                                             >  4:146719/1‑149 (MQ=255)

                                                                                        acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACt                                                           >  1:219250/1‑149 (MQ=255)

                                                                                          aaCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTc                                                         >  3:325438/1‑149 (MQ=255)

                                                                                               aCCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCa                                                    >  1:347280/1‑149 (MQ=255)

                                                                                               aCCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCa                                                    >  1:408986/1‑149 (MQ=255)

                                                                                                ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAt                                                   >  4:552029/1‑149 (MQ=255)

                                                                                                 cAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtt                                                  >  1:31051/1‑149 (MQ=255)

                                                                                                   aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                >  3:575085/1‑149 (MQ=255)

                                                                                                   aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                >  1:271941/1‑149 (MQ=255)

                                                                                                      aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                             <  3:339691/149‑1 (MQ=255)

                                                                                                      aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                             <  1:136040/149‑1 (MQ=255)

                                                                                                       gcgcGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtcg                                            >  4:717539/1‑149 (MQ=255)

                                                                                                        cgcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtcgt                                           <  4:339473/149‑1 (MQ=255)

                                                                                                            gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCa                                       >  1:365452/1‑149 (MQ=255)

                                                                                                              ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGa                                                                                                                       >  3:28523/1‑67 (MQ=255)

                                                                                                              ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGa                                                                                                                       <  4:28523/67‑1 (MQ=255)

                                                                                                                 ctgGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCa                                  <  3:622634/149‑1 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                >  1:122983/1‑149 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                                                                                  <  1:643462/83‑1 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                                                                                  >  2:643462/1‑83 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                <  4:427772/149‑1 (MQ=255)

                                                                                                                    gTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAggg                               >  3:449936/1‑149 (MQ=255)

                                                                                                                      cACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTa                             <  2:422823/149‑1 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                            >  3:686099/1‑149 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                            <  2:516497/149‑1 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                            <  2:460506/149‑1 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                            >  4:319958/1‑149 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                            <  2:288549/149‑1 (MQ=255)

                                                                                                                         cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                                                                                          >  1:60163/1‑69 (MQ=255)

                                                                                                                         cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTcc                                         >  1:528953/1‑134 (MQ=255)

                                                                                                                         cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTcc                                         <  2:528953/134‑1 (MQ=255)

                                                                                                                           gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                     <  3:194025/88‑1 (MQ=255)

                                                                                                                           gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGAt                        <  2:806325/149‑1 (MQ=255)

                                                                                                                           gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGAt                        <  1:244069/149‑1 (MQ=255)

                                                                                                                             gATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGt                      >  2:149610/1‑149 (MQ=255)

                                                                                                                              aTCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                     >  3:257226/1‑149 (MQ=255)

                                                                                                                              aTCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                     <  1:601537/149‑1 (MQ=255)

                                                                                                                              aTCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                     >  3:54963/1‑149 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                  >  4:552879/1‑149 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                  <  4:706931/149‑1 (MQ=255)

                                                                                                                                      gaTTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAacaaca             <  4:325438/149‑1 (MQ=255)

                                                                                                                                       aTTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACAt            >  4:214796/1‑149 (MQ=255)

                                                                                                                                       aTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACAt            >  4:376699/1‑149 (MQ=255)

                                                                                                                                        ttGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGATGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg           <  2:165358/149‑1 (MQ=255)

                                                                                                                                        ttGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg           <  4:70907/149‑1 (MQ=255)

                                                                                                                                          tgAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt         <  2:271941/148‑1 (MQ=255)

                                                                                                                                          ggAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt         <  4:575085/149‑1 (MQ=255)

                                                                                                                                          ggAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt         >  3:301623/1‑149 (MQ=255)

                                                                                                                                           gAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa        >  2:431351/1‑149 (MQ=255)

                                                                                                                                           gAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa        <  2:112604/149‑1 (MQ=255)

                                                                                                                                           gAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTACTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa        <  4:454628/149‑1 (MQ=255)

                                                                                                                                                cgctTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACgaga   <  3:425787/149‑1 (MQ=255)

                                                                                                                                                 gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGagag  <  2:408986/149‑1 (MQ=255)

TTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAG  >  NC_000913/4116369‑4116662

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