BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,005	T→A	5.5%	D73V (GAT→GTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,005	0	T	A	5.5%	367.8 / 10.0	109	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (49/54); new base A (3/3); total (52/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.68e-01

Only 100 of 110 total aligned reads displayed.

TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGT  >  NC_000913/4116857‑4117152

ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                     >  3:246366/1‑149 (MQ=255)

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

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

 tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATCTGATc                                                                                                                                                    <  3:39620/149‑1 (MQ=255)

   tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCgg                                                                                                                                                  <  4:437541/149‑1 (MQ=255)

     aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                                                                                                >  1:583279/1‑149 (MQ=255)

     aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                                                                                                >  2:115337/1‑149 (MQ=255)

       ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                                                                                              <  1:353319/149‑1 (MQ=255)

          tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                                                                                           >  1:95587/1‑149 (MQ=255)

          tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                                                                                           >  4:289655/1‑149 (MQ=255)

            aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGatat                                                                                                                                         >  1:531242/1‑149 (MQ=255)

              gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                                                                                       >  1:648538/1‑149 (MQ=255)

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

                 ttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCt                                                                                                                                    >  4:613487/1‑149 (MQ=255)

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

                         tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAg                                                                                                                            >  4:488621/1‑149 (MQ=255)

                           gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGcc                                                                                                                          >  1:135318/1‑148 (MQ=255)

                           gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                          <  2:28127/149‑1 (MQ=255)

                           gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                          >  3:194569/1‑149 (MQ=255)

                                  tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                   <  1:569867/149‑1 (MQ=255)

                                  tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                   <  3:613487/149‑1 (MQ=255)

                                        gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                                                                                             >  1:151757/1‑149 (MQ=255)

                                               ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCtt                                                                                                      <  1:251975/149‑1 (MQ=255)

                                               ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCtt                                                                                                      >  2:443811/1‑149 (MQ=255)

                                                gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCttt                                                                                                     <  3:255049/149‑1 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                   <  1:463973/149‑1 (MQ=255)

                                                   cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                                  <  1:328097/149‑1 (MQ=255)

                                                       cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                          >  2:531831/1‑121 (MQ=255)

                                                       cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                          <  1:531831/121‑1 (MQ=255)

                                                        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGcc                                                                                             >  2:137386/1‑149 (MQ=255)

                                                              tGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGAtt                                                                                       <  2:358325/149‑1 (MQ=255)

                                                                gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCttt                                                                                                     <  3:497717/133‑1 (MQ=255)

                                                                gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCttt                                                                                                     >  4:497717/1‑133 (MQ=255)

                                                                 cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTcc                                                                                    <  4:416057/149‑1 (MQ=255)

                                                                 cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTcc                                                                                    <  2:328433/149‑1 (MQ=255)

                                                                    tataGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCAtt                                                                                 <  3:509465/149‑1 (MQ=255)

                                                                         atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACTCTGGTTCGCAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                            <  4:431359/149‑1 (MQ=255)

                                                                         atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                            >  3:280372/1‑149 (MQ=255)

                                                                           agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                          >  1:152895/1‑149 (MQ=255)

                                                                           agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                          <  3:488621/149‑1 (MQ=255)

                                                                           agGCTTGCAGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGAGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                          <  2:698808/149‑1 (MQ=255)

                                                                             gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGTACTGATATCGGCTCTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg                                                                         <  4:558002/148‑1 (MQ=255)

                                                                             gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                        >  1:552384/1‑149 (MQ=255)

                                                                               ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                      >  3:407125/1‑149 (MQ=255)

                                                                                tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt                                                                     <  4:261498/149‑1 (MQ=255)

                                                                                     gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCa                                                                >  3:237844/1‑149 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                             >  1:86234/1‑135 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                             <  2:86234/135‑1 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                                               >  1:686569/1‑149 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                                               <  4:213103/149‑1 (MQ=255)

                                                                                          tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                           >  2:199725/1‑149 (MQ=255)

                                                                                          tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                           >  3:20859/1‑149 (MQ=255)

                                                                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAGCTgg                                                          >  3:660101/1‑149 (MQ=255)

                                                                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                                          >  1:303010/1‑149 (MQ=255)

                                                                                             ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                        <  4:111538/149‑1 (MQ=255)

                                                                                                ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     >  4:21990/1‑149 (MQ=255)

                                                                                                ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     <  3:45390/149‑1 (MQ=255)

                                                                                                ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     <  1:487864/149‑1 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    >  1:410223/1‑149 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    >  4:153582/1‑149 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    <  4:503662/149‑1 (MQ=255)

                                                                                                  aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGTTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    >  1:175302/1‑148 (MQ=255)

                                                                                                  aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                   >  1:245425/1‑149 (MQ=255)

                                                                                                  aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    <  2:175302/148‑1 (MQ=255)

                                                                                                    aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                                 >  2:387383/1‑149 (MQ=255)

                                                                                                        tAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAtt                                             >  1:188611/1‑149 (MQ=255)

                                                                                                           tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGAt                                                                                        <  3:70325/103‑1 (MQ=255)

                                                                                                           tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGAt                                                                                        >  4:70325/1‑103 (MQ=255)

                                                                                                            ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCt                                         <  4:569717/149‑1 (MQ=255)

                                                                                                             gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGAGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                        >  3:409690/1‑149 (MQ=255)

                                                                                                               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa                                      >  4:295290/1‑149 (MQ=255)

                                                                                                                tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                     <  2:648538/149‑1 (MQ=255)

                                                                                                                 cACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGCGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAt                                    <  3:322649/149‑1 (MQ=255)

                                                                                                                  aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                                   <  3:89610/149‑1 (MQ=255)

                                                                                                                    gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcg                                 <  2:449176/149‑1 (MQ=255)

                                                                                                                         ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                            >  3:308049/1‑149 (MQ=255)

                                                                                                                          tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                           <  1:215367/149‑1 (MQ=255)

                                                                                                                          tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCATGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                           <  3:465547/149‑1 (MQ=255)

                                                                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                        >  1:181329/1‑149 (MQ=255)

                                                                                                                             tAATACCGATAACTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                        >  3:441354/1‑149 (MQ=255)

                                                                                                                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    <  1:669864/117‑1 (MQ=255)

                                                                                                                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    >  2:669864/1‑117 (MQ=255)

                                                                                                                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGc                    <  1:199725/149‑1 (MQ=255)

                                                                                                                                   cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa                  >  3:625241/1‑149 (MQ=255)

                                                                                                                                    gATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGCTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                 <  4:409690/149‑1 (MQ=255)

                                                                                                                                    gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                 >  4:259/1‑149 (MQ=255)

                                                                                                                                     aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                <  2:361621/149‑1 (MQ=255)

                                                                                                                                     aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATATGCTCAAATTCGCGCTGCGACACGCTAAt                <  2:472238/149‑1 (MQ=255)

                                                                                                                                        gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGat             <  1:115337/149‑1 (MQ=255)

                                                                                                                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat           <  2:196068/149‑1 (MQ=255)

                                                                                                                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat           <  3:438360/149‑1 (MQ=255)

                                                                                                                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCCGCGACACGCTAATGatat           >  2:537068/1‑149 (MQ=255)

                                                                                                                                          tGCAACTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat           <  4:441354/149‑1 (MQ=255)

                                                                                                                                            cAATTTGATCGGAACTGATATCGGCTTTCGCCCGCACTTCTCCCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg         <  4:409047/149‑1 (MQ=255)

                                                                                                                                             aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTgg        <  4:299795/149‑1 (MQ=255)

                                                                                                                                              aTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGc       <  4:248176/149‑1 (MQ=255)

                                                                                                                                               tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa      >  4:147011/1‑149 (MQ=255)

                                                                                                                                                 tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc    <  2:181329/149‑1 (MQ=255)

                                                                                                                                                   aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt  >  2:403784/1‑149 (MQ=255)

                                                                                                                                                   aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt  <  2:526913/149‑1 (MQ=255)

TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGT  >  NC_000913/4116857‑4117152

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