BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,265,885	A→G	S13G (AGC→GGC)	yeiR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,265,885	0	A	G	98.6%	460.7 / ‑2.9	146	S13G (AGC→GGC)	yeiR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (75/69); total (76/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 171 total aligned reads displayed.

TTTTTTAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCAGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAGTC  >  NC_000913/2265795‑2265974

ttttttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCggcggc                                                                                         <  1:1063619/93‑1 (MQ=255)

ttttttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCggcggc                                                                                         <  1:4624404/93‑1 (MQ=255)

ttttttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCggcggc                                                                                         <  1:1964076/93‑1 (MQ=255)

  ttttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTcggcggcgg                                                                                       <  2:2197199/93‑1 (MQ=255)

  ttttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTcggcggcgg                                                                                       >  2:4601190/1‑93 (MQ=255)

   tttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCggg                                                                                      <  1:646159/93‑1 (MQ=255)

   tttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCggg                                                                                      >  2:1272602/1‑93 (MQ=255)

   tttAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCggg                                                                                      <  1:3467155/93‑1 (MQ=255)

      aGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGaaa                                                                                   <  2:3639039/93‑1 (MQ=255)

        ccAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAAc                                                                                 >  2:47766/1‑93 (MQ=255)

         cAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAAcc                                                                                <  2:185404/93‑1 (MQ=255)

          aGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCa                                                                               <  1:933636/93‑1 (MQ=255)

          aGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCa                                                                               >  1:4731733/1‑93 (MQ=255)

          aGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCa                                                                               <  1:4441051/93‑1 (MQ=255)

          aGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCa                                                                               <  2:4051258/93‑1 (MQ=255)

           ggATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCAc                                                                              <  2:1468533/93‑1 (MQ=255)

           ggATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCAc                                                                              >  1:4336605/1‑93 (MQ=255)

            gATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACg                                                                             >  1:3019411/1‑93 (MQ=255)

              ttAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTc                                                                           <  1:3207045/93‑1 (MQ=255)

              ttAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTATGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTc                                                                           >  1:3343305/1‑93 (MQ=255)

                aaCGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGa                                                                         >  2:2580088/1‑93 (MQ=255)

                 aCGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAt                                                                        <  1:2791253/93‑1 (MQ=255)

                  cGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAtt                                                                       <  1:1404794/93‑1 (MQ=255)

                  cGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAtt                                                                       <  2:648314/93‑1 (MQ=255)

                    cTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAttct                                                                     <  1:2155694/93‑1 (MQ=255)

                    cTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAttct                                                                     >  2:163375/1‑93 (MQ=255)

                      cAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGAttcttc                                                                   >  1:348169/1‑93 (MQ=255)

                      cAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCCGGAAAACCACGTCGAttcttc                                                                   >  2:2087591/1‑93 (MQ=255)

                            aCTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGt                                                             >  1:2580903/1‑93 (MQ=255)

                             cTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCAGCGGGAAAACCACGTCGATTCTTCATCTGtt                                                            >  2:3677879/1‑93 (MQ=255)

                             actACCAGAAGTTTGGCGGCGCCGGTTTGACCAGGACCAACCTCACCACGGGTTTTCCCGCGGACGGGAAACCCACGTCGATTCTTCATCTGt                                                             <  2:33920/91‑1 (MQ=255)

                                aCCAG‑AGTATTGCAGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGc                                                         <  2:114859/93‑1 (MQ=255)

                                 ccAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGcc                                                        >  1:2064881/1‑93 (MQ=255)

                                  cAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGccc                                                       <  2:3521663/93‑1 (MQ=255)

                                   ag‑agTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCa                                                      >  2:3052152/1‑93 (MQ=255)

                                   ag‑agTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCa                                                      >  1:1011086/1‑93 (MQ=255)

                                   ag‑agTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCa                                                      <  1:119005/93‑1 (MQ=255)

                                      agTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATa                                                    <  1:234240/93‑1 (MQ=255)

                                       gTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAACTCCACGTCGATTCTTCATCTGTTAGCCCATaa                                                   <  1:1676990/93‑1 (MQ=255)

                                       gTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATaa                                                   >  1:1422340/1‑93 (MQ=255)

                                       gTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATaa                                                   >  1:1131444/1‑93 (MQ=255)

                                        tATTGGGCAGCCGTGGTGACCAGCACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATaaa                                                  <  2:2615171/93‑1 (MQ=255)

                                        tATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATaaa                                                  <  2:969023/93‑1 (MQ=255)

                                        tATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGAAAAACCACGTCGATTCTTCATCTGTTAGCCCATaaa                                                  <  2:2615170/93‑1 (MQ=255)

                                          ttGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGa                                                <  1:3881130/93‑1 (MQ=255)

                                          ttGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGa                                                >  1:3657327/1‑93 (MQ=255)

                                          ttGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGa                                                >  2:4414480/1‑93 (MQ=255)

                                          ttGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGa                                                <  1:624959/93‑1 (MQ=255)

                                           tGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGAt                                               >  2:1096249/1‑93 (MQ=255)

                                            gCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATc                                              >  1:2551476/1‑93 (MQ=255)

                                             cGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATcc                                             <  1:3648710/93‑1 (MQ=255)

                                              ggAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATccc                                            >  2:584702/1‑93 (MQ=255)

                                               gAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCa                                           >  2:1964953/1‑93 (MQ=255)

                                                aGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCaa                                          >  2:2916898/1‑93 (MQ=255)

                                                aGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTGCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCaa                                          <  1:1878409/93‑1 (MQ=255)

                                                   cGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGa                                       <  1:892202/93‑1 (MQ=255)

                                                   cGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAAACACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGa                                       >  2:2463027/1‑93 (MQ=255)

                                                    gTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGaa                                      >  1:3909259/1‑93 (MQ=255)

                                                    gTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGaa                                      >  2:1700103/1‑93 (MQ=255)

                                                    cttGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGaa                                      <  2:1639868/92‑1 (MQ=255)

                                                     ttGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGaaa                                     >  2:1409572/1‑93 (MQ=255)

                                                        tgACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAAt                                  <  1:3971492/93‑1 (MQ=255)

                                                        tgACCAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAAt                                  >  1:2066360/1‑93 (MQ=255)

                                                           ccAGGACCAACCTCATCCCCGGTTTTCTCCGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATggg                               <  2:2349458/93‑1 (MQ=255)

                                                           ccAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATggg                               <  2:2398298/93‑1 (MQ=255)

                                                           ccAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATggg                               <  2:3242760/93‑1 (MQ=255)

                                                           ccAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATggg                               <  1:1011440/93‑1 (MQ=255)

                                                            cAGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGc                              <  2:3668188/93‑1 (MQ=255)

                                                             aGGACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCg                             <  1:3521117/93‑1 (MQ=255)

                                                              ggACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCgg                            >  2:3292700/1‑93 (MQ=255)

                                                               gACCAACCTCATCACCGGTTTTCTCGGCGGCTGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGt                           <  2:3188236/93‑1 (MQ=255)

                                                               gACCAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGt                           <  1:771312/93‑1 (MQ=255)

                                                               gACCAACCTCATCACCGGTTTTCTCGGCAGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGt                           <  1:3133960/93‑1 (MQ=255)

                                                                 ccAACCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTAc                         <  1:4103928/93‑1 (MQ=255)

                                                                  cAACCTCATCACGGTTTTTCTCGGCGGCGGGAAAACCCCGTCGATTCTTCATCTGTTCGCCCATAAAGATAACAACGAAAAATGGGCGGTACt                        <  1:842753/93‑1 (MQ=255)

                                                                   acccTCATCCCCGTTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTg                       <  2:929982/91‑1 (MQ=255)

                                                                   aaCCTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTg                       >  1:2890489/1‑93 (MQ=255)

                                                                      cTCATCACCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGtt                    >  1:4405312/1‑93 (MQ=255)

                                                                       tcatcaCCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTa                   >  1:3075198/1‑93 (MQ=255)

                                                                        cctcaCCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTGTTCATCTGTTAGGCCATAAAGATCCCAACGAAAATTGGGCGGTACTGGTTaa                  >  2:1339913/3‑93 (MQ=255)

                                                                         ctcaCGGGTTTTATCGGCGGCACGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGACCCCAACGAAACATGGGCGGTACTGGTTAAt                 <  2:559992/92‑1 (MQ=255)

                                                                         atcaCCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAAt                 >  1:4278783/1‑93 (MQ=255)

                                                                          tcaCCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATg                <  1:32131/93‑1 (MQ=255)

                                                                           caCCGGTGTTCTCGGCGGCGGGAAATCCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGa               <  1:1308846/93‑1 (MQ=255)

                                                                            aCCGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGaa              >  1:2504475/1‑93 (MQ=255)

                                                                             ccGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAAt             >  2:1157014/1‑93 (MQ=255)

                                                                              cGGTTTTCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAAtt            >  1:3621436/1‑93 (MQ=255)

                                                                                 ttttCTCGGCGGCGTGAAAGCTACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTgg         <  2:426320/93‑1 (MQ=255)

                                                                                 ttttCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTgg         <  2:426319/93‑1 (MQ=255)

                                                                                  tttCTCGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTggg        <  2:513089/93‑1 (MQ=255)

                                                                                    tctcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGa      >  1:3454594/1‑93 (MQ=255)

                                                                                    tctcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGa      <  2:3675739/93‑1 (MQ=255)

                                                                                    gctcGACGGCGGCAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAATATGGGCGGTACTGGTTAATGAATTTGGGGa      <  2:1988030/92‑1 (MQ=255)

                                                                                     ctcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGATAATGAATTTGGGGaa     <  1:2753961/93‑1 (MQ=255)

                                                                                      tcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGGCCATAAAGATCCCAACGAAAAATAGGCGGTACTGGTTAATGAATTTGGGGAAg    >  2:1913015/1‑93 (MQ=255)

                                                                                      tcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAg    <  1:3123082/93‑1 (MQ=255)

                                                                                      tcGGCGGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTAGGGAAg    >  2:1913016/1‑93 (MQ=255)

                                                                                       cggcggcggGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAGt   <  1:1586097/93‑1 (MQ=255)

                                                                                        ggcggcGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAGTc  >  1:385502/1‑93 (MQ=255)

                                                                                        ggcggcGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAGTc  <  1:1999330/93‑1 (MQ=255)

TTTTTTAGCCAGGATTAACGCTCAGTTAACTTACCAG‑AGTATTGCGGAGCCGTTGTGACCAGGACCAACCTCATCACCGGTTTTCTCGGCAGCGGGAAAACCACGTCGATTCTTCATCTGTTAGCCCATAAAGATCCCAACGAAAAATGGGCGGTACTGGTTAATGAATTTGGGGAAGTC  >  NC_000913/2265795‑2265974

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