BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,351	2 bp→CA	coding (44‑45/318 nt)	metJ ←	DNA‑binding transcriptional repressor MetJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,351	0	G	C	100.0%	2234.9 / NA	595	H15M (CAC→CTG)	metJ	DNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (270/325); total (270/325)
*	NC_000913	4,128,352	0	T	A	99.8%	2234.3 / NA	595	H15M (CAC→CTG)	metJ	DNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (270/324); minor base C (0/1); total (270/325)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 595 total aligned reads displayed.

CAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCG  >  NC_000913/4128207‑4128495

||

cAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTATGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTc                                                                                                                                             >  2:2719151/1‑150 (MQ=255)

cAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTc                                                                                                                                             >  1:1530773/1‑150 (MQ=255)

 aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCa                                                                                                                                            <  1:2378780/150‑1 (MQ=255)

 aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCa                                                                                                                                            <  2:1510214/150‑1 (MQ=255)

 aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCa                                                                                                                                            >  1:1376451/1‑150 (MQ=255)

   aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGc                                                                                                                                          >  2:2035754/1‑150 (MQ=255)

   aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGc                                                                                                                                          >  1:2654413/1‑150 (MQ=255)

     aCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGt                                                                                                                                        >  2:250338/1‑150 (MQ=255)

     aCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGt                                                                                                                                        >  2:406018/1‑150 (MQ=255)

        cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATg                                                                                                                                     <  2:1187321/150‑1 (MQ=255)

        cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATg                                                                                                                                     >  1:546314/1‑150 (MQ=255)

        cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATg                                                                                                                                     <  2:844246/150‑1 (MQ=255)

        cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATg                                                                                                                                     >  1:2205471/1‑150 (MQ=255)

          tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATggg                                                                                                                                   <  1:1344981/150‑1 (MQ=255)

            gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCt                                                                                                                                 <  1:1136340/150‑1 (MQ=255)

              acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGa                                                                                                                               >  1:1804728/1‑150 (MQ=255)

               cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGat                                                                                                                              <  2:1314619/150‑1 (MQ=255)

               cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGat                                                                                                                              <  1:906692/150‑1 (MQ=255)

                 gcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGatat                                                                                                                            <  1:1930409/150‑1 (MQ=255)

                     cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTc                                                                                                                        >  1:2089629/1‑150 (MQ=255)

                       cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTcgc                                                                                                                      <  2:1042056/150‑1 (MQ=255)

                       cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTcgc                                                                                                                      <  1:2124959/150‑1 (MQ=255)

                            ttGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTc                                                                                                                 >  1:336436/1‑150 (MQ=255)

                             tGGTAGCGTGACGCAGGTTGTTCACCTGTCGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCCGCTCTGCGTATGGGCTGATATATTCGCCGCTcc                                                                                                                <  2:648069/150‑1 (MQ=255)

                              ggTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                                                               >  2:1951605/1‑150 (MQ=255)

                               gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                                                                                                              <  2:413317/150‑1 (MQ=255)

                                   cGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                                                                          <  1:930371/150‑1 (MQ=255)

                                    gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                                                                         <  1:254257/150‑1 (MQ=255)

                                    gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                                                                         >  1:2055110/1‑150 (MQ=255)

                                     tGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCACCGCTCCATTCAGcc                                                                                                        <  1:628281/150‑1 (MQ=255)

                                      gACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa                                                                                                       <  1:1749395/150‑1 (MQ=255)

                                      gACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa                                                                                                       >  1:873254/1‑150 (MQ=255)

                                       aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                                                                      <  2:2342542/150‑1 (MQ=255)

                                       aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                                                                      <  2:1632484/150‑1 (MQ=255)

                                       aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                                                                      <  2:2179241/150‑1 (MQ=255)

                                         gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                                                                    >  1:1236588/1‑150 (MQ=255)

                                         gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                                                                    >  2:1743998/1‑150 (MQ=255)

                                         gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                                                                    >  1:1574717/1‑150 (MQ=255)

                                           aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                                                                  >  1:165504/1‑150 (MQ=255)

                                             gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                                                                >  2:1269622/1‑150 (MQ=255)

                                               tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                                                                              <  1:1307742/150‑1 (MQ=255)

                                                 ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa                                                                                            >  2:810559/1‑150 (MQ=255)

                                                  tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                                                           <  1:1593740/150‑1 (MQ=255)

                                                    aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                                                                         <  1:1455251/150‑1 (MQ=255)

                                                    aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                                                                         <  1:1706106/150‑1 (MQ=255)

                                                      cTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                                                       >  1:1081859/1‑150 (MQ=255)

                                                      cTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                                                       >  1:2474317/1‑150 (MQ=255)

                                                        gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                                                                     <  1:1834671/150‑1 (MQ=255)

                                                           gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                                                                  >  2:1625858/1‑150 (MQ=255)

                                                            acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                 <  2:456316/150‑1 (MQ=255)

                                                             cgcgTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                                                                <  1:2297305/150‑1 (MQ=255)

                                                              gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                                                               >  2:1804941/1‑150 (MQ=255)

                                                                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                             <  2:506188/150‑1 (MQ=255)

                                                                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                             <  1:1015854/150‑1 (MQ=255)

                                                                            gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacca                                                                 >  2:2073847/1‑150 (MQ=255)

                                                                              gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                                                               >  2:2145975/1‑150 (MQ=255)

                                                                               aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                                                              >  1:349794/1‑150 (MQ=255)

                                                                               aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                                                              >  2:337236/1‑150 (MQ=255)

                                                                               aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                                                              <  1:1173111/150‑1 (MQ=255)

                                                                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                                                             >  2:693495/1‑150 (MQ=255)

                                                                                 gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                                                            >  2:547160/1‑150 (MQ=255)

                                                                                 gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                                                            >  2:2448144/1‑150 (MQ=255)

                                                                                 gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                                                            <  2:2566410/150‑1 (MQ=255)

                                                                                   tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcaca                                                          >  2:1377134/1‑150 (MQ=255)

                                                                                    ttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                                                         <  2:2628550/150‑1 (MQ=255)

                                                                                      ttAACACCTTAAGTGGAATTGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTTCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTGCATTCAGCCATGAGATACTTAATCGTCTTCGTCAATAAATTGTGACCAGACCACACAg                                                       <  2:166589/150‑1 (MQ=255)

                                                                                      ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg                                                       >  1:2501716/1‑150 (MQ=255)

                                                                                      ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg                                                       <  2:1086016/150‑1 (MQ=255)

                                                                                        aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                                                     >  2:2003267/1‑150 (MQ=255)

                                                                                         acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                    <  1:823867/150‑1 (MQ=255)

                                                                                           acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt                                                  <  2:1501788/150‑1 (MQ=255)

                                                                                            ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                                                 <  1:1990031/150‑1 (MQ=255)

                                                                                             cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt                                                <  2:154491/150‑1 (MQ=255)

                                                                                               tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                              <  1:462102/150‑1 (MQ=255)

                                                                                               tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                              <  1:1406873/150‑1 (MQ=255)

                                                                                               tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                              <  1:1273273/150‑1 (MQ=255)

                                                                                                   agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAc                                          <  2:2272926/150‑1 (MQ=255)

                                                                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt                                         >  2:24019/1‑150 (MQ=255)

                                                                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt                                         >  2:2668509/1‑150 (MQ=255)

                                                                                                        tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                     <  1:500248/150‑1 (MQ=255)

                                                                                                          gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                   >  1:772129/1‑150 (MQ=255)

                                                                                                           aaaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTa                                  >  1:678766/1‑150 (MQ=255)

                                                                                                            aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTaa                                 <  2:2110637/150‑1 (MQ=255)

                                                                                                            aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTaa                                 >  2:1151090/1‑150 (MQ=255)

                                                                                                              ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAc                               <  2:882672/150‑1 (MQ=255)

                                                                                                                gTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAccc                             <  2:1890220/150‑1 (MQ=255)

                                                                                                                gTAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAccc                             >  2:101223/1‑150 (MQ=255)

                                                                                                                 tAATCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCg                            >  2:1238251/1‑150 (MQ=255)

                                                                                                                   aTCTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCa                          <  2:1470132/150‑1 (MQ=255)

                                                                                                                     cTTTTTGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATa                        <  1:2608850/150‑1 (MQ=255)

                                                                                                                       ttttGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGt                      >  1:2248185/1‑150 (MQ=255)

                                                                                                                         ttGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGttt                    <  2:834568/150‑1 (MQ=255)

                                                                                                                          tGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTa                   <  2:320652/150‑1 (MQ=255)

                                                                                                                          tGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTa                   <  1:2692141/150‑1 (MQ=255)

                                                                                                                          tGACTTGTTCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTa                   <  1:2423630/150‑1 (MQ=255)

                                                                                                                                 ttCACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACa            >  2:232205/1‑150 (MQ=255)

                                                                                                                                   cACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg          >  2:353114/1‑150 (MQ=255)

                                                                                                                                   cACTCTTCTTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg          >  1:2018268/1‑150 (MQ=255)

                                                                                                                                      tcttctTGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGt       >  1:285593/1‑150 (MQ=255)

                                                                                                                                         tcttGCCCAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAc    >  1:4073/1‑150 (MQ=255)

||

CAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCG  >  NC_000913/4128207‑4128495

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