BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,349,053	A→T	7.0%	I352N (ATT→AAT)	yciW ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,349,053	0	A	T	7.0%	289.3 / 6.1	101	I352N (ATT→AAT)	yciW	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (43/50); new base T (3/5); total (46/55)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCG  >  NC_000913/1348904‑1349202

cTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaa                                                                                                                                                       >  2:663949/1‑150 (MQ=255)

 tGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaat                                                                                                                                                      >  1:1099946/1‑150 (MQ=255)

  gAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaata                                                                                                                                                     >  2:1431839/1‑150 (MQ=255)

    aaCATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTaataat                                                                                                                                                   <  2:407332/150‑1 (MQ=255)

         aaTTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCg                                                                                                                                              <  1:706821/150‑1 (MQ=255)

           ttGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTAATCTGctgg                                                                                                                                            <  1:650801/150‑4 (MQ=255)

            tGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCAATACCAGCAGATTAATAATCTGCGTtg                                                                                                                                           >  2:1504436/1‑150 (MQ=255)

            tGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTATTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtg                                                                                                                                           <  1:113515/150‑1 (MQ=255)

               cTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTAATCTGctggtatg                                                                                                                                        <  2:762246/150‑8 (MQ=255)

               cTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtgatg                                                                                                                                        <  2:582934/150‑1 (MQ=255)

               cTAATTTGGTGGCCAATAGAGCGGCGGACAGTGTTACGAGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCAGGCAACGCGATTTTTTAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtgatg                                                                                                                                        <  1:157390/150‑1 (MQ=255)

                   tttGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGAttt                                                                                                                                    >  2:1117488/1‑150 (MQ=255)

                      ggtgtCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTc                                                                                                                                 <  1:738368/150‑1 (MQ=255)

                      ggtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTc                                                                                                                                 <  2:1733423/150‑1 (MQ=255)

                       gtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCg                                                                                                                                <  1:485674/150‑1 (MQ=255)

                        tggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCgt                                                                                                                               >  2:1290441/1‑150 (MQ=255)

                        tggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCgt                                                                                                                               >  1:267709/1‑150 (MQ=255)

                         ggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGAATTTACgtg                                                                                                                              >  1:261056/1‑150 (MQ=255)

                         ggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGAAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCgtg                                                                                                                              <  1:258709/150‑1 (MQ=255)

                          gCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAAAAATAATCTGCGTTGATGATGTTTCgtgt                                                                                                                             >  2:1165146/1‑150 (MQ=255)

                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGtt                                                                                                                            <  1:1747493/150‑1 (MQ=255)

                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGtt                                                                                                                            >  1:1292008/1‑150 (MQ=255)

                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGtt                                                                                                                            <  1:646202/150‑1 (MQ=255)

                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGtt                                                                                                                            >  1:1483387/1‑150 (MQ=255)

                            cATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTc                                                                                                                           <  1:1661893/150‑1 (MQ=255)

                            aataagagCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGGTc                                                                                                                           <  1:1661895/146‑1 (MQ=255)

                              ttAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGa                                                                                                                         <  1:1610771/150‑1 (MQ=255)

                               tAGAGCGGCTGACAGTTTTACGCGAATCTGGCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAg                                                                                                                        <  1:1184816/150‑1 (MQ=255)

                                       cTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGc                                                                                                                <  1:664965/150‑1 (MQ=255)

                                           cAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCCCCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGaa                                                                                                            >  1:1333331/1‑150 (MQ=255)

                                            aGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAAc                                                                                                           <  1:557443/150‑1 (MQ=255)

                                             gTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTGTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACt                                                                                                          >  2:1877925/1‑150 (MQ=255)

                                             gTTTAACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCACACCCAACCCCATATTTAAGCGATTTATACAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCAGCCTTGATGATTTATCGAGTTCGAGTAAAGGGATCAACt                                                                                                          >  2:505322/1‑150 (MQ=255)

                                                  aCGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGggcg                                                                                                     >  1:1406239/1‑150 (MQ=255)

                                                  aCGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAATAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGggcg                                                                                                     <  1:260179/150‑1 (MQ=255)

                                                      gAATCTGTCTGACGCGGAAAGGTTAATATTTCTCAACAAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCATCAGATAAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTTGGCGgcgc                                                                                                 <  2:906524/150‑1 (MQ=255)

                                                         tctgtctgACGCGGCAAGGTTAATATGTCTCATCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGCTTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTaa                                                                                              <  1:308042/150‑1 (MQ=255)

                                                              ctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGa                                                                                         >  1:1023678/1‑150 (MQ=255)

                                                                gACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATc                                                                                       <  1:933055/150‑1 (MQ=255)

                                                                gACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATc                                                                                       >  1:1475949/1‑150 (MQ=255)

                                                                 ccgcgGCAAGGTTAATATGGCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTTATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCg                                                                                      <  1:1252950/149‑1 (MQ=255)

                                                                  cgcgGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCgg                                                                                     >  1:1245321/1‑150 (MQ=255)

                                                                   gcgGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCggg                                                                                    <  1:817270/150‑1 (MQ=255)

                                                                   gcgGCAAGGGTAATATGACTCACCCACCGCGAATTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATAAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCggg                                                                                    <  1:144776/150‑1 (MQ=255)

                                                                    cgGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGt                                                                                   <  1:1637489/150‑1 (MQ=255)

                                                                     ggCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTg                                                                                  >  2:1320173/1‑150 (MQ=255)

                                                                       cAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTACGAGTAAAGTGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCt                                                                                >  2:902977/1‑150 (MQ=255)

                                                                          ggTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGTGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCgg                                                                             >  1:1204517/1‑150 (MQ=255)

                                                                           gTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCggg                                                                            <  1:952128/150‑1 (MQ=255)

                                                                           gTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCggg                                                                            >  1:725263/1‑150 (MQ=255)

                                                                             tAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTc                                                                          >  2:1307072/1‑150 (MQ=255)

                                                                              aaTATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCa                                                                         <  1:1021745/150‑1 (MQ=255)

                                                                                tatGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCAGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGc                                                                       <  2:1798076/150‑1 (MQ=255)

                                                                                    tctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCAtt                                                                   <  2:1781778/150‑1 (MQ=255)

                                                                                     ctcACCCAACGCGTTTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTATTCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGATCTGGGCGGGGCTAAAGCGATGGGGTGCTCGGGTCAGCCATTg                                                                  <  1:1504436/150‑1 (MQ=255)

                                                                                     ctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTg                                                                  >  1:1130542/1‑150 (MQ=255)

                                                                                     ctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTg                                                                  <  2:678295/150‑1 (MQ=255)

                                                                                      tcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGAAGATTAATAATCTGCGTTGAGGATTTTTAGTGTTAGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGt                                                                 >  2:323499/1‑150 (MQ=255)

                                                                                          ccAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCg                                                             >  1:1723069/1‑150 (MQ=255)

                                                                                           cAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGc                                                            >  1:1046940/1‑150 (MQ=255)

                                                                                             aCGCTATTATTCAGCGATTTATCCAGCCACACATCCCGCTGCATACCAGCAGATAAATAATCTGCGTTGATGATTTTTCGTGTTTGAGTAAAGGTCTGAACTGGGCGGCGCAAAAGCGATCGGGTGCCCGGGTCAGCCATTGTATCGCCt                                                          <  1:1212553/150‑1 (MQ=255)

                                                                                              cgcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTg                                                         <  2:1873202/150‑1 (MQ=255)

                                                                                                 gATTTTTCAGCGATCTCTCCAGCCCCACAGCCCGCTCCCTCCCCGCAGCTGAATAAGCTGCGATGATGATTTTTGGCGTGCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAAc                                                      <  1:306139/150‑1 (MQ=255)

                                                                                                   tttttAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGAGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACgg                                                    >  1:898639/1‑150 (MQ=255)

                                                                                                         aGCGATTTATCCAGCCACACAGCCCGCTCCATACCCGCAGATTAATAATCTGCGTTGATGATTTTTCGTGGGCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGACCGGTAgcgc                                              <  1:1565064/150‑1 (MQ=255)

                                                                                                          gCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAgcgcg                                             >  2:1594834/1‑150 (MQ=255)

                                                                                                           cGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAgcgcgc                                            <  1:1292504/150‑1 (MQ=255)

                                                                                                            gATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCt                                           <  1:543653/150‑1 (MQ=255)

                                                                                                               ttATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCaa                                        <  1:647442/150‑1 (MQ=255)

                                                                                                                tATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAAc                                       <  2:1118134/150‑1 (MQ=255)

                                                                                                                   ccAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGCACTGGGCGGCGCTTAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCGGAACGGTAGCGCGCTCAACAgc                                    >  2:1410241/1‑149 (MQ=255)

                                                                                                                   ccAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGa                                    >  1:1846278/1‑150 (MQ=255)

                                                                                                                    cAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGCGCGGGTGCTCGGGTCAGCCATTGTCTCGCGTGAACGGTAGCGCGCTCAACAGaa                                   >  2:1642978/1‑150 (MQ=255)

                                                                                                                    cAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGaa                                   <  1:1431839/150‑1 (MQ=255)

                                                                                                                      gCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATg                                 <  2:1569224/150‑1 (MQ=255)

                                                                                                                       ccACACAGCCCTCTCCATACCAGCATATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGGCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGa                                <  2:261056/150‑1 (MQ=255)

                                                                                                                       ccACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGa                                >  2:547977/1‑150 (MQ=255)

                                                                                                                        cacacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGAc                               <  2:791943/150‑1 (MQ=255)

                                                                                                                          cacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAATGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACgg                             <  2:1475949/150‑1 (MQ=255)

                                                                                                                          cacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACgg                             >  2:1831064/1‑150 (MQ=255)

                                                                                                                          cacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACgg                             >  1:955645/1‑150 (MQ=255)

                                                                                                                          cacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTAGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACgg                             >  1:23540/1‑150 (MQ=255)

                                                                                                                          cacaGCCCGATCCATACCAGCAGATTAATAAACTGCGTTGATGATTTTTCGTGATCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCAATGTATCGCCTGAACGGTAGCGCGCTCCACCGAATGACgg                             >  2:1596381/1‑150 (MQ=255)

                                                                                                                            caGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGCACTGGGCGGCTCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTAAACAGAATGACGGTc                           >  1:1365752/1‑150 (MQ=255)

                                                                                                                            caGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTc                           >  1:1638598/1‑150 (MQ=255)

                                                                                                                            caGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTc                           <  2:714998/150‑1 (MQ=255)

                                                                                                                              gCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAg                         >  1:356368/1‑150 (MQ=255)

                                                                                                                              gCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGAGTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAg                         >  2:284788/1‑150 (MQ=255)

                                                                                                                                ccGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGcc                       <  1:1729742/150‑1 (MQ=255)

                                                                                                                                     ccATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGACCGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTg                  >  1:1477059/1‑150 (MQ=255)

                                                                                                                                        tcaacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGtt               >  1:1572897/6‑150 (MQ=255)

                                                                                                                                          ccAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCAATTGGTTat             >  1:1594210/1‑149 (MQ=255)

                                                                                                                                          aacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTaa             >  1:1708711/4‑150 (MQ=255)

                                                                                                                                          aacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTaa             >  2:1064374/4‑150 (MQ=255)

                                                                                                                                             gcagATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAAttt          <  1:1575636/150‑1 (MQ=255)

                                                                                                                                               agATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGTGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCa        <  2:1483387/150‑1 (MQ=255)

                                                                                                                                                  ttAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGTGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGACATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCAttt     >  1:1402527/1‑150 (MQ=255)

                                                                                                                                                    aataatCTGCGTTGATGATTTTTCGTGTTCGAGTATAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTCTCGCCTGAACGGTAGCGCGCTAAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTc   >  2:87378/1‑150 (MQ=255)

                                                                                                                                                     ataatCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCg  <  1:1320173/150‑1 (MQ=255)

                                                                                                                                                     ataatCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCg  <  1:830274/150‑1 (MQ=255)

CTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCG  >  NC_000913/1348904‑1349202

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