BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,999,231	T→A	5.5%	T84S (ACT→TCT)	dcyD ←	D‑cysteine desulfhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,231	0	T	A	5.5%	435.8 / 9.7	145	T84S (ACT→TCT)	dcyD	D‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand): ref base T (66/69); new base A (5/3); total (72/73)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 146 total aligned reads displayed.

GGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAA  >  NC_000913/1999082‑1999377

ggTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTG‑CAATAGGATTTTGCAGCAGCGCCCCGGAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGtg                                                                                                                                                   >  1:1659546/1‑149 (MQ=255)

   cAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTg                                                                                                                                                 <  1:1884024/150‑1 (MQ=255)

        cGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcc                                                                                                                                            >  2:2108841/1‑149 (MQ=255)

            gCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATg                                                                                                                                        >  2:1494586/1‑150 (MQ=255)

               caTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                     <  1:262869/150‑1 (MQ=255)

               caTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                     >  2:1456755/1‑150 (MQ=255)

                    cAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACt                                                                                                                                >  2:736087/1‑150 (MQ=255)

                       tCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGa                                                                                                                             <  2:565034/150‑1 (MQ=255)

                        cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                            <  2:963367/150‑1 (MQ=255)

                        cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                            <  2:539826/150‑1 (MQ=255)

                        cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                            <  1:809279/150‑1 (MQ=255)

                         tGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGATTGGATc                                                                                                                           >  2:791872/1‑150 (MQ=255)

                            gTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcc                                                                                                                        >  1:2343174/1‑150 (MQ=255)

                             tATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                                                                       >  1:2062226/1‑150 (MQ=255)

                              aTTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                                                                      >  1:1260673/1‑150 (MQ=255)

                                 gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                   <  2:1617776/150‑1 (MQ=255)

                                 gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                   >  1:358271/1‑150 (MQ=255)

                                  aaCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCa                                                                                                                  >  1:1432530/1‑150 (MQ=255)

                                    cAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                                                                >  1:124167/1‑150 (MQ=255)

                                       aTCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                                                                             >  2:2009876/1‑150 (MQ=255)

                                        tCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                                                                            <  1:1423312/150‑1 (MQ=255)

                                          cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                          <  1:1633847/150‑1 (MQ=255)

                                          cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                          >  2:374296/1‑150 (MQ=255)

                                          cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                          >  1:2188523/1‑150 (MQ=255)

                                           aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                         >  1:1050511/1‑150 (MQ=255)

                                             caacaaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGt                                                                                                       >  1:910949/1‑150 (MQ=255)

                                              aacaaGCGATTGCCGTTGGTTAAATAGTTTTCTGCTGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTa                                                                                                      <  1:1636604/150‑1 (MQ=255)

                                                  aGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCgg                                                                                                  >  1:2463989/1‑150 (MQ=255)

                                                    cGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCa                                                                                                <  2:2066874/150‑1 (MQ=255)

                                                    cGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGGCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCa                                                                                                <  2:1666567/150‑1 (MQ=255)

                                                     gATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAc                                                                                               <  1:376843/150‑1 (MQ=255)

                                                      aTTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGGCTGCCGCACATGGTTAGACGGGAACGCCCCGGCAGTAATCAGCGTATCGGCAcc                                                                                              >  1:615285/1‑150 (MQ=255)

                                                         gCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTc                                                                                           >  2:1084093/1‑150 (MQ=255)

                                                          ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCa                                                                                          <  2:276170/150‑1 (MQ=255)

                                                          ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCa                                                                                          <  2:1432530/150‑1 (MQ=255)

                                                          ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTACAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATNGGCACCTTCa                                                                                          <  1:1525252/150‑1 (MQ=255)

                                                              tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCTGCAGTAATCAGCGTATCGGCACCTTCACGCa                                                                                      <  2:466223/150‑1 (MQ=255)

                                                              tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCa                                                                                      >  1:1281330/1‑150 (MQ=255)

                                                              tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCACTCCAGCCCGAGTTTGGCAGCTACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCACGAATCAGCGTATCGGCACCTTCACGCa                                                                                      >  2:1180426/1‑150 (MQ=255)

                                                               ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCag                                                                                     <  1:2070370/150‑1 (MQ=255)

                                                                   aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                                                                 <  1:1009454/150‑1 (MQ=255)

                                                                   aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                                                                 <  1:229365/150‑1 (MQ=255)

                                                                     aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                               <  2:2093467/150‑1 (MQ=255)

                                                                       aGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                                                                             >  2:1889237/1‑150 (MQ=255)

                                                                        gTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGc                                                                            >  1:202320/1‑150 (MQ=255)

                                                                              tGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACTCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTTGATCGCCCCGTCAGTAATCAGAGTATCGGCACCTTCACGCAGAGCATCTGCCGCgag                                                                      <  2:2076310/150‑1 (MQ=255)

                                                                              tGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCgag                                                                      >  2:1837499/1‑150 (MQ=255)

                                                                                 ggTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaaa                                                                   >  1:1804381/1‑150 (MQ=255)

                                                                                  gTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAt                                                                  <  2:1541150/150‑1 (MQ=255)

                                                                                   ttGTGCCAATAGGATTTTCCAGCAGCGCCACTCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAtt                                                                 <  1:1906832/150‑1 (MQ=255)

                                                                                     gtgCCAATAGGATTTTCCAGGAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                               <  2:1662869/150‑1 (MQ=255)

                                                                                     gtgCCAATAGGATTTTCCAGCAGCGCCAGGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                               <  2:511781/150‑1 (MQ=255)

                                                                                     gagCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                               >  2:1247487/3‑150 (MQ=255)

                                                                                       gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg                                                             <  1:283160/150‑1 (MQ=255)

                                                                                        ccAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGc                                                            <  1:1604126/150‑1 (MQ=255)

                                                                                         cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                               <  2:824522/130‑1 (MQ=255)

                                                                                           aTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTa                                                         <  2:2463989/150‑1 (MQ=255)

                                                                                            tAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc                                                        <  1:326373/150‑1 (MQ=255)

                                                                                              ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGTCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt                                                      >  2:1493783/1‑150 (MQ=255)

                                                                                              ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt                                                      >  2:633488/1‑150 (MQ=255)

                                                                                               gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTa                                                     >  2:1962289/1‑150 (MQ=255)

                                                                                                aTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTaa                                                    >  1:908844/1‑150 (MQ=255)

                                                                                                aTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTaa                                                    <  1:908369/150‑1 (MQ=255)

                                                                                                     ccAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAt                                               <  2:887243/150‑1 (MQ=255)

                                                                                                        gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgc                                            <  1:1955586/150‑1 (MQ=255)

                                                                                                        gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgc                                            <  1:1889237/150‑1 (MQ=255)

                                                                                                        gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCCGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgc                                            >  2:2204451/1‑150 (MQ=255)

                                                                                                         cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgcc                                           >  1:734787/1‑150 (MQ=255)

                                                                                                            cgcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgccgcc                                        <  2:1494675/150‑1 (MQ=255)

                                                                                                             gcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGccc                                       >  1:2384886/1‑150 (MQ=255)

                                                                                                              ccACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCa                                      >  1:419009/1‑150 (MQ=255)

                                                                                                               cACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgccgc                                         <  1:1910950/146‑1 (MQ=255)

                                                                                                               cACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgccgc                                         >  2:1910950/1‑146 (MQ=255)

                                                                                                                  gCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGc                                  <  1:852986/150‑1 (MQ=255)

                                                                                                                   cAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGcc                                 >  1:1106373/1‑150 (MQ=255)

                                                                                                                     gTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                               >  1:404559/1‑150 (MQ=255)

                                                                                                                      tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATg                              <  2:583877/150‑1 (MQ=255)

                                                                                                                       gCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTCCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCAGTGCCATgg                             >  1:1163610/1‑150 (MQ=255)

                                                                                                                        cAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggg                            <  2:1599561/150‑1 (MQ=255)

                                                                                                                          gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggggg                          >  2:1311607/1‑150 (MQ=255)

                                                                                                                          gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggggg                          <  2:631104/150‑1 (MQ=255)

                                                                                                                          gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTGACGTAAATTATTGCCGCCCATTGCCAAggggg                          >  1:2221189/1‑150 (MQ=255)

                                                                                                                            ccGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTg                        >  2:1563503/1‑150 (MQ=255)

                                                                                                                             cGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGa                       <  1:1061306/150‑1 (MQ=255)

                                                                                                                               aGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGaag                     <  1:732691/150‑3 (MQ=255)

                                                                                                                               aGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACg                     <  1:1837499/150‑1 (MQ=255)

                                                                                                                                gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGt                    >  1:677182/1‑150 (MQ=255)

                                                                                                                                  ttGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCa                  >  2:1928168/1‑150 (MQ=255)

                                                                                                                                      cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc              >  2:1244292/1‑150 (MQ=255)

                                                                                                                                        ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGt            >  1:1245342/6‑150 (MQ=255)

                                                                                                                                          gACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGttt          <  1:2009876/150‑1 (MQ=255)

                                                                                                                                           tgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg         >  2:2091523/3‑150 (MQ=255)

                                                                                                                                           aCTGCGGCAGTCTGCCGCACATGGTTAGATTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg         <  1:791872/150‑1 (MQ=255)

                                                                                                                                            cTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtga        <  2:747075/150‑1 (MQ=255)

                                                                                                                                              gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatg      <  1:721259/150‑1 (MQ=255)

                                                                                                                                               cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGCGCATCTGCCGCGAGAAATTCAAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatga     >  2:385716/1‑150 (MQ=255)

                                                                                                                                               cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatga     >  2:1726434/1‑150 (MQ=255)

                                                                                                                                               cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatga     >  2:1332316/1‑150 (MQ=255)

                                                                                                                                               cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatga     >  1:834447/1‑150 (MQ=255)

                                                                                                                                                 gCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaa   >  1:2426861/1‑150 (MQ=255)

GGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAA  >  NC_000913/1999082‑1999377

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