BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	908,737	T→C	Q512Q (CAA→CAG)	poxB ←	pyruvate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	908,737	0	T	C	99.6%	696.6 / ‑5.7	223	Q512Q (CAA→CAG)	poxB	pyruvate oxidase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (102/120); total (102/121)
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.95e-01

Only 100 of 227 total aligned reads displayed.

TGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTG  >  NC_000913/908648‑908826

tGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGc                                                                                        <  2:1406903/93‑1 (MQ=255)

 gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                       <  1:6063687/93‑1 (MQ=255)

 gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                       >  2:3558381/1‑93 (MQ=255)

 gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                       >  1:3336700/1‑93 (MQ=255)

 gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                       <  1:6792911/93‑1 (MQ=255)

  ttCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAg                                                                                      >  1:5331927/1‑93 (MQ=255)

  ttCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAg                                                                                      <  2:4581725/93‑1 (MQ=255)

   gcGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCACCACCGGACCGTCGATGGAGAAGGCGCGCTGCAgg                                                                                     <  2:6054273/92‑1 (MQ=255)

     gAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGc                                                                                   <  2:6822651/93‑1 (MQ=255)

     gAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGc                                                                                   >  1:2002256/1‑93 (MQ=255)

      aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  >  1:7367566/1‑93 (MQ=255)

      aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  <  1:1598705/93‑1 (MQ=255)

      aGTTTGATCTGCGGTGGAATGGCTAACTCTTCCTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  >  2:4435611/1‑93 (MQ=255)

        tttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtc                                                                                >  2:7678402/1‑93 (MQ=255)

        tttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtc                                                                                >  2:437101/1‑93 (MQ=255)

         ttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtca                                                                               <  2:5971960/93‑1 (MQ=255)

           gATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatc                                                                             >  1:6924534/1‑93 (MQ=255)

           gATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatc                                                                             >  2:1788364/1‑93 (MQ=255)

           gATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatc                                                                             <  2:2251766/93‑1 (MQ=255)

           gATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatc                                                                             <  1:5965456/93‑1 (MQ=255)

            aTCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatca                                                                            >  2:6119649/1‑93 (MQ=255)

            aTCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatca                                                                            <  1:5623425/93‑1 (MQ=255)

             tCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCaa                                                                           >  1:3113310/1‑93 (MQ=255)

             tCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCaa                                                                           >  1:3360652/1‑93 (MQ=255)

              cTGCGGTGGAATGGCTACCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                          <  2:3378809/93‑1 (MQ=255)

              cTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                          <  1:5179626/93‑1 (MQ=255)

                  gggggAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      <  1:4565727/90‑1 (MQ=255)

                  ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      >  2:4230740/1‑93 (MQ=255)

                  ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      <  1:2705924/93‑1 (MQ=255)

                  ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      <  1:3594904/93‑1 (MQ=255)

                    tGGAATGGCTTACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGa                                                                    >  2:5092359/1‑93 (MQ=255)

                    tGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGa                                                                    >  2:6342713/1‑93 (MQ=255)

                         tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                               >  1:4548966/1‑93 (MQ=255)

                          ggCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                              >  2:7934061/1‑93 (MQ=255)

                          ggCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                              >  1:5025453/1‑93 (MQ=255)

                            cTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttttt                                                            >  1:1042446/1‑93 (MQ=255)

                            cTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttttt                                                            >  2:2693904/1‑93 (MQ=255)

                             tAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTc                                                           >  1:1179541/1‑93 (MQ=255)

                             tAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTc                                                           <  1:4088764/93‑1 (MQ=255)

                              aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                          >  2:1639216/1‑93 (MQ=255)

                              aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                          >  1:3856401/1‑93 (MQ=255)

                               aCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                         <  2:4676380/93‑1 (MQ=255)

                               aCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                         >  1:5570800/1‑93 (MQ=255)

                               aCTCTTCTTTGGCGACCACCACAACCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                         <  1:7018264/93‑1 (MQ=255)

                                 tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                       >  1:2234604/1‑93 (MQ=255)

                                 tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                       >  1:2385634/1‑93 (MQ=255)

                                  cttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                                                      <  1:4929386/93‑1 (MQ=255)

                                   ttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGTTTCATCAACTTCTGTCGCTTTTTCTACACg                                                     >  2:5811118/1‑93 (MQ=255)

                                   ttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACg                                                     >  1:7164010/1‑93 (MQ=255)

                                   ttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACg                                                     <  1:1924556/93‑1 (MQ=255)

                                      tttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGAt                                                  >  1:3244292/1‑93 (MQ=255)

                                       ttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATa                                                 <  2:4352199/93‑1 (MQ=255)

                                       ttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTAAACGGATa                                                 >  1:2968743/1‑93 (MQ=255)

                                        tGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAc                                                <  2:2635721/93‑1 (MQ=255)

                                        tGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAc                                                >  2:4722373/1‑93 (MQ=255)

                                        tGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAc                                                <  1:5713886/93‑1 (MQ=255)

                                         ggCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAcc                                               <  2:5660560/93‑1 (MQ=255)

                                         ggCGACCACCACATCAACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATaac                                               >  2:1817686/1‑91 (MQ=255)

                                          gCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAccc                                              <  2:376314/93‑1 (MQ=255)

                                            gACCACCACATCCGCCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGt                                            <  1:3569809/93‑1 (MQ=255)

                                             accaccACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTa                                           >  1:2902875/1‑93 (MQ=255)

                                              ccaccaCATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTaa                                          <  2:1590064/93‑1 (MQ=255)

                                               caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                         >  2:4989545/1‑93 (MQ=255)

                                               caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                         >  2:4523740/1‑93 (MQ=255)

                                                cccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATg                                        <  2:5744543/92‑1 (MQ=255)

                                                accacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATg                                        >  1:4283318/1‑93 (MQ=255)

                                                accacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATg                                        >  2:6489355/1‑93 (MQ=255)

                                                 ccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGc                                       >  1:4747437/1‑93 (MQ=255)

                                                 ccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGc                                       >  1:1612551/1‑93 (MQ=255)

                                                  cacaTCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGcc                                      >  2:1971520/1‑93 (MQ=255)

                                                     aTCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGcg                                   >  2:2853963/1‑92 (MQ=255)

                                                      tCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGcgc                                  >  2:1024187/1‑91 (MQ=255)

                                                          ccaACACCGGACCGTCGATGGAGAAGGGGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCtt                              >  2:3455837/1‑93 (MQ=255)

                                                           caACACCGGACCGTCGATGGAGAAGGGGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTc                             >  2:1128726/1‑93 (MQ=255)

                                                           caACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTc                             <  2:1383438/93‑1 (MQ=255)

                                                              cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGATTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCgga                          >  2:7144116/1‑92 (MQ=255)

                                                               acCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCa                         >  2:4783757/1‑93 (MQ=255)

                                                                ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                        >  1:4628873/1‑93 (MQ=255)

                                                                ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                        <  1:6029522/93‑1 (MQ=255)

                                                                ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                        >  1:7645541/1‑93 (MQ=255)

                                                                 cGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAAt                       >  1:1815153/1‑93 (MQ=255)

                                                                   gACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGc                     >  2:3601185/1‑93 (MQ=255)

                                                                   gACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGc                     >  2:3366336/1‑93 (MQ=255)

                                                                    gccGTCGATGGAGAAGGGGCGTTGCAGCGCTTCATCAACTGCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCg                    <  1:6552918/92‑1 (MQ=255)

                                                                    aCCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCg                    >  2:7781413/1‑93 (MQ=255)

                                                                    aCCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCg                    <  2:3967677/93‑1 (MQ=255)

                                                                         cGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCaa               >  1:6537418/1‑93 (MQ=255)

                                                                          gATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCaaa              >  1:6259362/1‑93 (MQ=255)

                                                                            tGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGt            >  1:7817735/1‑93 (MQ=255)

                                                                             ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt           <  1:6989489/93‑1 (MQ=255)

                                                                             ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt           <  1:6162696/93‑1 (MQ=255)

                                                                             ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt           <  2:3310301/93‑1 (MQ=255)

                                                                              gagaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGttt          >  2:1193018/1‑93 (MQ=255)

                                                                               agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTtg         <  2:6513164/93‑1 (MQ=255)

                                                                               agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTtg         <  2:6993379/93‑1 (MQ=255)

                                                                                gaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTtgt        <  2:4910207/93‑1 (MQ=255)

                                                                                   ggCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTc     <  1:2414636/93‑1 (MQ=255)

                                                                                   ggCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTc     >  2:947924/1‑93 (MQ=255)

                                                                                    gcgcgcTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCg    <  2:5276807/93‑1 (MQ=255)

                                                                                     cgcgCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgt   <  1:5883992/93‑1 (MQ=255)

TGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTG  >  NC_000913/908648‑908826

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