BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,697	G→T	33.8%	Q176K (CAG→AAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,697	0	G	T	33.8%	434.7 / 296.5	278	Q176K (CAG→AAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (91/93); new base T (57/37); total (148/130)
Fisher's exact test for biased strand distribution p-value = 9.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-02

Only 100 of 280 total aligned reads displayed.

TATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTG  >  NC_000913/4116549‑4116845

 ataCTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                                                                                                                                                     <  4:1916273/149‑1 (MQ=255)

   aCTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGt                                                                                                                                                   <  2:265892/149‑1 (MQ=255)

       cGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCAtt                                                                                                                                               <  1:357460/149‑1 (MQ=255)

       cGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtt                                                                                                                                               <  4:1846214/149‑1 (MQ=255)

         gaagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                                                                                             >  1:1430030/1‑149 (MQ=255)

         gaagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                                                                                             <  2:1566598/149‑1 (MQ=255)

            gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                                                                          <  4:928302/149‑1 (MQ=255)

            gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGCAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCa                                                                                                                                          >  2:1018550/1‑149 (MQ=255)

                 acgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATaa                                                                                                                                     <  1:1734519/149‑1 (MQ=255)

                   gCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc                                                                                                                                   <  4:58138/149‑1 (MQ=255)

                      cTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCAc                                                                                                                                <  4:2418489/149‑1 (MQ=255)

                         cTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCCGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTAAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                                                                                             <  3:1034496/149‑1 (MQ=255)

                           ggCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATc                                                                                                                           <  3:24336/149‑1 (MQ=255)

                            gcagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa                                                                                                                          >  4:2177929/1‑149 (MQ=255)

                             cagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCaa                                                                                                                         <  1:244387/149‑1 (MQ=255)

                              agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAg                                                                                                                                         <  2:2341890/132‑1 (MQ=255)

                              agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAc                                                                                                                        >  2:3945/1‑149 (MQ=255)

                                caTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCg                                                                                                                      >  1:1979892/1‑149 (MQ=255)

                                 aTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGt                                                                                                                     >  4:424189/1‑149 (MQ=255)

                                 aTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGt                                                                                                                     <  2:2031043/149‑1 (MQ=255)

                                        gcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAc                                                                                                              >  2:280664/1‑149 (MQ=255)

                                         cgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAATCCACGTATCAACCGTaacaaaca                                                                                                             <  2:332013/149‑1 (MQ=255)

                                         cgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCTGATAAGCCACGTATCAACCGTACCAAACa                                                                                                             <  4:258844/149‑1 (MQ=255)

                                           gAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGc                                                                                                           >  1:645543/1‑149 (MQ=255)

                                           gAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGc                                                                                                           >  3:705360/1‑149 (MQ=255)

                                           gAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGc                                                                                                           >  1:2501810/1‑149 (MQ=255)

                                             atatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                                                                                         <  3:524909/149‑1 (MQ=255)

                                              tatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAt                                                                                                        <  4:293611/149‑1 (MQ=255)

                                               atCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAtt                                                                                                       <  4:691059/149‑1 (MQ=255)

                                                tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTc                                                                                                      <  2:2461522/149‑1 (MQ=255)

                                                tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTc                                                                                                      <  2:1487741/149‑1 (MQ=255)

                                                  cttaaCTTCGAGAATGTTGCCGTCCGAGTACAGGGTAAGGATCTTGAACAACATGGTAAGGGAGGCGTTGGTGTGATCGGTCACGTGGACACGGCACTTAGTCATTTTCCAGATAAGCCAGGTATCAACCGTACCAAACAACAATTcac                                                                                                    >  4:1364072/5‑149 (MQ=255)

                                                  cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                                                                                                    >  3:1403907/1‑149 (MQ=255)

                                                   aGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcacc                                                                                                   >  3:1731124/1‑149 (MQ=255)

                                                   aGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAg                                                                                                                                         >  3:2037237/1‑111 (MQ=255)

                                                   aGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAg                                                                                                                                         <  4:2037237/111‑1 (MQ=255)

                                                    gCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcacca                                                                                                  >  4:142634/1‑149 (MQ=255)

                                                    gCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                                                                                                                                                     >  4:2295574/1‑98 (MQ=255)

                                                      aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTGGAACAACATGGTACGAGAGGCGTTGGCGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacg                                                                                                >  1:292822/1‑149 (MQ=255)

                                                        ttCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacgac                                                                                              >  4:1978958/1‑149 (MQ=255)

                                                          ccAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                                                                                            <  4:1731124/149‑1 (MQ=255)

                                                           cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                           >  2:1201014/1‑149 (MQ=255)

                                                           cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                           >  3:1846555/1‑149 (MQ=255)

                                                            aGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGc                                                                                          >  1:1915385/1‑149 (MQ=255)

                                                            aGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTc                                                                                                      >  1:2515199/1‑137 (MQ=255)

                                                               aTTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACg                                                                                       >  1:1807251/1‑149 (MQ=255)

                                                                ttttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGc                                                                                      >  3:774505/1‑149 (MQ=255)

                                                                    gtcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgc                                                                                  >  3:2311355/1‑149 (MQ=255)

                                                                      cgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCga                                                                                <  2:1960580/149‑1 (MQ=255)

                                                                         cccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcg                                                                                   >  2:2219197/1‑143 (MQ=255)

                                                                         cccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcg                                                                                   <  1:2219197/143‑1 (MQ=255)

                                                                         cccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGagag                                                                             <  4:1266015/149‑1 (MQ=255)

                                                                          ccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGc                                                                            <  2:1807251/149‑1 (MQ=255)

                                                                           cAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCAc                                                                                        >  2:620783/1‑136 (MQ=255)

                                                                           cAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCAc                                                                                        <  1:620783/136‑1 (MQ=255)

                                                                               ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTcc                                                                       <  1:2330014/149‑1 (MQ=255)

                                                                                cAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCca                                                                      <  4:798876/149‑1 (MQ=255)

                                                                                    gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATg                                                                  <  2:250995/149‑1 (MQ=255)

                                                                                       tGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTc                                                               >  1:685068/1‑149 (MQ=255)

                                                                                         gATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTTATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                             >  3:2015055/1‑149 (MQ=255)

                                                                                          aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCcaca                                                                    >  3:1950303/1‑141 (MQ=255)

                                                                                           tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCa                                                                                                                                                 >  1:1472935/1‑63 (MQ=255)

                                                                                           tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAgg                                                           >  3:473857/1‑149 (MQ=255)

                                                                                            gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                          <  3:1200679/149‑1 (MQ=255)

                                                                                            gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                          >  1:208705/1‑149 (MQ=255)

                                                                                              tGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                        <  1:1691463/149‑1 (MQ=255)

                                                                                               gAACAATATGGTACGGGAGGCGTTGGTATCATCGGTCACATGGACACGGCCCTGAGTCATTTACTAGATCAGCCCCGTATGAACAGTACCAAACAGCAAGGCGCAACGACGTGCAGCCTCGCGAGAACCATCCACATGGACGAGGATcc                                                       >  4:1142571/1‑149 (MQ=255)

                                                                                                 acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                     <  4:1647582/149‑1 (MQ=255)

                                                                                                 acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                     <  4:1584395/149‑1 (MQ=255)

                                                                                                   aacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcactt                                                   >  4:451419/1‑149 (MQ=255)

                                                                                                   aacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcactt                                                   >  2:1051113/1‑149 (MQ=255)

                                                                                                          tACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTg                                            <  3:142634/149‑1 (MQ=255)

                                                                                                           aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                           >  3:19110/1‑149 (MQ=255)

                                                                                                            cGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGt                                          <  1:1002579/149‑1 (MQ=255)

                                                                                                                  gCGTTGGTGTAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCATGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGa                                    >  2:2314031/1‑149 (MQ=255)

                                                                                                                     ttGGTGTAAGCGGTCACATGGACACGGCCCAGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAGAGCAATTCCCCACCACGTGCACGGTCGCGAGAGCCTTCCACATGGTCGAAGATCCACTGCACTTTGGTCCCAGaaaa                                 <  2:2168688/149‑1 (MQ=255)

                                                                                                                        gtgtAATCGGTTACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCAACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGATAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTa                              >  2:911191/1‑149 (MQ=255)

                                                                                                                         tgtAATCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcact                                               >  4:1116587/1‑131 (MQ=255)

                                                                                                                         tgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTAc                             <  4:705360/149‑1 (MQ=255)

                                                                                                                          gtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                            <  3:424189/149‑1 (MQ=255)

                                                                                                                          gtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACATGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                            <  2:835674/149‑1 (MQ=255)

                                                                                                                           tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                           <  3:284068/149‑1 (MQ=255)

                                                                                                                             aTCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                               >  1:1037907/1‑143 (MQ=255)

                                                                                                                             aTCGGTCACATGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAAGGCAATTCACCACGACGGGCACGCTCGCGAGAGCCTCCAACATGGTCGAGGATCCACTTCACTTTGTTGCCAGAAAAGTACGGGt                         >  1:786234/1‑149 (MQ=255)

                                                                                                                             aTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                         <  4:1947524/117‑1 (MQ=255)

                                                                                                                                ggTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                        <  3:1248345/115‑1 (MQ=255)

                                                                                                                                    acaTGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                  >  3:663985/1‑149 (MQ=255)

                                                                                                                                    acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                  >  3:2280980/1‑149 (MQ=255)

                                                                                                                                     caTGGACACGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                            >  3:1790138/1‑106 (MQ=255)

                                                                                                                                     caTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAcc                 <  2:685068/149‑1 (MQ=255)

                                                                                                                                      aTGGACACGGCCCTTAGTCATTTTCCAGATAATCCACGTATCAACCGTACCACACAGCAATTCACCACGACGTGCACGCCCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                >  2:1422040/1‑149 (MQ=255)

                                                                                                                                      aTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                <  1:1273161/149‑1 (MQ=255)

                                                                                                                                       tGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                            <  2:1779986/104‑1 (MQ=255)

                                                                                                                                          acacGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTcc                                                                       >  1:707994/1‑90 (MQ=255)

                                                                                                                                          acacGGCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc            >  1:2344568/1‑149 (MQ=255)

                                                                                                                                          acacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                             <  3:1647124/100‑1 (MQ=255)

                                                                                                                                            acGGCCCTGAGTCATTTTCTAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcctg          >  4:2132267/1‑147 (MQ=255)

                                                                                                                                              ggCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                        <  2:699721/101‑1 (MQ=255)

                                                                                                                                              ggCCCTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTa        >  3:889178/1‑149 (MQ=255)

                                                                                                                                                 ccTTAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtg     >  3:304310/1‑149 (MQ=255)

TATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTG  >  NC_000913/4116549‑4116845

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