BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,950,954	Δ1 bp	8.3%	coding (123/567 nt)	yecD →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,954	0	G	.	8.3%	427.2 / 30.8	132	coding (123/567 nt)	yecD	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (60/59); new base . (6/5); total (66/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

Only 100 of 135 total aligned reads displayed.

ACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGG  >  NC_000913/1950806‑1951103

aCGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                                                                       <  1:2106144/150‑1 (MQ=255)

   aTTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGa                                                                                                                                                    >  2:1670756/1‑150 (MQ=255)

   aTTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGa                                                                                                                                                    <  2:1868811/150‑1 (MQ=255)

    ttttGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGaa                                                                                                                                                   <  2:823047/150‑1 (MQ=255)

      ttGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAt                                                                                                                                                 <  2:2328355/150‑1 (MQ=255)

      ttGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAt                                                                                                                                                 >  1:2265621/1‑150 (MQ=255)

        gTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAttt                                                                                                                                               >  1:1223456/1‑150 (MQ=255)

        gTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGAGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAttt                                                                                                                                               >  2:2011964/1‑150 (MQ=255)

          tttttCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg                                                                                                                                             <  1:449736/150‑1 (MQ=255)

            tttCAGCATTCATTATGCTTGAACTTAATGCTAAAAACACTCCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑cg                                                                                                                                          >  2:584178/1‑150 (MQ=255)

             ttCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCg‑cgc                                                                                                                                         >  2:169109/1‑150 (MQ=255)

             ttCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCg‑cgc                                                                                                                                         <  1:1335415/150‑1 (MQ=255)

              tCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTcg‑ccgcc                                                                                                                                       <  2:804430/150‑5 (MQ=255)

              tCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTcg‑ccgcc                                                                                                                                       <  2:796048/150‑5 (MQ=255)

                 tgAGTCATTATGCTTTAACTTAATGCTAAAACCACTGATCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTTCCTGTTGTCAACTTACTGCCGATGAGGTAGTTAATCTCGCCGGTAAGCTGGC‑GAGAAATTTCGCCGCCAGc                                                                                                                                     <  2:1649927/149‑1 (MQ=255)

                 ggAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                                                                     <  2:112086/150‑1 (MQ=255)

                   aGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                                                                                   >  1:1470223/1‑150 (MQ=255)

                      cATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGTTCa                                                                                                                                <  2:659310/150‑1 (MQ=255)

                         tATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                                                                             <  1:839778/150‑1 (MQ=255)

                           tccTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCg                                                                                                                           >  1:440912/3‑150 (MQ=255)

                            gCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                                                          >  1:1587992/1‑150 (MQ=255)

                            gCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                                                          >  1:996952/1‑150 (MQ=255)

                            gCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                                                          <  1:134213/150‑1 (MQ=255)

                             cTTGCACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGACGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgtg                                                                                                                         <  1:2306768/150‑1 (MQ=255)

                                 aaCTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                                                                     <  2:1054252/150‑1 (MQ=255)

                                   cTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAATTCG‑CGCCAGCGGTCAGCCCGTGTTTCTg                                                                                                                   >  2:2162468/1‑150 (MQ=255)

                                    ttAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTgg                                                                                                                  <  2:1969647/150‑1 (MQ=255)

                                     tAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGt                                                                                                                 >  2:1613996/1‑150 (MQ=255)

                                     tAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGt                                                                                                                 <  1:1663926/150‑1 (MQ=255)

                                        tGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                                                                              <  2:1754406/150‑1 (MQ=255)

                                        tGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGACTTGAAGCTGGCTGCGAAATTTCG‑CGCCAGCTGTCAGCCCGTGTTTCTGGTgcg                                                                                                              <  2:1641505/150‑1 (MQ=255)

                                         gCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                                                                             >  2:1871839/1‑150 (MQ=255)

                                          cTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGcgcg                                                                                                            >  1:1235310/1‑150 (MQ=255)

                                            aaaaCCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                                                          >  2:1112510/1‑150 (MQ=255)

                                            aaaaCCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                                                          <  1:531723/150‑1 (MQ=255)

                                               aCCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGc                                                                                                       <  1:214285/150‑1 (MQ=255)

                                                   cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                                                                                   <  2:363706/150‑1 (MQ=255)

                                                    tGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                                                                  <  2:63119/150‑1 (MQ=255)

                                                       gcTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGc                                                                                               >  1:292560/1‑150 (MQ=255)

                                                         tggtggTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                                                             <  2:288437/150‑1 (MQ=255)

                                                         tggtggTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                                                             >  1:139831/1‑150 (MQ=255)

                                                          ggtggtGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                                                                            >  1:2204525/1‑150 (MQ=255)

                                                            tggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAtt                                                                                          <  1:1766557/150‑1 (MQ=255)

                                                             ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                                                         >  2:943392/1‑150 (MQ=255)

                                                             ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                                                         >  2:1532433/1‑150 (MQ=255)

                                                               tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                                                                       >  2:469492/1‑150 (MQ=255)

                                                                gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTTCCGATGAGGTAGTTATTCGCGCCGGGAAGCTGGTGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                                                                      <  1:675433/150‑1 (MQ=255)

                                                                gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                                                                      <  2:1523890/150‑1 (MQ=255)

                                                                gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                                                                      <  2:1243248/150‑1 (MQ=255)

                                                                    gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGaa                                                                                  >  1:1204829/1‑150 (MQ=255)

                                                                     aTTTACAAGAAGCCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTAGGCTGGTCTGCCGATTACGCTGAAg                                                                                 >  2:1959499/1‑150 (MQ=255)

                                                                      tttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                                                                <  1:938369/150‑1 (MQ=255)

                                                                      tttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                                                                >  1:827619/1‑150 (MQ=255)

                                                                        tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACAGCCGATGAGGTAGTTAATAGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGAGTTGGCTGGTCTGCCGATTACGCTGAAGCAt                                                                              >  1:1804116/1‑150 (MQ=255)

                                                                        tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACAGACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAt                                                                              >  1:106700/1‑150 (MQ=255)

                                                                         aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                                                             <  2:2445909/150‑1 (MQ=255)

                                                                          cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                                                            <  1:2169115/150‑1 (MQ=255)

                                                                            agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                                                                          <  2:2081298/150‑1 (MQ=255)

                                                                              aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAAc                                                                        >  1:58811/1‑150 (MQ=255)

                                                                               agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                                                                       <  1:1039915/150‑1 (MQ=255)

                                                                                     cTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGt                                                                 <  2:1908977/150‑1 (MQ=255)

                                                                                      ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc                                                                >  2:956015/1‑150 (MQ=255)

                                                                                       tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                                                               <  1:1895009/150‑1 (MQ=255)

                                                                                       tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                                                               <  1:1885869/150‑1 (MQ=255)

                                                                                               gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc                                                       <  1:469492/150‑1 (MQ=255)

                                                                                                 cGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTc                                                     <  1:1572297/150‑1 (MQ=255)

                                                                                                     ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCg                                                 >  2:1322293/1‑150 (MQ=255)

                                                                                                     ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCg                                                 >  2:41616/1‑150 (MQ=255)

                                                                                                     ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCg                                                 >  1:1491888/1‑150 (MQ=255)

                                                                                                       tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc                                               <  1:584178/150‑1 (MQ=255)

                                                                                                       tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc                                               <  1:1789461/150‑1 (MQ=255)

                                                                                                         cacaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTa                                             >  1:2043928/1‑150 (MQ=255)

                                                                                                          acaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTaa                                            <  2:1235310/150‑1 (MQ=255)

                                                                                                            aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg                                          <  1:2314775/150‑1 (MQ=255)

                                                                                                              aCTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAgtg                                        >  1:2480719/1‑150 (MQ=255)

                                                                                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGtt                                      >  1:2363339/1‑150 (MQ=255)

                                                                                                                 gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTg                                     >  2:335365/1‑150 (MQ=255)

                                                                                                                  ccGATGAGGTAGTTAATCGCGCCCGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGc                                    >  2:1875513/1‑150 (MQ=255)

                                                                                                                   cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGcc                                   <  2:1491888/150‑1 (MQ=255)

                                                                                                                        aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGaaa                              >  1:2415292/1‑150 (MQ=255)

                                                                                                                        aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGaaa                              >  1:750405/1‑150 (MQ=255)

                                                                                                                           tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAAtt                           >  1:544964/1‑150 (MQ=255)

                                                                                                                            aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑AGCCAGCGGTCAGACCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtg                          >  1:1643629/1‑150 (MQ=255)

                                                                                                                             gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                         <  1:2342458/150‑1 (MQ=255)

                                                                                                                              ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggt                        <  2:1388399/150‑1 (MQ=255)

                                                                                                                               tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                       >  2:107732/1‑150 (MQ=255)

                                                                                                                               tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                       <  2:457755/150‑1 (MQ=255)

                                                                                                                               tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                       >  2:1240128/1‑150 (MQ=255)

                                                                                                                                aaTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtgg                      <  1:856215/150‑1 (MQ=255)

                                                                                                                                 aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGc                     >  1:334212/1‑150 (MQ=255)

                                                                                                                                  tCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCa                    <  2:1466755/150‑1 (MQ=255)

                                                                                                                                   cgcgCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCaa                   <  2:139831/150‑1 (MQ=255)

                                                                                                                                      gcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACAt                <  1:1052438/150‑1 (MQ=255)

                                                                                                                                       ccGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATc               >  1:2460376/1‑150 (MQ=255)

                                                                                                                                           accgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgc           >  1:903392/4‑150 (MQ=255)

                                                                                                                                            ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgct          >  2:644230/3‑150 (MQ=255)

                                                                                                                                             cgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgctg         >  2:2158773/2‑150 (MQ=255)

                                                                                                                                              gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCtgctgc        >  1:2371624/1‑150 (MQ=255)

                                                                                                                                                tGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGAGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCAt      >  2:1450607/1‑150 (MQ=255)

                                                                                                                                                    ggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAgg  >  2:2232624/1‑150 (MQ=255)

ACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGG  >  NC_000913/1950806‑1951103

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