BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,586,311	T→A	75.0%	H59L (CAC→CTC)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,586,311	0	T	A	75.0%	212.4 / 82.3	124	H59L (CAC→CTC)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (13/18); new base A (47/46); total (60/64)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 125 total aligned reads displayed.

CGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTA  >  NC_000913/1586214‑1586409

cGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGc                                                                                                  <  1:2669094/100‑1 (MQ=255)

 gTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCt                                                                                                 <  1:4021782/100‑1 (MQ=255)

   aTTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTa                                                                                               >  1:3139458/1‑100 (MQ=255)

   aTTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTa                                                                                               >  2:2135359/1‑100 (MQ=255)

    ttCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAg                                                                                              >  1:2452521/1‑100 (MQ=255)

    ttCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAg                                                                                              >  1:211188/1‑100 (MQ=255)

     tctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAgg                                                                                             <  1:901738/100‑1 (MQ=255)

      ctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGa                                                                                            <  2:903169/100‑1 (MQ=255)

        cAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATc                                                                                          >  2:3041232/1‑100 (MQ=255)

         ngCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCt                                                                                         >  2:966238/2‑100 (MQ=255)

         aGCAAAGAAAGAGGCGTTTACGTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCTAATGACGCACTGAGCTTAGGAGCt                                                                                         <  2:1882438/100‑1 (MQ=255)

         aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCt                                                                                         >  1:1129676/1‑100 (MQ=255)

          gCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTg                                                                                        >  1:3572131/1‑100 (MQ=255)

           cAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTgg                                                                                       >  1:3579620/1‑100 (MQ=255)

           cAAAGAAAGAGGCGTTTACCTGCTTATCCCTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTgg                                                                                       <  1:4175307/100‑1 (MQ=255)

            aaagaaagAGGCGTTTACTTGCTTATTCGTTAGTTCCCTGGCGATTGCTTTTGCGCCGTTTTCACATTTGTCAAATGACGCACTGTGCTTTGGATCTTGc                                                                                      <  2:1129017/100‑1 (MQ=255)

            aaagaaagAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                                                                      >  2:1341068/1‑100 (MQ=255)

            aaagaaagAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                                                                      >  2:2302977/1‑100 (MQ=255)

            aaagaaagAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                                                                      <  1:238501/100‑1 (MQ=255)

             aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGcc                                                                                     <  1:1063227/100‑1 (MQ=255)

              agaaagaGGCGTTTACCTGCTTATCCATGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCa                                                                                    <  2:609975/100‑1 (MQ=255)

                aaagaGGCGTTTACCTGCTAATCCGTGAATTCCCATGCGATTTCTTTTTAGCCTTTTGCACAAATGACGAATGAGTCTCTGAGCTTAGTATCTCTCCATg                                                                                  >  2:1825141/1‑100 (MQ=255)

                 aagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGc                                                                                 <  2:2282778/100‑1 (MQ=255)

                   gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGcac                                                                               <  2:3950715/100‑1 (MQ=255)

                   cagccGTGTACCTGAGGATCCGAGACTTCCCATGCTATTTCTGTTGCGCCGTTTTCACAAATTTCTAATGACGCACTTAGCTGAGGATCTGGCCATGcac                                                                               >  2:3479080/5‑100 (MQ=255)

                       cGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAAcc                                                                           >  1:1220142/1‑100 (MQ=255)

                        gTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAAccc                                                                          >  2:2834108/1‑100 (MQ=255)

                          ttACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCgg                                                                        <  1:1028660/100‑1 (MQ=255)

                           tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGa                                                                       >  1:366956/1‑100 (MQ=255)

                            aacTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGAc                                                                      >  1:2815503/3‑100 (MQ=255)

                            aCCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGAc                                                                      <  1:3495865/100‑1 (MQ=255)

                             ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     >  2:1953538/1‑100 (MQ=255)

                             ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     >  1:1136264/1‑100 (MQ=255)

                             ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     <  1:2724777/100‑1 (MQ=255)

                              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     >  2:593705/1‑99 (MQ=255)

                              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAg                                                                    <  1:230963/100‑1 (MQ=255)

                              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAg                                                                    <  2:3253311/100‑1 (MQ=255)

                               tGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGt                                                                   >  2:2144907/1‑100 (MQ=255)

                                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  <  1:1940878/100‑1 (MQ=255)

                                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  <  2:3387716/100‑1 (MQ=255)

                                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  >  2:4166887/1‑100 (MQ=255)

                                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  <  1:1562624/100‑1 (MQ=255)

                                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  >  2:3656737/1‑100 (MQ=255)

                                    aTCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAg                                                              <  2:3151971/100‑1 (MQ=255)

                                     tCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGc                                                             >  2:604851/1‑100 (MQ=255)

                                      ccGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGcc                                                            >  2:420717/1‑100 (MQ=255)

                                       cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCATAGGATCTGGCCATGCACAACCCGGCCAGTCAAAGcc                                                            >  2:2262641/1‑99 (MQ=255)

                                         tGACTTCCCAGGCGATTGCTTTGGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCctc                                                         <  2:1210396/100‑1 (MQ=255)

                                           aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  2:738623/100‑1 (MQ=255)

                                           aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  1:3656764/100‑1 (MQ=255)

                                            cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTgg                                                      >  2:1936621/1‑100 (MQ=255)

                                            cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTgg                                                      >  1:380/1‑100 (MQ=255)

                                            cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTgg                                                      >  1:2435816/1‑100 (MQ=255)

                                            cTTCCCAGGAGATTGCTTTTGCGCTGTTTTCACAAGTGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCctct                                                        <  2:1778142/98‑1 (MQ=255)

                                              tCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCt                                                    <  2:1981436/100‑1 (MQ=255)

                                              tCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCt                                                    >  1:1249663/1‑100 (MQ=255)

                                               cccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCtt                                                   >  1:3562755/1‑100 (MQ=255)

                                               cccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCtt                                                   >  1:1295860/1‑100 (MQ=255)

                                                ccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTa                                                  <  1:1702596/100‑1 (MQ=255)

                                                ccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTa                                                  >  1:3442978/1‑100 (MQ=255)

                                                 cAGGCGATTTCTTTTGCTCTTTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTAt                                                 >  2:1706190/1‑100 (MQ=255)

                                                 cAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTAt                                                 <  1:2623666/100‑1 (MQ=255)

                                                    gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCa                                              <  1:4039458/100‑1 (MQ=255)

                                                    gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCa                                              <  1:1590752/100‑1 (MQ=255)

                                                      gATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCCTACACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATg                                            <  1:2911199/100‑1 (MQ=255)

                                                       aTTTCTTTTGCGCGGTTTTCAGAAATGTCGAAGGACGCAGTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATg                                            <  2:2510228/99‑1 (MQ=255)

                                                         tGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGGTCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCa                                         <  1:3905174/100‑1 (MQ=255)

                                                           cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCaaa                                       >  2:994015/1‑100 (MQ=255)

                                                             tttGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACa                                     <  2:3779462/100‑1 (MQ=255)

                                                                gcgcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATg                                   <  2:1982673/99‑1 (MQ=255)

                                                                 cgcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGc                                 <  1:165814/100‑1 (MQ=255)

                                                                 cgcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGc                                 >  1:515721/1‑100 (MQ=255)

                                                                  gcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCa                                <  1:2600848/100‑1 (MQ=255)

                                                                  gcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCa                                >  1:2256266/1‑100 (MQ=255)

                                                                   ccGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCaa                               >  1:455870/1‑100 (MQ=255)

                                                                    cGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAAt                              <  2:3474308/100‑1 (MQ=255)

                                                                     gTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATa                             <  2:3788998/100‑1 (MQ=255)

                                                                       tttCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAAc                           <  1:195144/100‑1 (MQ=255)

                                                                        ttCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACa                          >  1:2720274/1‑100 (MQ=255)

                                                                           acaAATGTCGAATGACGCACTGAGCTTAGGATCAGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCGTATTCATGTCAAACATGGCAATAACATcc                       >  2:2506681/1‑100 (MQ=255)

                                                                            caAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCt                      >  2:286083/1‑100 (MQ=255)

                                                                            caAATGTCGAATGACGCACTGAGCTTAGGATCTCGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATGTATGTCAAACATGGGAATAACATCCt                      >  2:2856998/1‑100 (MQ=255)

                                                                            caAATGTAGAATGACGCACAGTGCTTAGGATCTGGCCATGCACAACAAGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCt                      >  2:995339/1‑100 (MQ=255)

                                                                               aTGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGc                   <  2:3695707/100‑1 (MQ=255)

                                                                                tGTCGACTGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  2:2823793/100‑1 (MQ=255)

                                                                                tGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  1:4081464/100‑1 (MQ=255)

                                                                                   cGAATGACGCACTGAGCTTAGTATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGtat               <  2:1021436/100‑1 (MQ=255)

                                                                                   cGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGtat               >  1:3211338/1‑100 (MQ=255)

                                                                                    gAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGtata              <  2:3377119/100‑1 (MQ=255)

                                                                                    gAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGtata              >  2:3994599/1‑100 (MQ=255)

                                                                                     aaTGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTatat             <  2:1325209/100‑1 (MQ=255)

                                                                                       tGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATcc           >  1:3676622/1‑100 (MQ=255)

                                                                                        gACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTCTATCCa          >  2:2703735/1‑100 (MQ=255)

                                                                                          cGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATc        <  2:1551255/100‑1 (MQ=255)

                                                                                           gCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCt       <  2:1126076/100‑1 (MQ=255)

                                                                                            cACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTg      >  1:1502546/1‑100 (MQ=255)

                                                                                            cACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTg      <  2:3153450/100‑1 (MQ=255)

                                                                                             aCTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGc     <  1:605181/100‑1 (MQ=255)

                                                                                              cTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCt    >  1:3338099/1‑100 (MQ=255)

                                                                                                gtgCTTAGGATCTGGCCATGCACAACCCGGACAGTCGAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTa  <  1:2925524/100‑1 (MQ=255)

CGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTA  >  NC_000913/1586214‑1586409

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